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SUMMARY 


This  report  describes  the  results  of  a  five-year  research  program  to  provide  the  basic 
information  needed  to  design  and  develop  high  performance  lasers  systems,  operating  in  the  visible 
and  near  ultraviolet,  for  use  in  blue-green  optical  communication  systems  and  in  weapon  systems. 
All  such  systems  require  a  reliable,  scalable,  and  efficient  laser.  In  several  applications,  operation 
within  a  specific  wavelength  range  is  required,  along  with  narrow  spectral  output  to  improve 
optical  conversion  efficiency  or  to  match  fixed-frequency  narrow-band  receivers. 

Over  the  past  twenty  years,  the  search  for  and  development  of  efficient,  high  energy, 
visible  and  ultraviolet  lasers  have  been  extremely  productive.  This  phenomenally  rapid  progress 
was  made  possible  by  exploiting  the  previously  established  base  of  molecular  concepts  and 
experimental  data.  Detailed  understanding  of  the  kinetic  and  optical  processes  operating  within  the 
laser  medium  is  needed  to  identify  new  molecules,  pumping  mechanisms,  or  optical  conversion 
schemes  that  could  lead  to  high  performance  lasers  and  to  guide  the  achievement  of  the  best 
possible  performance.  Sustaining  the  accelerated  pace  of  development  requires  a  continuously 
improving  description  of  the  microscopic  processes  within  the  laser  medium,  which  in  turn 
requires  continuing  research  programs  to  investigate  these  processes. 

During  the  period  covered  by  this  contract,  our  research  program  emphasized  the 
quantitative  characterization  of  optical  and  kinetic  processes  important  for  XeCl  and  XeF  excimer 
laser  systems.  The  documents  included  in  this  report  as  Appendices  describe  our  major  research 
accomplishments.  These  may  be  divided  into  four  areas,  as  briefly  outlined  below. 

EXCIMER  LASER  SPECTROSCOPY  AND  KINETICS 

Excimer  lasers  are  among  the  most  complicated  of  gas  lasers  because  many  chemical 
reactions  are  involved  in  production,  relaxation,  and  quenching  of  the  upper  and  lower  laser  levels. 
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Nevertheless,  rare  gas  halide  excimer  lasers  are  remarkably  efficient  The  complicaied  chemical 
pathways  are,  in  fact,  very  selective. 

The  kinetics  of  the  XeCl  laser  are  especially  worthy  of  investigation.  Even  after  years  of 
study,  significant  inconsistencies  and  contradictions  remain.  We  have  used  vacuum  ultraviolet 
synchrotron  radiation  from  the  Spear  storage  ring  at  the  Stanford  Synchrotron  Radiation 
Laboratory  to  excite  Xe  atoms  and  halogen-containing  molecules.  The  resulting  fluorescence  was 
time  resolved  and  analyzed  to  derive  rate  coefficients  for  production  and  quenching  of  XeCl  and 
Xe2Cl.  These  experiments  provided  branching  ratios  for  production  of  XeCl  from  various  Xe* 
states  reacting  with  a  number  of  chlorine  donors.  These  studies  show  that  pathways  that  might 
have  been  expected  to  dominate,  e.g.,  Xe*(6s)  +  HCl,  in  fact  have  a  low  probability  for  producing 
XeCl*.  These  results  strongly  support  the  need  for  developing  a  consistent  understanding  of  the 
role  of  nKxlerate  energy  electrons  in  mixing  the  excited  states  of  the  Xe  atom  and  in  exciting  the 
higher  vibrational  levels  of  HQ. 

The  XeF  laser  is  unusual  among  tiie  rare  gas  halide  lasers  in  that  the  ground  state  of  XeF  is 
bound.  Nonterminating  laser  action  requires  that  the  lower  laser  level  be  removed,  in  this  case  by 
collisional  vibrational  relaxation  and  dissociation.  We  measured  the  rates  of  collision-induced 
dissociation  of  XeF(X)  by  He,  Ne,  Ar,  Kr,  Xe,  N2,  SF6,  and  XeF2,  and  of  vibrational  relaxation 
by  XeF2.  In  our  experiment  XeF(Xv")  was  produced  by  KrF  or  ArF  excimer  laser 
photodissociation  of  XeF2  and  then  excited  by  a  doubled  dye  laser  to  XeF(Bv')>  whose 
fluorescence  was  detected.  The  time-evolution  of  the  XeFfXv")  population  was  monitored  by 
varying  the  time  delay  between  the  excimer  and  dye  lasers. 

Our  work  on  excimCT  laser  spectroscopy  and  kinetics  is  described  in  Appendices  A 
through  G. 

STIMULATED  VIBRATIONAL  AND  ROTATIONAL  RAMAN  SCATTERING 

Stimulated  Raman  scattering  (SRS)  is  a  convenient  means  of  frequency  shifting  of 
ultraviolet  excimer  lasers  to  permit  transmission  through  the  atmosphere.  Stimulated  Raman 
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scattering  is  also  useful  for  improving  the  optical  beam  quality  of  excimer  lasers  by  combining 
several  excimer  laser  beams  at  small  angles  in  a  Raman  conversion  cell.  Finally,  the  atmosphere 
itself  can  lead  to  stimulated  Raman  scattering,  which  can  limit  the  power  density  of  lasers  beams 
propagating  through  the  atmosphere. 

Our  research  program  was  intended  to  obtain  the  detailed  quantitative  experimental  data 
required  to  design  practical  frequency  converters  and  beam  combiners  and  to  quantify  the  effects  of 
Raman  scattering  in  the  atmosphere.  Our  work  on  vibrational  Raman  scattering  in  H2  included 
measurement  and  modeling  of  the  angular  dependence  of  the  Raman  linewidth,  measurement  and 
ab  initio  calculation  of  the  Raman  scattering  cross  sections  of  the  fundamental  and  overtone 
transitions,  and  determination  of  the  refractive  index  or  polarizibility  of  H2  in  its  first  vibrationally 
excited  level.  We  investigated  rotational  Raman  scattering  in  considerable  detail.  Our  nearly 
complete  determination  of  the  Raman  scattering  parameters  allowed  us  to  develop  a  comprehensive 
model  of  Raman  gain  coefficients  for  N2  in  the  atmosphere.  This  work  had  a  substantial  influence 
on  development  of  strategies  for  using  power  lasers  in  the  atmosphere. 

Our  work  on  stimulated  vibrational  and  rotational  Raman  scattering  is  described  in 
Appendices  H  through  M. 

STIMULATED  RAYLEIGH-BRILLOUIN  SCATTERING 

Stimulated  Brillouin  scattering  (SBS)  is  used  in  optical  phase-conjugation  mirrors  to 
compensate  for  refractive  distortions  in  the  excimer  laser  medium.  We  developed  a  high-resolution 
(10-MHz)  coherent  Brillouin  spectrometer  that,  for  the  first  time,  clearly  resolves  the  Rayleigh  and 
BriUouin  components  of  the  low  frequency  scattering  spectrum  of  gases.  We  measured  the  SBS 
gain  coefficients,  linewidths,  and  line  shifts  for  Xe  and  SFe,  improved  the  accuracy  of  previous 
measurements,  and  illustrated  the  power  of  the  technique  for  characterizing  potential  SBS  media. 
Our  work  is  described  in  Appendices  N  and  O. 
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OTHER  TOPICS  IN  LASERS  AND  NONLINEAR  OPTICS 


The  SRS  and  SBS  measurements  described  above  require  lasers  with  a  combination  of 
narrow  spectral  width  and  relatively  high  power.  To  perform  such  measurements,  we  developed  a 
Nd:YAG  laser  with  exceptional  performance  characteristics;  80-ns  pulse  duration,  8-MHz 
bandwidth,  300-mJ  pulse  energy  and  10-Hz  pulse  repetition  frequency.  The  laser  system  has 
formed  the  basis  of  numerous  state-of-the-art  experiments  in  nonlinear  optics.  It  is  described  in 
Appendix  P. 

The  high  efficiency  of  the  KrF  laser  has  led  to  numerous  schemes  to  convert  its  wavelength 
from  248  nm  to  the  visible,  which  would  allow  it  to  propagate  through  the  atmosphere.  Such  a 
conversion  of  a  248-nm  photon  to  a  400  -  600-nm  photon  results  in  a  loss  of  at  least  40%-60%  of 
the  photon's  energy.  We  have  investigated  a  new  scheme  by  which  a  single  KrF  photon  would  be 
converted  into  two  visible  photons,  by  four-wave  mixing  with  a  10-^m  photon  from  a  CO2  laser: 

248  nm  +  10  pm  484  nm  +  484  nm . 

This  process,  which  we  called  frequency  halving,  could  accomplish  the  conversion  to  the  visible 
with,  in  principle,  no  loss  of  energy.  Generation  of  two  waves  at  the  same  frequency  is  known  to 
occur  in  a  number  of  parametric  processes  in  solids.  Appendix  Q  describes  our  first  attempt  at 
frequency  halving-in  barium  vapor. 
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Appendix  A 


EXCITED  STATE  KINETICS  FOR  XeCI 


International  Conference  on  Lasers  '84 
26-30  November  1984,  San  Francisco,  CA 


reprinted  from  proceedings  of  the  international  conference  on  lasers  ’84.  NOVEMBER  26-30,  1984 


EXCITED  STATE  KINETICS  FOR  XeCI^ 


Donald  C.  Lorents 
Chemical  Physics  Laboratory 
SRI  International 
Menlo  Park,  California 


Although  the  XeCl  laser  has  operated  with  e-beam  and  discharge  pumping  for  nearly  10 
years  now  several  Intrlqulng  mysteries  remain  In  both  the  Ion  and  excited  state  kinetic 
behavior  of  Its  medium.  In  Che  positive  Ion  chemistry  the  pathways  to  the  formation  of 
xenon  Ions  In  a  neon  buffer  are  complex  and  111  defined.  In  mixtures  containing  HCl  Che 
negative  Ion  formation  processes  are  still  confronted  with  apparently  contradictory 
experimental  results. 

Although  the  excited  state  kinetics  of  XeCl  are  understood  on  a  qualitative  basis, 
many  of  the  Important  kinetic  rate  constants  remain  uncertain  or  unmeasured.  As  we  look  In 
more  detail  at  the  various  kinetic  pathways  more  questions  arise  In  the  overall  model  of 
XeCl  formation,  radiation,  and  quenching  particularly  as  it  pertains  Co  laser  performance 
and  modeling.  Such  issues  as  the  colllslonal  mixing  and  vibrational  relaxation  of  Che  B 
and  C  states  remain  highly  uncertain  In  spite  of  their  Importance  to  Che  understanding  of 
laser  performance.  Indeed  the  spectroscopic  location  of  Che  C  state  relative  to  the  B 
state  remains  uncertain  by  about  200  cm  Evidence  of  'the  poor  understanding  of  the  XeCl 
laser  medium  Is  particularly  clear  for  the  case  where  HCl  is  used  as  the  Cl  donor.  In  e- 
beam  pumped  XeCl  a  vexing  question  Is  why  Che  laser  turns  off  before  Che  pump  pulse  does 
while  Che  fluorescence  continues.  Indicating  chat  the  cause  is  not  depletion  of  Che 
donor.  The  laser  behavior  at  low  HCl  concentrations  giving  high  initial  peak  power  and 
rapid  termination  Is  not  understood.  These  problems  are  frequently  assigned  to  unidenti¬ 
fied  absorbers  formed  through  some  unknown  kinetic  reaction  scheme.  These  characteristics 
pose  difficult  problems  which  the  current  models  find  difficult  to  reproduce. 

At  LASERS  '80  Che  first  results  of  our  synchrotron  studies  of  the  kinetics  of  XeCl 
were  presented.  In  Che  Intervening  time  we  have  extended,  and  improved  the  measurements, 
Che  data  analysis,  and  ability  to  extract  rate  constants  based  on  an  Improved  understanding 
of  Che  overall  excited  state  kinetics  of  XeCl. 

We  present  here  results  of  our  extended  sec  of  excited  state  kinetics  measurements 
using  Che  tunable  uv  synchrotron  source  at  SSRL  which  has  Ideal  temporal  and  Intensity 
characteristics  for  such  studies.  These  experiments  provide  a  means  for  obtaining  In  a 
simple  reliable  and  readily  Interpretable  manner  many  of  the  important  kinetic  para- 
■eters.  Very  briefly,  uv  photons  of  specified  energy  are  used  to  excite  an  atom  or 
Bolecule  Co  a  desired  electronic  state  In  a  gas  mixture  containing  the  desired  reactants 
snd  or  quenchers.  The  formation  and  decay  of  a  fluorescing  exclmer  Is  followed  In  time  and 
Intensity  as  a  function  of  Che  concentration  of  Individual  components  of  the  gas  as  well  as 
Che  total  pressure.  As  far  as  possible  an  Individual  reaction  is  isolated  by  arranging  the 
gas  mixtures  and  excitation  conditions  to  make  It  the  dominant  process  in  which  case  Its 
reaction  rate  can  be  determined  uniquely.  When  chat  proves  Impossible  It  becomes  necessary 
to  model  Che  several  processes  Involved  using  Initial  estimated  rates  and  a  large  set  mea¬ 
surements  taken  over  a  large  range  of  gas  mixtures  and  pressures. 

These  measurements  were  done  on  the  8  degree  port  of  the  SSRL  facility  which  delivers 
0.3  ns  pulses  at  1.28  MHz  tunable  from  40  to  300  nm.  The  gas  cell  Is  equlped  with  entrance 
and  exit  windows  Co  transmit  the  uv  radiation  and  a  third  window  looking  perpendicular  to 
Che  beam  to  monitor  the  fluorescence.  Filtered  photomultipliers  were  used  to  detect  the 
fluorescent  photons.  MgFl  windows  used  on  the  gas  cell  limit  the  short  wavelength  trans¬ 
mission  to  115  nm.  Gas  mixtures  are  made  up  In  advance  of  a  measurement  using  a  ballast 
tank  and  a  well  passivated  gas  mixing  apparatus.  Gas  pressures  are  measured  with  cali¬ 
brated  capacitance  manometers. 

An  absorption  scan  of  a  He/Xe  mixture  at  1027/1  Torr  Is  shown  In  Figure  1  for  the 
wavelength  range  of  115-150  nm  taken  with  a  resolution  of  0.25  nm.  The  intensity  scale  is 
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proporclonal  to  the  absorption  cross  sec¬ 
tion  having  been  normalized  to  the  beam 
Intensity  profile  obtained  from  a  scan  at 
zero  pressure.  All  of  the  absorption  lines 
observed  are  Identifiable  Xe  lines  colli¬ 
sion  broadened  by  the  high  density  He.  It 
Is  noteworthy  that  not  only  are  the  allowed 
transitions  broadened  but  also  some  transi¬ 
tions  normally  disallowed  in  atomic  Xenon 
are  observed  to  absorb  In  the  presence  of 
high  density  buffer  (e.g.  ^*^1/2 

^‘*5/2^'  ^  similar  scan  for  an  Ar/Xe  mix  at 

1021/1  Torr  shows  that  the  broadened  struc¬ 
ture  of  the  lines  differs  from  those 
observed  In  He  but  otherwise  Is  similar. 
These  collision  induced  absorption  features 
have  been  observed  and  analyzed  In  higher 
resolution  by  Castex.^ 


Figure  1  Absorption  scan  of  1  Torr  Ze  In 
1027  Torr  He. 


In  Figure  2,  0.2  torr  of  CI2  has  been 
added  to  an  Ar/Xe  mix  of  1000/1  Torr  pro¬ 
ducing  a  very  complex  absorption  spectrum. 

In  addition  to  the  Xe  lines,  strong  absorp¬ 
tion  features  are  observed  that  have  been 
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Identified  in  higher  resolution  studies  as 
transitions  to  various  Rydberg  states  of 
CI2.  In  this  measurement  the  fluorescence 

at  308  nm  was  monitored  and  the  intensity  was  normalized  to  the  absorption  spectrum  at  the 
Xe(7s2^2^  state.  It  is  clear  that  the  major  XeCI  308  nm  fluorescence  is  due  to  excitation 
of  the  Xe  states  and  very  little  comes  from  the  CI2  excitation.  However  it  is  important  to 
note  that  there  Is  some  fluorescence  signal  at  136  nm  where  the  major  CI2  absorption  lies. 
This  signal  Increases  as  the  Xe^  concentration  increases  and  Is  clearly  due  to  the  reaction 
of  CI2  with  Xe  CO  produce  XeCl  .  The  chlorine  excitation  at  136  nm  Is  rapidly  quenched  by 
Ar  to  lower  states  so  Chat  this  mechanism  of  XeCl  formation  Is  klnecically  fast  and  thus  a 
very  useful  means  of  producing  It  for  race  measurements  on  ocher  kinetic  reactions.  In 
Figure  3  the  XeCl  308  signal  Is  plotted  together  with  an  absorption  scan  of  Ar/Xe  containing 
no  CI2  and  shows  that  Che  Xe  excited  states  give  high  but  not  lOOZ  yields  for  Che  lower 
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Figure  2  AbsorpClou/Fluorescence  scan  of 
an  Ar/Xc/Clj  mix. 
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Figure  3  Fluorescence  data  of  Figure  2 
normalised  to  an  absorption  scan  in  Ar/Ze. 
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_  lying  states.  The  yield  for  the 

state  which  is  53  ±  52  is  surpislng  in  view 
of  the  1002  yield  reported^  for  the  nearby 
but  lower  lying  ^^3/2  metastable 

level  (J  “  2).  In  Figure  4  an  absorp- 

tlon/f luorescence  scan  is  shown  for  a  mixt¬ 
ure  of  He/Xe/HCl  of  1000/3/1  Torr.  This 
shows  that  HCl  which  is  commonly  used  for 
the  Cl  donor  in  XeCl  lasers  gives  good 
fluorescence  yields  only  for  the  5d  and  7s 
and  presumably  higher  states  lying  above 
about  10  eV.  At  this  energy  the  reaction 
is  about  1.5  eV  exothermic  which  seems  to 
be  a  typical  requirement  for  a  high  yield 
harpoon  reaction.  This  reaction  is  thermo- 
neutral  for  Xe(6s2^2^  where  the  yield  is 

1150  1220  1290  1360  1430  1500 

WAVELENGTH  (A)  Flourescence  yields  obtained  by  the 

JA-7123.JI 

technique  described  above  for  several 
fluorine  and  chlorine  donors  are  shown  in 
Table  1.  The  quanitltles  in  parentheses 
indicate  where  the  yields  have  been 
Figure  4  AbsorpCloa/Fluoreaceace  scan  of  normalized.  In  the  cases  of  CIF  and  HCl 

a  He/Xe/HCl  mix.  XeCl  yields  were  normalized  to  the 

Xe(7s2y2^  yield  in  Cl^.  It  is  Interesting 
that  CIF  reacting  with  Xe  has  a  strong 
preference  to  form  XeCl  relative  to  XeF 

indicating  that  the  harpoon  reaction  is  strongly  Influenced  by  Che  higher  electron  affinity 
of  the  Cl  atom. 

Kinetic  information  is  obtained  from  measurements  of  the  time  history  of  the 

fluorescence  signal  as  a  function  of  total  and  component  gas  pressures  for  each  gas 
mixture.  To  obtain  unique  kinetic  information  on  a  specific  reaction  it  is  necessary  to 


Figure  4  AbsorpCion/Fluorescence  scan  of 
a  He/Xe/HCl  mix. 


Table  1 


QUANTUM  YIELDS  FOR  Xe*  +  RX  XeX*  +  R 


XeF 

XeCl 

* 

Xe 

^2 

NF3 

HU 

CIF 

CIF 

HCl 

CI2 

683/2 

0.45 

0.05 

0.21 

0.02 

0.29 

0.0 

0.53 

6si/2 

0.75 

0.50 

0.58 

0.04 

0.46 

0.05 

0.95 

5dj/2 

1  .0 

0.70 

1  .0 

? 

0.60 

0.30 

1  .00 

5d5/2 

1  .0 

1  .0 

1 .0 

0.05 

0.60 

0.7 

1.00 

0.90 

0.90 

0.90 

0.04 

0.50 

0.65 

0.95 

’'*3/2 

(1.0) 

(1.0) 

(1.0) 

0.05 

0.60 

0.72 

(1.0) 
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Isolate  that  reaction  as  much  as  possible  by  arranging  the  gas  mixture  and  pressures  such  as 
to  make  It  the  dominant  process  controlling  either  the  fast  or  slow  component  of  the 
signal.  It  Is  also  Important  that  the  rates  of  the  fast  and  slow  components  are  well 
separated  If  a  realistic  exponential  fit  Is  to  be  obtained  from  a  least  squares  fitting 
routine . 

In  Figure  5  Is  shown  the  results  of  a  measurement  of  the  rate  of  formation  of  XeCl* 
from  Xe(5dj^2)  reacting  with  HCl.  This  rate  Is  found  to  be  dependent  on  both  the  HCl  and  Xe 
concentrations.  Evidently  the  co Ills lonally  deactivated  by  Xe  to  a  lower 

state  that  does  not  prodttce  XeCl  giving  a  Xe  pressure  dependence  which  appears  very  clearly 
by  using  two  different  Xe/HCl  ratios.  The  decay  frequencies  as  a  function  of  HCl  pressure 
shown  for  100/1  and  10/1  Xe/HCl  ratios  are  observed  to  be  linearly  dependent  on  pressure 
within  the  error  limits.  Thus  we  can  readily  assign  the  slope  of  these  as  the  sum  of  the 
HCl  and  Xe  reaction  rates  and  using  the  two  equations,  derive  the  two  rate  constants  shown 
on  the  figure. 

To  obtain  the  XeCl  quenching  rate  due 

to  HCl  we  used  a  constant  mixture  of  Xe/Cl2 

to  which  HCl  was  added.  In  this  case  the 

XeCl  was  formed  from  the  Cl^  reaction  with  v 

*  *  ^ 

Xe  In  order  to  produce  the  XeCl  fast  com-  ^ 

pared  to  Its  radiative  and  quenched  life¬ 
time.  The  decay  rates  as  a  function  of  HCl 
are  displayed  In  Figure  6  for  Xe/Cl2  pres¬ 
sures  of  200/0.3  and  SO/0.3  and  In  each 
case  a  linear  dependence  on  HCl  Is  observed 
with  a  common  slope  but  differing 
asymtotes.  An  accurate  HCl  quenching  rate  Hi 

constant  Is  easily  derived  from  Che  slope 
of  the  two  straight  line  least  squares 
fits.  In  addition  the  two  asymtoclc  values 
give  Che  XeCl  decay  race  for  the  respec¬ 
tive  Xe/Cl2  mixtures  Indicating  a  strong 
dependence  of  this  rate  on  Che  Xe  concen¬ 
tration. 

Figure  5  Formation  rate  of  XeCl  from 
Xe(Sd3^2)  BCL  and  removal  rate  of 

Xe(5d3^2)  ** • 


JA-71J3-33 


Figure  6  Quenching  rate  of  XeCl  by  HCl. 


HO.  PRESSURE  aORR) 


Similar  measurements  using  a  450  nm 
filter  were  carried  out  to  obtain  the  life¬ 
time  and  quenching  of  Xe2Cl.  In  this  case 
it  was  found  that  the  decay  rate  of  Xe2Cl 
was  Independent  of  Che  Xe  pressure  and  only 
dependent  on  the  CI2  and  HCl  concentrations 
permitting  a  direct  measure  of  the  halogen 
quenching.  The  measurements  shown  In 
Figure  7  carried  out  In  Xe/Cl2  mixtures  at 
various  Xe  pressures  from  50  to  2000  torr 
gives  Che  CI2  quenching  rate.  The 

assymCoCe  gives  the  radiative  decay  rate  of 
Xe2Cl  of  3.8  X  10^  s  corresponding  to  a 
lifetime  of  about  260  ns,  a  value  longer 
than  any  previously  reported.  The  HCl 
quenching  rate  was  derived  from  a  measure¬ 
ment  of  Che  450  nm  decay  in  a  Xe/Cl2  mix¬ 
ture  of212/0.3  Torr  co  which  HCl  was  added. 

Finally  Che  decay  and  quenching  of 
XeCl  itself  Is  more  complex  because  of  the 
Involvement  of  the  two  close  lying 
electronic  B  and  C  states.  Since  the  B  and 
C  states  of  XeCl  are  known  co  lie  within 
~  200  cm~*'  of  each  ocher  they  are  strongly 
coupled  by  collisions  and  thus  at  most 
pressures  act  as  a  single  state  with  an 
Intermediate  lifetime.  Our  measurements 
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reflect  this  by  the  observation  that  the 
C-A  and  B-X  transitions  at  340  no  and 
308  no  respectively  have  Identical  decay 
rates.  The  radiative  lifetimes  of  each  of 
these  states  have  been  recently  measured 

4 

under  collision  free  conditions  but  the 
mixed  state  lifetime  and  the  equilibrium 
constant  are  not  yet  well  established  for 
this  exclmer;  the  reported  energy  separa¬ 
tions  vary  from  0  to  -200  cm  Because 

the  mixed  state  Is  quenched  by  both  Xe  and 
Cl^,  and  the  dependence  on  Xe  could  be  both 
linear  and  quadratic  corresponding  to 
quenching  and  to  Xe201  formation  respec¬ 
tively,  It  was  necessary  to  obtain  a  large 
data  set  with  a  large  range  of  pressures 
and  mixture  ratios.  The  decay  rate,  R,  was 
assumed  to  obey  a  four  parameter  function  R 
given  as  equation  (1). 

R  -  Kj  +  K2CI2  +  KjXe  +  K^Xe^  (1) 

Some  47  sets  of  data  were  fit  to  R  with  a 
multi-parameter  least  squares  fitting  rou¬ 
tine  to  determine  the  K^'s.  This  fit  for  the  Xe/Cl2  case  Indicated  that  the  linear  Xe  term 
Is  small  but  chat  the  three-body  quencjilng  Is  predominant  even  at  low  pressures.  This  Is 
Illustrated  In  Figure  8  where  the  XeCl  decay  frequency  measurements  corrected  for  the  CI2 
quenching  rates  are  plotted  versus  the  square  of  the  Xe  pressure.  The  plot  clearly  exhibits 
a  linear  dependence  that  Indicates  chat  Che  two-body  Xe  quenching  Is  very  small  even  at  low 
pressures.  Additional  dace  not  shown,  extending  a  factor  of  6  high  in  (pressure)^  fits  the 
same  function.  Nevertheless  Che  least  squares  fitting  does  give  a  small  positive  value  for 
Che  linear  Xe  term  with  a  large  (factor  of  2)  uncertainty.  The  value  we  quote  is  the  upper 
limit  Chat  Che  statistical  analysis  on  our  data  allows. 

Finally  an  additional  large  set  of  data  for  mixtures  of  Ar/Xe/Cl2  was  combined  with  the 
Xe/Cl2  data  sec  and  a  least  squares  fitting  of  Che  total  set  was  performed  to  an  R  function 
with  an  additional  term  in  Ar  Xe .  The  resulting  fit  gave  essentially  the  same  parameter 
values  as  found  for  Che  Xe  data  with  Che  addition  of  the  term  for  Ar  Xe .  The  fit  for  this 
total  data  sec  of  87  points  is  excellent  and  clearly  indicates  chat  Che  overall  model 
describes  Che  decay  processes  of  XeCl  very 
well.  This  model  chat  prescribes  slow 
2-body  and  fast  3-body  removal  of  XeCl  by 
Xe  differs  from  the  recent  measurements  of 
the  Kansas  group^  who  determined  a  2-body 
race  a  factor  of  ten  larger  chan  our  mea¬ 
surements  indicate.  However  a  recent  sec 
of  measurements  carried  out  by  a  European 
group^  using  synchrotron  radiation  In  an 
experiment  very  similar  to  ours  has  pro¬ 
duced  Xe  and  CI2  removal  races  nearly 
identical  to  chose  reported  here. 

An  additional  important  parameter  chat 
results  from  these  measurements  is  Che 
mixed  state  radiative  race  of  XeCl  .  This 
quantity  is  given  by  Che  pressure  indepen¬ 
dent  term  in  the  least  squares  fitting  and 
is  interpreted  as  the  mixed  state  life¬ 
time.  This  lifetime  obtained  from  Che 
fitting  of  the  total  data  set  described 
above  taken  together  with  the  individual  B 
and  C  state  lifetimes^  can  be  used  to 
determine  energy  separation  of  Che  two 


Figure  8  Three  body  removal  rate  of  ZeCl* 
by  Xe.  Rote  that  the  CI2  quenching  rate 
haa  been  eubstracted  from  the  total  decay 
rate. 
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Figure  7  Radiative  and  quenching  rate  of 
Xa2^^  by  HCl  and  Cl2. 
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electronic  states.  The  mixed  state  radiative  rate  A„  is  given  in  terms  of  the  B  and  C  state 
rates  and  the  equilibrium  constant  in  equation  (2) 


+  Ac)/(Ke  +  1) 


Using  the  recently  measured  rates  for  Ag  and  A^,  and  our  value  for  A^^^  of  5.0  x  10  /s  we  find 
that  K  •  ~  1  indicating  an  energy  separation  near  sero.  This  finding  is  consistent  with 
the  recent  work  of  the  European  group'  but  differs  from  other  observations. 

The  rate  constants  and  lifetimes  derived  from  the  present  analysis  of  our  synchrotron 
data  are  presented  in  the  Table  II.  In  these  table  KSU  refers  to  Kansas  State  University, 
SAC  refers  to  reference,^  and  R.ICE  to  reference  8. 


Table  II 

RATE  CONSTANTS  FOR  XeCl  KINETICS 
FORMATION 


REACTION 

Xef^Pj)  +  Clj  -►  XeCl* 

Xe(^Pj)  +  Clj  5CeCl^ 

XeC^Pp  +  Clj  -►  XeCl  ^ 

Xe(5d3/2)  +  HCl  *  XeCl 

CI2*  *  Xe  ♦  XeCl 


RATE  CONSTANT  REF 

7.2  X  10"^°  CM^/S  KSU 
7.9  +  0.9  X  10~^°CM^/S  SRI 
7.6  +  0.7  X  10~^°CM^/S  SRI 
7.5  ±  0.7  X  10~^°CM^/S  SRI 

1.1  +  0.1  X  10~^°CM^/S  SRI 


RADIATION  AND  QUENCHING 


XeCl*(B)  ■*  hv(308) 
XeCl*(C)  -►  hv(340) 
XeCl‘(M)  ♦  hv 


XeCl* 

+ 

CI2  ♦  QUENCH 

XeCl* 

+ 

HCl  *  QUENCH 

XeCl* 

+ 

Xe  QUENCH 

XeCl* 

+ 

2Xe  ♦  Xe2Cl* 

XeCl* 

+ 

Xe  +  Ar  ♦  Xe2Cl 

9.0  + 

0.2  X  lO^/S 

KSU 

7.6  ± 

0.6  X  10^/S 

KSU 

5.0  + 

0.6xl0'/S 

SRI, SAC 

5.6  ± 

0.25  X  lO'^'^CM-’/S 

SRI 

5.8  X 

10"^°  CM^/S 

SAC 

7.3  ± 

0.1  X  10“*°CM^/S 

SRI 

<  4  X 

10"^^CM^/S 

SRI, SAC 

2.3  X 

10"^^CM^/S 

KSU 

1.53 

i  0.1  X  10"^°CM®/S 

SRI 

1.3  X 

10"^°  CM^/S 

SAC 

1 .01 

+  0.05  X  10"^°CM®/S 

SRI 

Xe  QUENCHING 


Xe(5d3/2)  +  Xe  -»  Xe  +  Xe 
XeC^Pj)  +  HCl  PRODUCTS 

Xe(^Pl)  +  HCl  PRODUCTS 

Xe,*  +  HCl  ♦  PRODUCTS 
Xej  +  CI2  ♦  PRODUCTS 


8.1  +  0.8  X  10~^^CM^/S  SRI 
5.6  X  10"1°CM^/S  KSU 
(7  X  10"^°CM^/S)  SRI 

8.2  +  0.8  X  10"‘°CM^/S  SRI 
5.0  +  0.6  X  10~^°CM^/S  SRI 


XejCl*  RADIATION  & 

XejCl*  ♦  hv(470nm) 

Xe2Cl*  +  CI2  -►  PRODUCTS 
Xe2Cl*  +  HCl  -►  PRODUCTS 
Xe2Cl*  +  Xe  PRODUCTS 


QUENCHING 


3.8  +  0.4  X  10^/S  SRI 

7.4  X  lO^/S  RICE 

5.2  +  0.2  X  10"‘°CM^/S  SRI 

6.1  ±  0.2  X  10~^°CM^/S  SRI 

<  4  X  10~^^CM^/S  SRI 
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Appendix  B 


THE  ROLE  OF  ELECTRONS  IN  EXCIMER  LASERS 


International  Conference  on  Lasers  *85 
26  December  1985,  Las  Vegas,  NV 
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MANY  GAS  LASERS  ARE  ELECTRON-EXCITATION  OOMINATEO 
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needed*  and  usually  only  for  electron/ground-state 


EXCIMER  LASERS  ARE  DOMINATED  BY  HIGH-PRESSURE  CHEMICAL  KINETICS 
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Figure  2.  ' - •  ‘ - ' 

THIS  ILLUSTRATION  DEPICTS  THE  MAJOR  ENERGY  FLOW  PATHWAYS  IN 
E-BEAM  PUMPED  Ne/XeINF,  MIXTURES.  ENERGY  INPUT  FROM  THE 
ELECTRONS  COMES  IN  ALONG  THE  PUMP  ARROWS.  THE  SPECIE  IN  EACH 
BOX.  PLUS  THE  COMPONENT  ALONG  THE  REACTION  ARROW.  YIELDS  THE 
RESULT  IN  THE  NEXT  BOX.  THE  WAVELENGTHS  OF  THE  EMITTING  SPECIES 
ARE  INDICATED.  THEY  MAY  BE  SUBJECT  TO  QUENCHING  BY  NF,  OR  Xe  (NOT 
SHOWN).  NeF*  IS  PRESUMED  TO  PREDISSOCIATE. 
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RADIATIVE  LIFETIMES  OF  THE  O’  SUBLEVELS 
OF  THE  RARE  GAS  EXCIMERS 

David  L.  Hues t is 
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SRI  International 
Menlo  Park,  CA  94025 
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INTRODUCTION 


The  prototype  "excimer"  molecules. 

Still  interesting  as  efficient  light  sources, 
and  still  incompletely  tinders tood. 

The  weakly  bound  ground  state  complicates 
conventional  spectroscopic  investigation. 

In  the  present  study  we  are  investigating 
theoretically  the  fine  structure  of  lowest  of 
the  excimer  excited  states, 

We  are  especially  interested  in  the  radiative 
properties  of  the  0^  sublevel,  which  previous 
studies  have  supposed  to  be  completely 
metastable . 
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1900A 


A 


ORIGIN  OF  TRANSITION  PROBABILITY 

X  -  X 

0^  weak  V21  and  R2^  branches 
If  -♦  0^  weak  Q3^  branch 
Of  0^  very  weak  branch,  -  J(J+1) 

and  If  mix  with  by  spin- orbit  coupling 
Of  mixes  with  If  by  Coriolis  coupling 


C-4 


PERTURBATION  TREATMENT  IN  HUND'S  (a) 


A  -  E(^Sj)  -  E(^So)  is  large 

IA3  “  l/r2  ~  lAl  -  IPqI^ 

Pi  -  <^Si|p|%>  -  <^S|Hg^|%  <%|/»|  W(^n-3s) 
PO  “  <%|Pl^S0>  “  <^Sol«rotl^Si> 

<3rol“rotl^2l>  “  -2B.7[J(J+1)]*‘ 

«7  -  -  J2  ? 


raAi-  2(BVA)^  J(J+1) 
<r2/ri>  -  ~  ABkT/A^ 
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Rg2  CORIOLIS  COUPLING  IN  HUND'S  CASE  (d) 
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CALCtJLATED  Rg?  ENERGIES  (cm‘^) 
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RADIATIVE  LIFETIMES  (/is) 
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0.24 

<r> 

20 

9.4 
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0.35 

<expt> 
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3.2 
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0.099 
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2.4 
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NUCLEAR  SEPARATION  (A) 


FIGURE  1  XENON  INTERMOLECULAR  POTENTIAL  CURVES 
FROM  THE  WORK  OF  MULLIKAN 


CONCLUSIONS 


The  0^  sublevels  of  the  rare  gas  exclmers  have 
non-negligible  radiative  transitions  induced  by 
rotational  coupling. 

In  high  pressure  rare  gases  the  sublevels  of 
the  state  should  be  well  mixed.  The 
observed  radiative  rates  should  be  averages  of 
the  rates  from  the  1^  and  0^  sublevels. 

There  remains  a  factor  of  2-3  disagreement 
between  theory  and  experiment. 

The  lifetimes  of  the  0^  sublevel  in  Kr^  and  Xe^ 
are  in  a  favorable  range  for  molecular  beam 
investigation . 
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Three-body  Reaction  of  XeCl*  to  Form  Xe2Cl,* 

D.  LORI' NTS,  R.  L.  SHARPLESS,  D.  L.  HUESTIS,  SRI  INTER- 
NATIONAL  —  Synchrotron  mensurcments  of  the  decay  rate  of 

'  '  ~  ^  ‘  '  it 

XeCl  and  the  formation  rate  of  Xe2Cl  In  Ne/Xe/Cl2  and 

Kr/Xe/Cl2  mixtures  is  romposod  of  three  components: 

radiation,  two-body  quenching  and  three-body  reaction. 

The  measurements  at  pressures  up  to  5000  Torr  show  that 

the  reaction,  XeCl*  +  Xe  +  Ne  +  Xe2Cl  +  Ne,  is  not  a 

Linear  function  of  the  [XellNel  product  and  therefore  is 

not  a  simple  three  body  reaction.  Measurements  at  Xe 

pressures  of  50  and  150  Torr  and  Ne  pressures  ranging 

from  0  to  5000  Torr  can  be  fit  with  the  simple  pressure 

-3 1 

dependent  three-body  rate  constant  given  by  kp  =  7x10 
[Xe](Ne]/|l  +  10~^*^[Ne)[.  This  result  can  be  quali¬ 
tatively  understood  In  terms  of  the  Undemann  mechanism 
that  Implies  that  the  reaction  proceeds  through  an  Inter¬ 
mediate  complex  that  Is  relatively  long-lived. 

♦  Work  supported  by  ONR 
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FIG.  2.  Comparison  of  S|x>ntancous 
and  laser  emission  spectra  from  XeCl. 
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Number  of  quasibound  levels  behind  a  rotational  barrier 
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The  number  of  quasibound  or  rotationally  predissociating  levels  of  diatomic  molecules  is  investi¬ 
gated  using  a  long-range  approximation  of  the  difference  between  the  Jeffreys-Wentzel-Kramers- 
Brillouin  estimates  of  the  vibrational  quantum  numbers  at  the  maximum  of  the  rotational  barrier 
and  at  the  dissociation  limit.  It  is  found  that  for  a  fixed  J  the  expected  number  of  quasibound  levels 
is  approximately  7/20,  essentially  independent  of  the  magnitude  and  form  of  the  potential  and  in¬ 
dependent  of  the  mass  of  the  nuclei. 


The  effective  interatomic  potential  for  a  rotating  dia¬ 
tomic  molecule,  Uj(R),  is  typically  represented  at  large R 
by  adding  the  rotational  kinetic  energy  to  a  single-term 
multipole  expansion,  giving  an  expression  of  the  form 

Uj(R)^1^JU-iri)/2iiR'^-C„/R''  .  (1) 

Such  a  potential  has  a  maximum  value  of 

at 

1/(b-2) 

,  (3) 

derived  from  the  condition  U}lRi,)=0.  We  are  interest¬ 
ed  in  estimating  the  number  of  quasibound  vibrational 
levels  that  can  exist  behind  this  barrier.  To  do  this  we 
follow  Stogryn  and  Hirschfelder'  and  Dickinson  and 
Bernstein^  and  construct  the  Jeffreys-Wentzel-Kramers- 
Brillouin  (JWKB)  estimate  for  the  vibrational  quantum 
number  at  the  barrier  maximum  E=E^,  and  compare 
it  with  the  vibrational  quantum  number  Uq  Ibe  dissoci¬ 
ation  limit  £=0.  Suppose  that  R^iEf,)  and  Ri(Ef,)—Rf, 
are  the  left-  and  right-hand  turning  points  for  £  =£(,  as 
shown  in  Fig.  1(a);  then, 

,  (4) 

where  a=V (2n) /(irU).  Similarly,  if  £i(0)  and  ^2^0) 
~£o  tbe  left-  and  right-hand  turning  points  for 
£  =s0,  we  have 


n/iC, 


Thus 

-11,(0) 

u, -Uo =a )]'^  d£ 

-il,(0) 

-[-UjlR)y^]dR 

+“/jt2<0)  [Ei,-Uj{R)y'^dR  .  (6) 

Following  the  approach  of  LeRoy  and  Bernstein,^  we 
plan  to  ignore  the  first  term  [supposing  that 
£|(£^)— J{|(0)  is  small  and  that  the  repulsive  wall  is 
steep]  and  to  replace  C/y(/{)  by  its  asymptotic  expansion. 
This  gives 

•'  0 

+  [-Uj{R)/E^y^]-'dy 

+pf^'’^'^'’[l-Uj{R)/Et]''^dy  .  11) 

where  we  have  substituted  and  y  =R/Rq. 

See  Fig.  Kb)  for  a  comparison  of  the  exact  [Eq.  (6)]  and 
long-range  approximation  [Eq.  (7)]  to  the  JWKB 
difference  integrand.  Deriving 


£o  — 


1/(11 -2) 


(8) 


from  the  condition  Uj(Rq)=0,  we  are  ready  to  perform  a 
number  of  back  substitutions. 

The  remarkable  result  of  these  substitutions  is  that  the 
molecular  parameters  n  and  C„  disappear  entirely,  and 
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FIG.  1.  (a)  Upper  curve:  Effective  potential  UjiR)  based  on 
a  Leonard-Jones  (12,6)  form,  (b)  Lower  curves:  JWKB 
difference  integrands  for  the  exact  potential  and  for  the  long- 
range  approximation. 


only  the  dependencies  on  n  and  J  remain.  We  first  note 
that  then  we  evaluate 

P=[J(J  + 1  )/y„  Y^/ir,  where 


(n/2-U 


(9) 


and  further  observe  that  —Vj(R)/Ei,=y„(y~''—y~^). 
Thus  we  have 

+  '  UO) 

where  ' 


/„  = 


I 

m/Vn 


HY„(y~''-y~'^)V'^\~^dy 


+  f"''[i+rn(y-''-y'^n'^dy 


(ID 


We  have  just  shown  that  the  number  of  quasibound  levels 
does  not  depend  on  fi  and  C„ . 

An  even  more  surprising  result  is  obtained  when  we 
evaluate  the  integral  /„.  Integrating  numerically  we  find 
the  results  given  in  Table  I,  and  we  conclude  that  the 
number  of  quasibound  vibrational  levels  is  effectively  in¬ 
dependent  of  n  as  well.  Also  included  in  Table  I  are  the 
numbers  derived  by  Dickinson  and  Bernstein^  for  n  =4, 
6,  and  8.  They  considered  potentials  for  which  the 
JWKB  integrals  for  Vi,  and  Uq  could  be  evaluated  explicit¬ 
ly  in  terms  of  elliptic  integrals  and  extracted  the  leading 
term  in  /  of  the  difference  between  them.  They  noticed 
the  near  independence  on  n,  but  their  formulation  did  not 
make  clear  that  the  expected  number  of  quasibound  lev¬ 
els  for  low  J, 


TABLE  I.  Numerical  values  for  the  integral  /,. 


n 

Present 

Ref.  2 

3 

0.045 1 

4 

0.0498 

0.0498 

5 

0.0498 

6 

0.0482 

0.0482 

7 

0.0462 

8 

0.0441 

0.044 

9 

0.0420 

lO 

0.0401 

A^qb(-^)  =  i’6  -''o=0-048±0.02[y(y  + 1  )]'^^ 

=//20  ,  (12) 

is  entirely  independent  of  the  molecule  under  considera¬ 
tion. 

To  assess  the  applicability  of  the  expressions  derived 
above,  we  have  examined  the  published  compilations  of 
quasibound  levels  in  which  it  is  claimed  that  all  the  quasi¬ 
bound  levels  have  been  calculated.  The  formula  for 
Nq^{J)  above  indicates  the  average  number  of  quasi¬ 
bound  levels  expected.  For  any  specific  J,  we  expect  to 
find  at  least  J/20— I  and  at  most  7/20-I-1  levels.  To 
evaluate  whether  the  number  of  levels  found  matches  our 
expectation  we  count  all  the  levels  expected  for  angular 
momenta  up  to  the  specific  value  of  J, 

Tqb<^J)=  "2  Nq^(L)c^J(J+2)/40  .  (13) 

L=0 

Figure  2  shows  a  graphical  comparison  of  the  numbers 


FIG.  2.  Graphical  comparison  of  expected  number  of  quasi¬ 
bound  levels  with  the  numbers  found  in  exact  calculations,  x ’s, 
HeH+  (Ref.  4);  +  's,  CH+  (Ref.  5);  diamonds,  Hj,  HD,  Dj  (Ref. 
6)  and  HT,  DT,  T2  (Ref.  7);  squares,  HgH  (Refs.  8  and  9);  trian¬ 
gles,  H2'''  (Refs.  lO  and  ID;  filled  circles,  more  than  one  super¬ 
imposed  symbol. 
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TABLE  II.  Numbers  of  quasibound  levels  in  ranges  of  J. 

•^ou* 

Tqb 

obs 

Tqb 

calc 

0-9 

10-19 

JVqb  in  ranges  of  Ts,  Tqb(J) 
20-29  30-39 

-TqsU 

40-49 

-10) 

50-59 

60-69 

^He'H+  * 

26 

20 

19 

2 

8 

10 

^He'H^  * 

27 

21 

20 

2 

7 

14 

’He'H+* 

34 

32 

31 

1 

8 

14 

7 

CH+'* 

36 

37 

34 

2 

6 

15 

14 

36 

39 

34 

2 

8 

17 

12 

38 

47 

38 

2 

8 

19 

18 

HgH" 

39 

43 

40 

5 

9 

10 

10 

41 

46 

44 

2 

6 

12 

23 

2 

HD' 

44 

65 

51 

3 

8 

18 

19 

7 

HT' 

47 

73 

58 

2 

8 

16 

30 

17 

Dj' 

54 

96 

76 

2 

7 

15 

17 

36 

8 

Dr 

59 

113 

90 

1 

8 

15 

25 

36 

28 

T/ 

67 

144 

116 

1 

7 

16 

24 

34 

45 

17 

Eq.  (13) 

2.5 

7.5 

12.5 

17.5 

22.5 

27.5 

32.5 

*  Reference  4. 
'’Reference  5. 
'  Reference  6. 


References  8  and  9. 

‘  References  10  and  1 1 . 
'Reference  7. 


of  quasibound  levels  found  in  exact  calculations  from  the 
literature  with  the  above  simple  formula  (we  have  used  a 
square-root  scale  on  the  y  axis  to  make  the  low  J  values 
more  discernible).  For  low  J  (say  <30),  the  present 
theory  is  clearly  consistent  with  the  exact  results.  The 
only  exceptional  case  is  HgH,  for  which  the  potential  is 
known  to  have  an  abrupt  change  of  form  at  R  =:4  A  (Ref. 
12)  and  thus  cannot  be  represented  by  a  single  long-range 
multipole  term.  For  higher  values  of  J,  the  present 
theory  consistently  underestimates  the  number  of  quasi¬ 
bound  levels,  because  of  the  increasing  contribution  of 
the  left-hand  turning  points  to  the  JWKB  difference  in¬ 


tegral;  however,  the  exact  calculations  continue  to  exhibit 
a  common  dependence  on  J. 

Table  II  shows  numerical  summaries  of  these  same 
comparisons.  The  levels  have  been  grouped  in  ranges  of 
ten  values  of  J.  Again,  we  see  that  essentially  perfect 
agreement  for  low  /  and  rather  good  agreement  even  for 
the  total  number  of  quasibound  levels  up  to  the  highest  J 
for  which  quasibound  levels  are  supported  by  the  poten¬ 
tial. 

This  work  was  supported  by  the  U.S.  Office  of  Naval 
Research. 
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COLLISIONAL  PROCESSES  IN  XeF(X) 


G.  Black,  L.  E.  Jusinski,  and  D.  L.  Huestis 
Molecular  Physics  Laboratory 
SRI  International 
Menlo  Paik,  CA  94025 

ABSTRACT 

Collision-induced  dissociation  of  XeF(X)  has  been  studied  with  He,  Ne,  Ar,  Kr,  Xe,  N2, 
SFg,  and  XeF2  as  collision  partners,  giving  dissociation  rate  constants  of  0.58, 0.62, 0.76, 0.68, 
0.75, 1.13,  1.05,  and  7.3  (±10%)  x  10-12  cm^  molec  -1  s'l  respectively.  The  values  for  He  and 
Ne  are  in  reasonable  agreement  with  those  found  previously  [Fulghum  et  al.  Appl.  Phys.  Lett.  35, 
247  (1979)].  Except  for  XeF2,  there  is  only  slight  dependence  on  the  nature  of  the  collision 
partner.  By  following  the  approach  to  association/dissociation  equilibrium  in  Xe,  we  are  able  to 
determine  that  the  yields  of  bound  XeF(X)  from  photodissociation  of  XeF2  at  193  and  248  nm  are 
only  0.4%  and  0.08%  respectively,  with  fragmentation  to  give  three  atoms  the  dominant  pathway. 
Vibrational  relaxation  is  observed  and  studied  in  detail  for  collisions  with  XeF2,  with  which  it  is 
found  to  be  very  fast. 
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INTRODUCTION 


The  XeF(B-X)  laser  operates  in  the  near  ultraviolet  on  a  number  of  bound-bound 
transitions,  as  listed  in  Table  I  and  illustrated  in  Figure  1,  terminating  on  various  vibrational  levels 
of  the  ground  X22+  state.  The  ground  state  of  XeF  is  weakly  bound  (Dg  =  1 175  cm'l), 
supporting  a  total  of  15  ±  2  closely  spaced  vitffational  levels  (tOt  =  226  cm  cm-l).  [2,3]  In  the 
absence  of  collisional  processes  depleting  these  levels,  the  ground-state  population  would 
accumulate  and  eventually  exceed  that  in  the  upper  laser  level,  resulting  in  termination  of  laser 
action.  Vibrational  relaxation  in  the  lower  levels  is  insufficient  to  adequately  deplete  the 
population,  since  the  small  vibrational  spacing  leads  to  estimates  that  9.6%  and  4.4%  of  the  ground 
state  molecules  will  be  in  v"  =  2  and  v"  =  3,  respectively,  in  equilibrium  at  room  temperature  [4]. 
On  the  other  hand,  the  binding  energy  of  the  ground  state  of  XeF  is  sufficiently  weak  that  the 
lower  level  population  can  be  depleted  by  collision-induced  dissociation,  which  removes  the 
ground  state  entirely  (rather  than  simply  changing  its  vibrational  distribution). 

Determination  of  the  rates  of  collision-induced  dissociation  is  a  complicated  task  with  an 
extensive  literature  (e.g.,  see  Dove  et  al.  [5]).  Experimental  measurements  exist  on  a  number  of 
systems,  usually  from  shock  tube  studies.  These  have  been  modeled  by  master  equation 
treatments  using  semiempirical  microscopic  rate  coefficients.  In  a  number  of  cases,  ab  initio 
collision  dynamics  studies  helped  supply  these  rate  coefficients.  For  XeF(X),  the  previous  work 
consists  of  (1)  modeling  studies  [8,9,10],  (2)  inferences  from  laser  gain  and  fluorescence 
measurements  [4,6,7],  and  (3)  direct  measurements  of  vibrationally  resolved  collision-induced 
dissociation  by  He  and  Ne  [1 1,12,13].  In  addition  to  these  gases,  collision-induced  dissociation 
has  been  studied  with  At,  Kr,  Xe,  N2,  SFg,  and  XeF2  in  the  work  described  here.  We  also  report 
some  measurements  of  vibrational  relaxation  in  XeF(X)  by  XeF2. 
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EXPERIMENTAL 


The  source  of  the  XeF(X)  for  these  experiments  was  photodissociation  of  XeF2. 
Photodissociation  was  carried  out  at  both  193  nm  (ArF)  and  248  nm  (KrF)  using  an  excimer  laser 
(Lambda  Physik  EMG  102  in  the  unstable  resonator  configuration).  A  schematic  of  the  apparatus 
is  shown  in  Figure  2.  A  similar  approach  has  been  used  previously  in  our  laboratory  [14,15]  and 
elsewhere  [12,13].  At  193  nm,  the  absorption  cross  section  of  XeF2  is  6  x  10  cm^  and  at 
248  nm,  a  factor  of  3-4  smaller  [16,17].  The  apertured  excimer  beam  (6-mm  diameter )  irradiated 
the  sample  with  typically  10  mJ  at  KrF  and  2  mJ  at  ArF  in  a  10-ns  pulse,  dissociating  <1%  of  the 
XeF2.  The  path  length  through  the  cell  was  *30  cm  and  the  cell  volume  *300  cm^.  Gas  mixtures 
were  made  up  in  a  stainless  steel  mixing  tank  (*  4  liter  volume)  and  then  slowly  flowed  through 
the  cell  to  avoid  buildup  of  photolysis  produas. 

The  XeF(X)  ground  state  was  detected  by  LIF  on  the  B  <-  X  transition.  Observations 
were  made  on  v"  =  0-5,  in  all  cases  pumping  up  to  the  v*  =  0  level  and  monitoring  either  (0,2)  or 
(0,3)  emission.  The  pumping  source  was  a  Quanta-Ray  Nd:Yag-dye  laser  system  that,  after 
doubling  provided  *5  mJ  in  the  cell  in  the  340  -  360-nm  region  of  the  B-X  (0,0-5)  transitions  in  a 
pulse  of  *5-ns  duration  with  a  linewidth  of  *0.5  cm'L  For  these  measurements,  the  lasers  were 
operated  at  10  Hz  and  were  adjusted  to  overlap  spatially  along  the  length  of  the  cell.  The  time 
delay  between  the  photodissociation  and  probe  lasers  was  set  by  a  digital  delay  generator 
programmed  in  a  nonlinear  ramp  from  a  VAX  1 1/750  computer  in  another  room.  [DID  WE 
CALffiRATE  THE  TIME  BASE?] 

The  detection  system,  situated  perpendicular  to  the  laser  beams,  consisted  of  a  0.35-m 
monochromator  equipped  with  an  RCA  C31034A  photomultiplier.  For  most  measurements,  the 
monochromator  slits  were  set  at  200  pm  providing  a  spectral  bandwidth  of  0.4  nm.  The 

photomultiplier  output  passed  to  a  boxcar  averager  (Stanford  Research  Systems  Model  SR250)  and 


G-3 


then  to  the  computer  for  storage  and  subsequent  analysis.  [WHAT  WAS  THE  SETTING  ON 
THE  BOXCAR?  WINDOW  POSITION  AND  WIDTH?] 

Most  measurements  were  made  with  the  dye  laser  wavelengths  close  to  the  band  origins,  in 
which  case  a  number  of  rotational  lines  overlapped  the  dye  output  and  the  largest  signal  was 
obtained.  Because  this  approach  favors  low  values  of  J",  many  measurements  were  made  with  the 
dye  pump  wavelength  shifted  several  A  to  the  red  of  the  band  origin,  where  higher  J"  levels  would 
be  excited.  Since  no  significant  difference  could  be  found  in  the  results,  rotational  energy 
apparently  was  rapidly  equilibrated  within  each  vibrational  level  (i.e.  faster  than  vibrational 
relaxation  or  collision-induced  dissociation).  Because  the  dye  laser  power  is  sufficient  to  saturate 
the  B<-X  transitions  within  its  »0.5cm*l  bandwidth,  the  intensity  of  the  LIE  signal  had  only  very 
slight  dependence  on  dye  laser  energy.  Hence,  the  signal  intensity  is  relatively  independent  of  the 
absorption  cross  section  for  the  pumped  transition  and  effectively  proportional  to  the  ground-state 
population. 


G-4 


RESULTS  AND  DISCUSSION 


COLLISION-INDUCED  DISSOCIATION 

Figure  3  shows  the  temporal  behavior  of  the  XeF(v"=0)  level  when  XeF2  at  a  pressure  of 
0.5  torr  is  photodissociated  at  248  nm.  The  solid  line  is  a  two-exponential  least-squares  fit.  The 
monochromator  is  set  to  observe  the  (0,3)  B<— X  emission,  and  the  spike  at  the  origin  originates 
from  the  production  of  this  emission  feature  by  the  KrF  laser  alone.  This  interfering  emission  was 
observed  to  exhibit  a  squared  dependence  on  the  KrF  intensity,  as  expected,  since  direct 
production  of  XeF(B)  is  not  energetically  possible  at  248-nm.  In  previous  work  in  this  laboratory 
[15]  the  production  of  B— >X  emission  in  248  nm  photodissociation  of  XeF2  was  attributed  to  the 
sequence  of  reactions 


XeF2  +  hv  (248  nm)  XeF(X)  +  F 

(1) 

XeF(X)  +  hv  (248  nm)  ^  XeF(D)v'=7.8.9 

(2) 

XeF(D)  +  M  ^  XeF(B,C)  +  M 

(3) 

It  was  also  shown  previously  [14]  that  photodissociation  of  XeF2  by  the  KrF  laser  leads  to  a 
strongly  nonthermal  vibrational  distribution  in  XeF(X).  The  behavior  shown  in  Figure  3,  in  which 
little  v"  =  0  is  produced  directly  but  builds  up  rapidly  by  vibrational  cascade  from  higher  levels 
followed  by  a  slower  removal  by  collision-induced  dissociation,  supports  this  conclusion. 

Figure  4  shows  a  similar  experiment  using  193  nm  to  photodissociate  the  XeF2. 
Interference  from  the  ArF  laser  is  also  apparent.  In  this  case,  the  B->X  (0,3)  emission  produced 
by  the  193-nm  radiation  exhibited  a  first-power  dependence  on  excimer  energy,  since  direct 
production  of  XeF(B)  is  possible  at  this  wavelength.  Figure  4  shows  that,  in  this  case,  any  rise  of 
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XeF(v"=0)  due  to  vibrational  relaxation  was  small  compared  to  the  direct  production  by 
photodissociation  at  193  nm.  This  finding  conflicts  with  earlier  work  [12]  on  this  same  system. 

As  with  248  nm,  the  subsequent  slow  decay  reflects  collision-induced  dissociation. 

Confirmation  of  this  view  was  obtained  from  observations  of  the  other  v"  levels.  Figure  5 
shows  temporal  profiles  of  v"  =  0-3  using  KrF  photodissociation.  Up  to  v"  =  3,  effectively  the 
same  slow  final  decay  could  be  seen,  contirming  collision-induced  dissociation  from  a  relaxed  or 
steady-state  distribution  of  vibrational  levels.  For  v"  =  4  and  5,  only  the  much  faster  initial 
disappearance  could  be  detected.  Thus,  we  interpret  the  fast  initial  decays  of  v"  >1  and  build  up  of 
v"  =  0  in  terms  of  vibrational  relaxation,  (see  below).  The  steady-state  vibrational  distribution  is 
expected  to  be  Boltzmann-like  for  low  values  of  v,  with  the  higher  levels  (within  kT  of  the 
dissociation  limit)  suppressed  owing  to  more  rapid  collision-induced  dissociation  [18],  The 
experimentally  observed  rate  of  collision-induced  dissociation  is  then  an  average  of  the  state- 
specific  rates  weighted  by  the  steady-state  vibrational  populations. 

Figure  6  shows  the  rate  of  collision-induced  dissociation  as  a  function  of  XeF2  pressure  for 
experiments  with  both  248-  and  193-nm  photodissociation  using  observation  on  only  the  v"  =  0 
level.  With  248-nm  photodissociation,  where  both  the  rise  and  faU  of  the  XeF(X)v"=o  signal  could 
be  observed,  the  addition  of  other  gases  was  observed  to  increase  both  the  rate  of  v"  =  0 
production  by  vibrational  relaxation  and  its  rate  of  removal  by  collision-induced  dissociation.  The 
effect  of  He  on  the  rate  of  collision-induced  dissociation  (after  correcting  for  the  effect  of  XeF2)  is 
shown  in  Figure  7.  The  other  gases  studied  gave  similar  effects,  and  the  resulting  rate  coefficients 
are  shown  in  Table  II. 

As  can  be  seen  from  the  results  in  Table  11,  except  for  XeF2,  there  is  very  little  dependence 
of  the  rate  coefficients  on  the  nature  of  the  collision  partner.  This  finding  agrees  with  theoretical 
expectations  [4].  Qearly,  specific  chemical  effects  are  associated  with  the  much  higher  rates  for 
XeF2.  Within  the  respective  error  bars,  the  value  for  He  is  in  agreement  with  the  value  of  (0.44 
0. 10)  X  10-12  cm^  molec'l  si  determined  by  Fulghum  et  al.  [12],  but  it  is  a  factor  of  2  larger  than 
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their  subsequent  value  of  0.31  x  10‘12.  The  intercept  value  in  Figure  2  of  this  earlier  work  [12], 
which  should  be  due  to  collision-induced  dissociation  by  XeF2,  gives  a  very  large  rate  coefficient, 
»2  X  10'*l  cm3  molec  *  S'*  when  interpreted  in  this  way.  Perhaps  the  accumulation  of  photolysis 
products,  resulting  from  the  lack  of  a  flow  system  in  this  earlier  work  is  the  source  of  the 
disagreement  (despite  the  use  of  very  low  photolysis  intensities  to  minimize  the  extent  of 
photodecomposition).  Our  rate  coefficient  for  Ne  is  somewhat  larger  than  the  0.33  0.2  x  10' *2 
found  by  Fughum  et  al.  and  more  than  three  times  the  value  of  0.18  x  10*^2  inferred  independently 
by  Tang  et  al.  [8,10]  and  Rokni  et  al.  [9]  from  laser  gain  and  from  sidelight-fluorescence 
suppression  measurements,  respectively. 

XENON  EXPERIMENTS 

With  xenon  as  the  added  gas,  the  observations  were  somewhat  different,  as  Figure  8 
shows.  A  rise  is  seen  to  a  steady-state  intensity  lasting  several  milliscads.  The  steady-state  signal 
arises  from  the  steady-state  concentration  of  XeF(X)v"=o  established  in  the  equilibrium 

Xe  +  F+  Xe< — » XeF+  Xe  (4) 

and  the  rate  of  approach  to  this  steady  state  will  reflect  the  rate  of  collision-induced  dissociation  (if 
vibrational  relaxation  is  the  faster  process,  as  found  in  the  other  systems).  On  this  basis,  the  value 
for  Xe  shown  in  Table  II  was  obtained.  The  Xe  dependence  of  the  rise  frequency  is  shown  in 
Figure  9.  The  steady-state  signal,  whose  magnitude  is  roughly  proportional  to  the  Xe  density  as 
expected  [at  high  density,  Xe  quenches  XeF(B)],  eventually  decays  as  F  atoms  and  XeF  molecules 
diffuse  out  from  the  common  volume  illuminated  by  the  two  laser  beams. 

It  is  interesting  to  compare  our  results  to  those  of  Appelman  and  Clyne  [19],  who  studied 
the  disappearance  of  F  atoms  in  a  mixture  of  Xe  and  Ar  buffer  gases.  They  attributed  this  loss  to 
three-body  formation  of  XeF  in  the  reaction 

Xe  +  F  +  Ar-^XeF  +  Ar  (5) 
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with  subsequent  loss  of  XeF  resulting  in  formation  of  XeF2.  They  performed  these  measurements 
under  a  single  set  of  conditions,  [Xe]  =  3  x  lO^^  cm-3  and  [Ar]  =  5  x  IQl^  cm'3,  and  observed  a 
decay  frequency  of  5.6*s*l.  Our  results  suggest  that  under  these  conditions  the  three-body 
formation  of  XeF  will  reach  equilibrium  with  collision-induced  dissociation  on  a  much  faster  time 
scale  (6  x  l(yl-s*l).  Thus,  the  net  loss  of  F  atoms  Appleman  and  Clyne  observed  must  be 
controlled  by  other  reactions.  Possible  candidates  are 

XeF  +  F2  — ^  XeF2  +  F  (6) 

and 

XeF  +  F  — >  F2  +  Xe  (7) 

Our  Xe  studies  also  allow  us  to  estimate  the  branching  ratio  fw  the  process 

XeF2  +  hv  -4  {/xcF  +  +  F  +  F/  (8) 

We  can  use  steady-state  concentration  of  XeF  at  long  time  to  calibrate  the  XeF  produced  directly  by 
the  photolysis  laser.  Figure  10  shows  temporal  profiles  of  XeF(X)v"=o  resulting  from  193-nm 
photolysis  at  two  relatively  low  Xe  pressures  (25  and  50  torr).  In  both  cases,  the  amount  of 
XeF(X)  produced  by  the  laser  is  comparable  to  that  in  the  eventual  steady  state.  For  pressures  >50 
torr,  the  signal  rises  to  the  steady-state  value,  where  as  for  [Xe]  <25  torr,  the  signal  falls.  A 
similar  relationship  holds  for  the  higher  vibrational  levels,  v"  =  1-4,  confirming  that  vibrational 
relaxation  is  much  faster  than  collision-induced  dissociation.  [COULD  THIS  JUST  MEAN  THAT 
THE  INITIAL  DISTRIBUTION  IS  SIMILAR  TO  THE  STEAD-STATE  DISTRIBUTION?] 

Thus,  we  conclude  that  for  193-nm  photolysis  at  Xe  pressure  of  about  35  torr. 


[XeF(X)]initial  “  [XeF(X)]eqQiiibriuin 


(9) 


when  summed  over  the  first  five  vibrational  levels.  Our  previous  calculation  of  the  equibrium 
constant  [4]  allows  us  to  infer  that  at  35  torr 

_  Q  ^®^^^^eqiiilibriuin  *  0.002  f^equilibriiun  (10) 

Hence,  we  can  conclude  that  the  production  of  XeF(X)  is  a  minor  product  of  XeF2 
photolysis  at  193  nm,  with  a  quantum  yield  of  0.4  0.1%  (note  that  two  F  atoms  are  produced  for 
each  dissociation).  Furthermore,  the  direct  yield  of  XeF  (B)  at  193  nm  must  be  very  small 
(<0.01  %). 

A  similar  investigation  with  248-nm  photodissociation  gave  an  even  smaller  yield  of  the 
v"  =  0  4  bound  levels  (quantum  yield  ==0.08%).  Within  the  experimental  errors,  the  bound  levels 
above  v"  =  can  be  neglected,  and  the  above  numbers  represent  the  total  yield  of  bound  levels. 

VIBRATIONAL  RELAXATION 

As  mentioned  above  and  illustrated  in  Figure  5,  with  XeF2  alone  and  with  the  KrF  laser, 
the  v"  =  0  rise  time  and  the  initial  decays  of  v"  =  2-5  were  sufficiently  slow  that  they  could  be 
measured  despite  interference  from  excimer-produced  fluorescence.  Figure  1 1  shows  the  results 
as  functions  of  XeF2  density,  and  Table  III  lists  the  slopes.  The  rate  of  formation  of  v"  =  0  should 
correspond  to  an  "effective"  rate  coefficient  for  vibrational  relaxation  from  all  the  higher  levels. 
Because  of  the  suggestive  pattern  exhibited  by  the  effective  rate  coefficients  listed  in  Table  HI,  it  is 
tempting  to  infer  that  vibrational  relaxation  occurs  in  Av"  =  1  steps,  for  which  the  numbers  listed  in 
Table  HI  would  then  be  the  state-to-state  rate  coefficients  (the  rise  of  v"  =  0  being  the  relaxation  of 
v"=l).  This  inference  poses  two  problems  .  First,  the  rates  observed  are  very  high,  and  could 
result  from  a  "chemical"  rather  than  "physical"  relaxation  mechanism.  A  possibility  is  F  atom 
exchange. 
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XeF  +  XeF2  — >  XeF2  +  XeF 


(11) 


which  might  be  expected  to  accomplish  energy  scrambling.  Second,  in  each  case,  the  separation 
between  the  neighboring  vibrational  levels  is  less  than  kT,  so  that  in  addition  to  cascade  filling 
from  above,  any  given  vibrational  level  is  also  significantly  repopulated  by  vibrations',  excitation 
from  the  next  lower  level.  The  net  result  is  that  multiple-exponential  decays  should  always  be 
observed  and  the  rate  coefficients  derived  from  fitting  the  data  should  depend  on  the  initial 
conditions,  on  the  rate  of  decay  of  the  resulting  equilibrium  distribution,  and  on  the  signal-to-noise 
ratio. 

In  an  attempt  to  quantify  these  arguments,  we  tried  to  simulate  the  data  shown  in  Figure  5 
using  a  model  similar  to  that  developed  previously  [4,7]  (The  vibration  and  dissociation 
parameters,  Q  and  C,  respectively,  were  increased  by  a  factor  of  17,  and  a  constant,  4.4  x  10'12, 
was  added  to  each  state-specific  dissociation  rate  coefficient).  The  result,  assuming  the  population 
is  initially  in  v"  =  3  is  shown  in  Figure  12.  Table  ID  lists  the  model  rate  coefficients  for  single- 
step  vibrational  relaxation  and  excitation,  along  with  the  total-removal  rate  coefficient,  for  each 
level.  We  are  encouraged  in  our  interpretation  of  the  experiments  by  the  agreement  between  the 
observed  effective  rate  coefficients  and  the  Av"  =  1  values  from  the  model,  despite  the  fact  that  the 
riKxlel  predicts  that  the  total  rates  of  removal  from  each  level  are  a  factor  of  2  to  3  larger. 
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Table  1 

XeF  LASER  TRANSITIONS  [1] 


348.8  nm 

v’ 

=  2  -4  v"  =  5 

351.1  nm 

v’ 

=  1  v"  =  4 

-351.3  nm 

v* 

=  0  v"  =  2 

353.2  nm 

v’ 

=  0  v"  =  3 

353.2  nm 

V* 

=  1  -4  v”  =  6 

Gas 

Table  2 

Rate  Coefficients  for  Collision-Induced 

Dissociation  of  XeF(X)  at  300  K 

Rate  Coefficients  (10*'' 2  cm3  ■nolec*’!  s*"*) 

Present  Work 

Previous  Work 

He 

.58  ±  0.06 

0.44±0.ia  0.3lb 

Ne 

.62  ±  0.06 

0.33  ±0.23  0.1 8C 

Ar 

.76  ±  0.08 

Kr 

.68  ±  0.07 

Xe 

.75  ±  0.08 

N2 

1.13  ±0.11 

SF6 

1.0510.15 

XeF2 

7.3  ±  0.8 

CM 

t 

®  Reference  [12] 

^  Reference  [13] 

®  References  [8,9,10] 

^  Estimated  from  intercept  in  Figure  2  from  Reference  [12] 
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Table  3 


Effective  Rate  Coefficients  for  Vibrational  Relaxation 
of  XeF(X)  by  XeF2  at  300  K 


Rate  Coefficients  from  Model^ 


Process 

Observed 

Rate^ 

v" 

> 

II 

+ 

II 

> 

< 

Total 

0 

0.4 

0.5 

Rise  of  v"=0 

0.78 

1 

1.1 

0.6 

1.8 

Decay  of  v"=2 

1.2 

2 

1.4 

0.9 

2.5 

Decay  of  v"=3 

1.7 

3 

1.9 

1.4 

3.7 

Decay  of  v"=4 

2.8 

4 

2.6 

2.1 

5.8 

Decay  of  v"=5 

3.4 

5 

3.6 

3.3 

10.2 

3  In  units  of  10*"*®  cm^  molec'*  s'"* 
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FIGURE  CAPTIONS 


1 .  Schematic  XeF  energy  level  diagram. 

2 .  Schematic  of  the  experimental  apparatus  for  XeF  production  by  photodissociation  of  XeF2 
and  for  monitoring  its  decay  by  LIF. 

3 .  XeF(X)v"  =  0  production  by  248-nm  photodissociation  of  XeF2  at  0.5 1  torr  and  its 
subsequent  disappearance  by  collision-induced  dissociation. 

4.  XeF(X)v"  =  0  production  by  193-nm  photodissociation  of  XeF2  at  0.42  torr  and  its 
subsequent  disappearance  by  collision-induced  dissociation. 

5 .  Comparison  of  time  histories  of  XeF(X)v"=o,i.2,3  populations  following  KrF  laser  photolysis 
of  0.66  torr  of  XeF2. 

6.  XeF(X)v"=o  decay  rate  by  collision-induced  dissociation  versus  concentration  of  XeF2. 

7 .  XeF(X)v"=o  decay  rate  by  collision-induced  dissociation  versus  concentration  of  He. 

8 .  Temporal  profile  of  XeF(X)v"=o  in  the  248-nm  photodissociation  of  XeF2  at  1.0  torr  in  the 
presence  of  50  torr  of  Xe. 

9 .  Rate  of  collision-induced  dissociation  of  XeF(X)  by  Xe,  derived  from  the  rate  of  rise  of 
XeF(X)v"=o- 

1 0.  Temporal  profile  of  XeF(X)v"=o  in  the  193-nm  photodissociation  of  XeF2  at  0.33  torr  in  the 
presence  of  50  torr  of  Xe  (upper  curve)  and  at  0.50  torr  in  the  presence  of  25  torr  of  Xe 
(lower  curve). 


G-15 


1 1 .  Fast  decay  rates  of  XeF(X)v"=2-5  and  vibrational  relaxation  rate  into  XeF(X)v"=o  versus 
concentration  of  XeF2. 

1 2.  Comparison  of  model  calculations  with  observed  XeF(X)v"=o,l^,3  populations  following 
KrF  laser  photolysis  of  0.66  torr  of  XeF2. 
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Figure  1 .  Schematic  XeF  energy  level  diagram. 
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Figure  2.  Schematic  of  the  experimental  apparatus  for  XeF  production  by 
photodissociation  of  XeFa  and  for  monitoring  its  decay  by  LIF. 
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Figure  3.  XeF(X)v*.o  production  by  248'nm  photodissociation  of  XeF2  at 
0.51  torr  and  its  subsequent  disappearance  by  collision-induced 
dissociation. 
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Figure  4.  XeF(X)v-o  production  by  193-nm  photodissociation  of  XeF2  at 
0.42  torr  and  its  subsequent  disappearance  by  collision-induced 
dissociation. 
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Figure  5.  Comparison  of  time  histories  of  XeF(X)y.^  ^  2, 3  populations  following 
KrF  laser  photolysis  of  0.66  torr  of  XeF2. 
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Figure  6.  XeF(X)v'.o  decay  rate  by  collision-induced  dissociation  versus 
concentration  of  XeF2. 
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Figure  7.  XeF(X)v*.o  decay  rate  by  collision-induced  dissociation  versus 
concentration  of  He. 
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Figure  8.  Temporal  profile  of  XeF(X)v*.o  *n  the  248-nm  photodissociation  of 
XeF2  at  1 .0  torr  in  the  presence  of  50  torr  of  Xe. 


G-24 


Figure  9.  Rate  of  collision-induced  dissociation  of  XeF(X)  by  Xe, 
derived  from  the  rate  of  rise  of  XeF(X)v’^). 
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Figure  10.  Temporal  profile  of  XeF  {X)y-  ^  in  *he  193-nm  photodissociation  of 
XeF2  at  0.33  torr  in  the  presence  of  50  torr  of  Xe  (upper  curve)  and  at 
0.50  torr  in  the  presence  of  25  torr  of  Xe  (lower  curve). 
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Figure  1 1 .  Fast  decay  rates  of  XeF(X)y-„  2-5  ^^d  vibrational  relaxation  rate  into 
XeF(X)y-.o  versus  concentration  of  XeF2. 
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Figure  12.  Comparison  of  model  calculations  with  observed  XeF(X)  y- ^  o  i  2  3 

populations  following  KrF  laser  photolysis  of  0.66  torr  of  XeF2. 
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ABSTRACT 

We  have  measured  the  angular  and  density  dependence  of  the  Dicke 
narrowed  vibrational  line  shapes  for  stimulated  Raman  scattering  in 
H2.  Measurements  were  made  on  the  Q(l)  transition  for  angles  of 
0-165  deg  and  for  densities  of  1-25  amagat.  We  have  found  an  empiri¬ 
cal  formula  that  agrees  with  the  results  of  these  linewidth  measure¬ 
ments  to  within  5%.  This  empirical  formula  will  be  useful  for 
modeling  the  Reiman  gain  coefficient  as  a  function  of  angle  for  Raman 
beam  clean-up  and  aperture  combining. 

INTRODUCTION 

Raman  beam  clean-up  and  aperture  combining  is  currently  a  topic 
of  intense  research  within  the  laser  community.  This  technique 
improves  the  spatial  mode  of  coherent  light  sources  by  using  stimu¬ 
lated  Raman  scattering  in  H2  to  amplify  a  high  quality  stokes  beam. 
Energy  is  simultaneously  transferred  from  several  high  intensity, 
poor  quality  pump  beams  to  the  high  quality  Stokes  seed  beam  by 
crossing  the  pump  beams  at  small  angles  to  the  stokes  beam.  Thus  the 
angular  dependence  of  the  Raman  gain  needs  to  be  well  characterized 
before  this  type  of  beam  clean-up  process  can  be  accurately 
modeled.  For  pump  and  stokes  beams  that  are  linearly  polarized  in 
the  same  directions,  the  angular  dependence  of  the  gain  comes  only 
from  the  angular  dependence  of  the  linewidth  (the  peak  value  of  the 
steady-state  gain  coefficient  is  inversely  proportional  to  line- 
width)  .  Previous  experimental  work  reports  linewidth  measurements 
for  only  the  forward  and  backward  directions.^  This  work  presents 
experimental  results  for  the  angular  dependence  of  the  vibrational 
Raman  linewidth  in  gas  phase  H2. 

EXPERIMENTAL  DETAILS 

Using  a  standard  quasl-cw  stimulated  Raman  spectrometer,  we 
have  measured  the  linewidth  of  the  Q(l)  line  for  angles  of  0-165 
Jegrees  and  densities  of  1-25  amagat.  Both  the  cw  probe  (683  nm)  and 
the  pulsed  pump  (532  nm)  are  single  mode  lasers,  with  the  10  ns 
pulsed  laser  limiting  the  spectral  resolution  to  about  100  MHz.  This 
laser  linewidth  was  determined  from  the  excess  linewidth  not 
accounted  for  by  pressure  or  Doppler  broadening  in  a  copropagating 
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geometry,  where  the  lineshape  is  well  known,  and  from  assuming  that 
the  laser  line  shape  is  Gaussian.  For  small  angles  (less  than 
30  deg)  and  small  densities  (less  than  10  amagat) ,  this  laser  line- 
width  was  accounted  for  by  fitting  the  data  to  Voigt  line  shapes. 

For  large  angles  or  large  densities  the  100  MHz  width  of  the  laser  is 
negligible.  These  data  were  fit  to  pure  Lorentzlan  line  shapes. 

EMPIRICAL  MODEL 

The  density  dependence  of  these  measured  linewidths  was  compared 
to  an  empirical  model  given  by 


Ai/ 


+  Bp 


(1) 


3 

and  to  the  diffusion  model  given  by 

Ai/  -  ^  +  Bp  (2) 

where  ii/p  is  the  Doppler  broadened  linewidth,  B  is  the  density 
broadening  coefficient,  p  is  the  density.  The  quantity  A  is  given  by 

where  0  is  the  crossing  angle,  D  is  the  self-diffusion  coefficient, 
and  kp,  k^,  and  kj^  are  the  pump,  stokes  (probe),  and  Raman  shift  wave 
numbers  respectively. 

We  find  that  Eq.  1,  with  n  -  1.5,  gives  the  best  fit  to  our  data 
and  predicts  the  linewidths  to  within  5%  of  our  measurements  for  all 
angles  and  densities.  At  densities  below  1  amagat,  the  differences 
between  the  data  and  the  model  become  larger.  Because  Eq.  1  is 
empirical,  no  physical  significance  should  be  attached  to  the  value 
of  1.5  for  n.  The  experimental  data  and  the  empirical  fit  of  Eq.  1 
for  four  different  crossing  angles  is  illustrated  in  Fig.  1. 

Although  useful  for  modeling  purposes,  the  empirical  Eq.  1  does 
little  to  help  the  understanding  of  the  physics  of  the  collisional 
contribution  to  the  line  shape.  However,  line  shape  measurements 
with  uncertainties  of  less  than  1%  are  necessary  to  compare  with  more 
physically  meaningful  models  such  as  the  Galatry  and  related  line 
shapes.^  We  plan  to  make  improved  measurements  (modifying  our  pulsed 
pump  laser  with  a  single  mode  seed  laser)  to  make  this  comparison. 
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Figure  1.  Comparison  of  the  experimental  data  (symbols)  and  the 

empirical  model  (solid  and  dashed  lines)  of  Eq.  1,  with 
n  -  1.5. 
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VIBRATIONAL  RAMAN  CROSS  SECTIONS  OF  THE  FUNDAMENTAL  AND 
OVERTONE  TRANSITIONS  IN  H2,  D2,  AND  N2 


VIBRATIONAL  RAMAN  CROSS  SECTIONS 
OF  THE  FUNDAMENTAL  AND  OVERTONE 
TRANSITIONS  IN  H2,  D2,  and  N2 

G.  C.  Herring, t  and  William  K.  Bischeltt 
Molecular  Physics  Labcaatory 
SRI  International,  Menlo  Park,  CA  94025 

and 

Winifred  M.  Huo 
NASA  Ames  Research  Center 
Moffett  Field,  CA  94035 

ABSTRACT 

Theoretical  and  experimental  results  are  presented  for  Raman  scattering  cross 
sections  of  the  vibrational  fundamental  overtone  transitions  in  gas  phase  H2,  D2,  and  N2. 
Cross  sections  for  the  first  overtones  of  H2,  D2,  and  N2  have  been  measured  at  488  nm, 
while  cross  sections  for  the  fundamental  and  first  two  overtones  of  H2  and  D2  have  been 
calculated  as  a  fiinction  of  wavelength  from  ab  initio  theory.  Fw  H2,  the  first  overtone 
experimental  result  of  5.1  ±  1.5  x  10‘33  cm2  sr*,  is  in  reasonable  agreement  with  the 
theoretical  value  of  4.0  x  10*33.  xhe  measured  value  for  N2  of  1.3  ±  0.5  x  10*33  is  to  be 
compared  with  the  theoretical  value  of  0.64  x  10*33  of  Langhoff  ct  al.  The  role  of 
nonadiabatic  effects,  arising  from  the  breakdown  of  the  Bom-Oppenheimer  approximation, 
is  discussed  for  the  overtone  transitions. 


t  Current  Address:  NASA  Langley  Research  Center,  Hampton,  VA  23665-S22S 
tt  Current  Address:  Coherent,  Inc.  3210  Porter  Drive,  Palo  Alto,  CA  94304 
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INTRODUCTION 


In  the  vibrational  Raman  spectroscopy  of  diatomic  gases,  rigorous  selection  rules 
governing  the  change  of  vibrational  quantum  numbers  do  not  exist.  However,  most  of  the 
past  studies  have  been  restricted  to  the  O.  Q,  and  S  branches  of  the  fundamental 
transition.  This  limitation  is  mainly  due  to  the  fact  that  the  overtone  transitions  are 
orders  of  magnitude  weaker,  making  quantitative  measurements  difficult  Although  the 
Raman  overtone  frequencies  do  not  provide  any  new  information  about  the  vibrational 
constants  of  diatomic  molecules,  they  can  furnish  an  independent  check  of  the  constants 
derived  from  previous  quadrupole  experiments.^  More  importantly,  experimental  overtone 
cross  sections  can  be  used  to  verify  or  improve  the  approximations  employed  in  the 
calculation  of  Raman  transition  probabilities.  Since  overtone  transitions  are  related  to  the 
second  and  higher  derivatives  of  the  dynamic  polarizability,  they  generally  provide  a  more 
stringent  test  of  the  quality  of  the  calculation  than  the  fundamental  transition.  Overtone 
cross  sections  will  also  be  useful  for  experiments^  that  require  populating  the  higher  lying 
vibrational  levels  of  molecules. 

The  first  overtone  (Av  =  2)  has  been  measured^  in  solid  H2,  and  theoretical  woiit 
for  the  first  overtone  of  gas  phase  hJ  has  also  been  published,^  but  very  little  work  has 
been  reported  on  overtones  of  neutral  diatomic  molecules  in  the  gas  phase.  In  this  paper, 
measurements  of  Raman  scattering  cross  sections  for  the  first  vibrational  overtone  of  H2, 
D2.  and  N2  gases  are  reported.  Theoretical  calculations  of  the  transition  polarizabilities  of 
H2  and  D2  based  on  a  sum-over-states  method,  are  also  presented  for  transitions  from  the 
ground  state  (v  =  0-^  v'  =  0,1,2, 3)  and  the  first  vibration  state  (v  =  1  ->  v’  =  1,2,3, 4). 
These  results  are  compared  with  the  experiment  for  Av  =  1,2.  For  N2,  an  overtone  cross 
section  is  determined  from  the  dependence  of  thetnetical  static  polarizabilities^  on  the 
intemuclear  separation.  The  present  experimental  results  are  compared  with  Knippers  et 
al.9  and  Veirs,*®  who  have  also  observed  Raman  overtones  in  diatomic  gases.  No 
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previous  calculations  of  Raman  overtones  of  H2  and  N2  have  been  reported,  even  though 
accurate  calculations  of  the  dynamic  polarizability*!'^^  and  Raman  fundamental  cross 
sections*^**5  have  been  reported  for  H2.  A  number  of  theoretical  calculations  on  N2  have 
also  appeared  in  the  literature.  Good  agreements  with  experimental  data  are  found  in  the 
recent  calculation  of  Langhoff  et  al.^  on  the  static  polarizability  and  polarizability 


derivatives  of  N2. 


EXPERIMENTAL  DETAILS 


The  experimental  apparatus  used  for  the  overtone  measurements  is  similar  to  those 
used  routinely  in  spontaneous  Raman  spectroscopy.^  An  Ar  ion  laser  provided  from  one  to 
five  watts  of  linearly  polarized  light  at  wavelengths  of  458  nm,  488  nm,  or  515  nm.  This 
beam  was  focussed  into  a  cell  containing  H2  at  295  K  and  760  torr  or  into  the  room  air  for 
the  N2  experiments.  In  both  cases,  a  25-cm  focal  length  lens  was  used  to  produce  a  beam 
diameter  of  50  pm  at  the  focal  plane.  The  Raman-scattered  light  was  collected  over  a  solid 
angle  of  1  steradian  at  right  angles  to  both  the  propagation  direction  and  the  electric  field 
direction  of  the  pump  beam.  This  light  was  then  directed  through  a  polarization  scrambler 
and  into  a  0.7 5-m  Spex  1402  double  monochromator  (f/7).  Slit  widths  were  set  at  800 
pm,  giving  the  monochromator  a  resolution  of  5  A  at  visible  wavelengths.  The  Raman 
signal  was  detected  with  a  RCA  31034A  photomultiplier  tube  which  was  cooled  to  -30“  C 
and  biased  at  1300  volts.  Output  pulses  from  the  photomultiplier  tube  were  amplified  by  a 
factor  of  250  and  a  pulse  height  discriminator  was  set  to  reject  pulses  smaller  than  200  mv. 
Under  these  conditions,  the  dark  count  was  typically  75  counts/s,  while  the  Raman 
overtone  signals  were  about  5  and  20  counts/s  for  N2  and  H2  respectively.  Time  constants 
were  2  s  and  40  s  for  the  fundamental  and  overtone  scans,  respectively.  The  combined 
relative  efficiency  of  the  monochromator  and  photomultiplier  tube  was  calibrated  with  an 
Optronics  220C  spectral  standard  which  is  traceable  to  the  N6S. 

Overtone  cross  sections  were  experimentally  determined  by  measuring  the  ratio  of 
the  overtone  signal  to  the  fundamental  signal  and  using  this  ratio  with  the  previously 
measured  fundamental  cross  sections.^  The  Raman  signal  ratio  of  the  overtone  to  the 
fundamental  was  measured  by  sequentially  tuning  the  monochromator  over  these  two 
transitions  and  comparing  the  integrated  areas  of  the  two  transitions.  Corrections  in  all  of 
these  ratios  were  made  for  changes  in  laser  power,  monochromator  transmission  and 
dispersion,  and  photomultiplier  tube  sensitivity.  In  N2,  all  of  the  individual  rotational 
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transitions  of  the  Q  branch  are  so  closely  spaced  that  only  a  single  spectral  line  was 
observed  under  the  experimental  conditions  described  above.  Thus  the  nitrogen  cross 
section  measured  here  is  the  integrated  cross  section  over  the  entire  Q  branch.  In  H2,  the 
rotational  transitions  were  resolved  in  the  overtone,  but  were  not  resolved  in  the 
fundamental  transition.  Because  of  the  long  time  constants  used  in  the  overtone  scans, 
only  the  Q(0)  and  Q(l)  transitions  were  typically  recorded,  while  the  entire  integrated  Q 
branch  was  recorded  for  the  fundamental  scans.  Thus  an  additional  correction  (ratio  of 
Q(0)  plus  Q(l)  to  the  integated  Q  branch)  was  made  for  the  measured  ratio  in  H2  to  obtain 
the  total  integrated  Q  branch  ratio. 

Results  of  the  present  measurements  of  the  Raman  overtone  cross  sections  in  H2 
and  N2  are  shown  in  Table  I,  where  the  reported  values  are  the  average  of  five  separate 
measurements.  Uncertainties  for  each  individual  measurement  were  estimated  to  be  ±  20  % 
and  40  %  for  H2  and  N2  respectively.  The  uncerminties  shown  in  Table  I  for  this  work  are 
one  standard  deviation  of  the  average  for  the  five  different  measurements.  Previous 
measurements  of  the  fundamental  cross  sections  and  present  theoretical  results  are  also 
shown  in  Table  I. 

Our  current  overtone  measurements  have  varying  degrees  of  agreement  with  two 
previous  overtone  measurements.  Our  H2  result  is  about  30%  larger  than  the  result  of 
Veirs.  However,  this  difference  is  well  within  the  error  bars  of  both  experiments.  In 
contrast,  our  N2  and  D2  results  are  a  factor  of  3  greater  than  the  results  of  Knippers  et  al.^ 
Since  we  have  carefully  calibrated  the  wavelength  response  of  our  Raman  spectrometer, 
and  our  results  are  consistent  with  those  of  Veirs,  we  believe  our  determination  of  the 
cross  section  is  the  most  accurate. 

An  attempt  to  observe  the  second  overtone  (Av=3)  in  both  H2  and  N2  was  made 
using  the  364  nm  line  of  the  Ar**  ion  laser.  No  signal  was  found  in  either  case.  This  was 


not  surprising  since  calculations  show  that  the  signal  exacted  for  a  2: 1  signal-to-noise  ratio 
was  40  times  smaller  than  the  detection  limit. 
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THEORETICAL  CALCULATIONS 


In  the  case  of  a  linearly  polarized  pump  beam,  with  the  scattering  angle  90°  to  both 
the  polarization  vector  and  the  direction  of  incidence  (without  polarization  analyzer),  the 
Raman  cross  section  for  a  Q(J)  branch  transition  of  a  diatomic  molecule  in  the  state  is 
given  by^ 


aovjv-j 

30 


=  (2tc)4(vp-VR)4(2J  +  l) 


X 


OvVvj 


(1) 


where  v  and  v'  are  the  vibrational  quantum  numbers  of  the  initial  and  final  states,  J  is  the 

rotational  quantum  number,  Vp  is  the  frequency  (in  cm*^)  of  the  incident  beam  and  vr  is  the 

Raman  transition  frequency.  Also,  avj.v'j  and  Yvjyj  are  the  isotropic  and  anisotropic 

components  of  the  transition  polarizability.  Fot  v=v',  they  reduce  to  the  standard  deHnition 

of  the  isotropic  and  anisotropic  components  of  the  polarizabiliQr  of  a  molecule  in  the  (v  J) 

state.  Since  this  paper  deals  with  Q  branch  transition  only,  the  simpler  notation  aw'  and 

Yw'  will  be  used  and  their  J  dependence  will  be  implicitly  assumed. 

The  quantities  Ow'  and  Yw'  are  related  to  the  parallel  and  perpendicular  components  of  the 

transition  polarizabilities,  a„  ,  and  a , 

II  ,vv  ± ,  vv’ 


a 


vv 


5  (“ll.w 


+  2a 


(2) 


Yyv' 


“ll  .  -  “  1 

II  .vv’  1 


vv 
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Using  the  Kramers-Heisenberg  formula,  a 


,  can  be  expressed  as 


i!  .w 


X  <  VoVJ  I  m  I  Vivjjj  X  Vivjjj  I  H||  I  VovJ 
i  Vi 


■  ^ovj 


E: 


iviJj 


^ovj  ''s 


(4) 


where  |i||  denotes  the  parallel  component  of  the  dipole  operate,  y  is  the  rovibnmic 
wavefunction  of  the  molecule,  o  denotes  the  initial  (=rmal)  electronic  state,  and  Vs  is  the 
frequency  of  the  scattered  beam.  The  index  i  sums  over  all  electronic  states,  vj  the 
vibradonal  levels  of  i,  and  Ji  its  rotational  quantum  number.  The  summations  over  i  and  vi 
include  integration  over  the  continuum  portion  of  the  electronic  and  vibrational  manifold. 
The  summation  over  Jj,  (Ji  =  J,  J  ±  1),  has  been  carried  out  implicitly  in  the  derivation  of 
the  3-j  symbols  in  Eq.  (1),  with  the  assumption  that  the  dipole  matrix  elements  of  the  three 
possible  Ji  are  the  same.  Notice  that  Eq.  (4)  is  applicable  to  a  noolecule  without  a 
permanent  dipole  moment. 

In  Eq.  (4,)  nuclear  motions  are  included  explicitly  in  the  summation.  On  the  other 
hand,  in  the  Bom-Oppenheimer  approximatitMi.  the  electronic  polarizability  is  first 
calculated  in  the  fixed-nuclei  approximation  and  then  integrated  over  the  initial  and  final 
rovibronic  states: 

a||  (Vp,R)  =  I  i  <  <I>o(r,R)  I  mi  I  Oj  (r.R)  >r  <  (Dj  (r,R)  I  U||  I  <l>o  (r,R)  >!• 

X  ( - ^ ^ (5) 

l,Ei  (R)  -  Eo  (R)  -  Vp  Ei  (R)  -  Eo  (R)  +  VsJ  ’ 
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and 


“II.  vv’  =  (R)  I  “11  (Vp  R)  I  Xv(R)>- 


(6) 


In  Eq.  (5),  r  and  R  denote  the  electronic  and  nuclear  coordinates,  and  the  subscript  r 
attached  to  the  bracket  denotes  integration  over  electronic  coordinates  alone.  Also,  is  the 
fixed-nuclei  electronic  wavefuncrion  and  x  is  the  vibrational  wavefuncdon  of  the  electronic 
state  o.  The  difference  between  Bom-OppenheinKr  and  nonadiabadc  treatments  have  been 
discussed  in  detail  by  Bishop  and  Cheung.  1  *  In  their  calculadon  on  H2,  they  recast  the  ro- 
vibronic  problem  into  a  secular  equadon  problem.  TTieir  results  are  considered  intermediate 
between  the  fixed-nuclei  and  nonadiabatic  treatments.  The  calculadons  of  Rychlewsld.l^ 
Huo  and  Jaffe,^^  Ford  and  Browne,!^  and  Langhoff  et  al,®  are  all  based  on  the  Bom- 
Oppenheimer  approximadon,  even  though  different  methods  are  used  in  the  treatment  of  the 
electronic  part  of  the  problem. 

In  the  present  calculadon  on  H2,  a  method  which  is  a  combinadon  of  Eq.  (4)  and 
(5)  was  employed. f 6  From  a  previous  sum-over-state  calculation,^^  it  was  found  that 
approximately  75%  of  ooo  and  90%  of  aoi  come  from  the  contribudons  of  B^lJ  and 
C^riu  states,  the  two  lowest  excited  singlet  states  of  H2.  Thus,  in  this  study,  the 
nonadiabadc  effects  from  these  two  states  were  accounted  for  by  carrying  out  explicit 
summadon  over  their  ro-vibrational  levels,  as  given  in  Eq.  (4).  The  contributions  from  the 
remainder  of  the  electronic  states  were  calculated  using  the  Bom-Oppenheimer 
approximation,  as  in  Eq.  (5).  In  this  manner,  a  major  portion  of  the  ntmadiabatic  effects 
has  been  included.  The  electronic  wavefunedons  used  are  configuration-interaction  (Cl) 
wavefunctions  constructed  to  satisfy  the  Dalgamo-Epstein  condition,^^  which  enables  us  to 
replace  the  infinite  sum  over  electronic  states  by  a  finite  sum.  Tlie  details  of  the  electronic 
calculadons  have  been  described  previously.!^  The  vibrational  wavefuncdmis  were  the 
numerical  solutions  of  the  one-dimensional  Schrbdinger  equation  of  vibrational  motion. 
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obtained  using  the  Numerov -Cooley  technique.  The  potential  functions  for  the  X,  B,  and 
C  state  vibrations  were  taken  from  the  work  of  Kolos  and  Wolniewciz.*®'20 

A  summary  of  the  H2  Rayleigh  and  Raman  cross  sections  from  the  present 
calculation  is  given  in  Tables  II  and  in.  Average  polarizabilities  and  polarizability 
anisotropies  are  also  given  for  each  transition.  Table  II  contains  results  for  transitions 
originiating  in  v  =  O,  while  Table  HI  contains  results  for  transitions  originating  in  v  =  1, 
The  numbers  given  in  Table  11  are  used  in  Rg.  1  to  show  the  relative  wavelength 
dependence  of  the  square  of  the  polarizabilities  (proportional  to  the  scattering  cross  section 
for  the  different  scattering  processes  in  H2.  In  addition  to  illustrating  the  difference  in 
the  wavelength  dependence  of  the  different  processes,  Rg.  1  also  shows  the  preresonance 
enhancement  for  the  fundamental  with  the  2pa  state  at  90,(XX)  cm*^  All  cross 
sections  given  in  Tables  II  and  m  are  averaged  over  the  thermal  distribution  of  initial 
rotational  levels.  One  of  the  interesting  results  of  these  calculations  is  that  the  polarizability 
anisotropy  for  the  first  overtone  crosses  zero  between  250  and  300  nm.  Thus  the  cross 
sections  for  the  S  branch  of  the  v  =  0->2  transition  will  vanish  for  pump  lasers  in  this 
wavelength  region. 

Theoretical  cross  sections  at  488  nm  are  compared  to  experimental  cross  sections  in 
Table  I.  For  H2,  the  theoretical  overtone  values  of  Table  I  are  taken  directly  finran  Table  n. 
The  corresponding  quantity  for  N2,  calculated  from  the  theoretical  static  polarizability  of 
Langhoff  et  al.  using  Eq.  (6),  is  also  presented.  An  RKR  potential  is  used  for  the  N2 
vibration.^^  For  comparison  purposes,  previous  calculatitms  of  the  fundamental  cross 
sections  are  also  listed.  Since  Rychlewski,  and  Ford  and  Browne,  did  not  include  488  nm 
in  their  calculation  of  H2  fundamental,  the  tabulated  cross  sectitxis  are  obtained  by  spline 
fitting  the  frequency  dependence  of  their  published  ooi  and  Yol- 
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DISCUSSION 


It  is  illustrative  to  compare  the  experimental  cross  section  of  the  fundamental 
transition  with  theory  Hrst.  Overall,  good  agreement  between  theory  and  experiment  is 
obtained  for  both  molecules.  In  particular,  the  three  calculations  on  H2  using  the  adiabatic 
approximation  are  in  excellent  agreement  with  each  other.  Among  the  three,  the 
wavefunction  employed  by  Rychlewski  is  of  the  highest  quality.  The  nonadiabatic  result, 
from  an  unpublished  calculation  of  Huo,^^  is  »  10%  smaller,  and  is  in  somewhat  closer 
agreement  with  experiment.  The  electronic  wavefuncdons  used  in  the  ntxiadiabatic 
calculation  are  determined  in  a  manner  similar  to  those  used  in  Ref.  14.  The  difference 
between  the  results  in  Ref.  14  and  16  is  partially  due  to  the  difference  in  the  basis  sets,  and 
partially  due  to  nonadiabatic  effects.  The  calculation  of  Qieung  et  al.^^  also  includes 
nonadiabatic  effects.  Their  cross  section  is  not  included  in  Table  I  because  their  ooi  is  not 
tabulated  in  Ref.  1 3.  There  is  no  reponed  nonadiabatic  calculations  in  the  literature  for  N2. 
However,  N2  is  a  heavier  molecule.  It  is  possible  that,  at  least  for  the  fundamental,  the 
nonadiabatic  effect  in  N2  is  smaller  than  in  H2.  Indeed  the  theoretical  cross  section  of  the 
N2  fundamental  agrees  well  with  experiment.  It  should  be  noted  that  the  N2  cross  section 
is  deduced  from  a  static  field  calculation.  However,  ooi  and  Yqi  of  N2  are  weakly 
dependent  on  Vp  at  488  nm,  and  only  a  small  error  is  expected  to  result  from  the  neglect  of 
frequency  dependence  in  these  two  quantities. 

For  the  overtone  transitions,  the  calculated  cross  section  of  H2  is  approximately 
20%  smaller  than  the  experimental  value,  but  lies  inside  the  experimental  error  limit  Fot 
N2,  the  theoretical  value  is  also  smaller  than  our  experiment  by  a  factor  of  two.  There  are 
two  possible  sources  for  the  discrepancies:  the  quality  of  the  electronic  wavefiinctions 
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used,  and  incomplete  treatment  of  nonadiabatic  effects.  In  the  simple  harmonic  oscillator 
approximation,  aw*  and  Yw’  for  the  first  overtone  transition  are  expressed  in  terms  of  the 
second  derivadves  of  the  dynamic  polarizabilides  a(Vp,  R)  and  yfVp,  R).  Thus,  a 
calculation  of  these  quantitites  generally  requires  wavefunctions  of  higher  quality  than  the 
fundamental  transition.  Even  though  exact  vibrational  wavefunctions  are  used  in  our 
calculation,  the  stringent  demand  of  high-quality  wavefunctions  probably  still  holds. 
Certainly  it  is  one  possible  source  of  the  observed  deviation. 

Nonadiabatic  effects  is  another  possible  source  of  error  in  the  theoretical  treatment 
of  the  N2  overtone  cross  section.  To  understand  the  importance  of  vibronic  coupling  in 
overtone  transitions,  it  is  worthwhile  to  analyze  the  two  nonadiabatic  calculatimis  of  H2. 
For  the  fundamental,  91%  of  Ooi  and  1 12%  of  Yoi  comes  from  the  contributions  of  B  and 
C  states.  (For  Yoi,  the  contribution  from  the  remainder  of  the  electronic  states  are  opposite 
in  sign  from  the  B  and  C  contribution,  resulting  in  a  final  value  smaller  than  the  B  and  C 
contribution  alone.)  It  is  reasonable  to  assume  that  most  of  tiie  nonadiabatic  effects  are 
accounted  for  if  Eq.  (4)  is  used  for  those  two  states.  On  the  other  hand,  only  72%  of  002 
and  34%  of  Y02  come  from  the  B  and  C  contribution.  Thus,  noindiabatic  treatment  of 
those  two  states  alone  is  probably  insufficient  It  is  even  more  illustrative  to  compare  the 
contributions  from  the  bound  and  continuum  vibrational  levels  of  those  two  states.  For 
C(oi  the  contribution  from  the  continuum  portion  of  B  state  vibrational  manifold  is  4.1%  of 
the  bound  level  contribution,  and  it  is  14.5%  for  the  C  state.  In  the  case  of  002,  the  ratios 
are  77.4%  and  74.2%,  respectively.  The  heavy  participation  of  the  B  and  C  continuum  in 
the  overtone  is  a  reflection  of  the  vastly  different  spatial  extent  of  the  v=0  and  v=2 
vibrational  wavefunctions  of  the  ground  state,  much  more  so  than  the  case  between  v=0 
and  1.  This  finding  strongly  suggests  that  contributions  from  other  electronic  states  should 
also  be  calculated  using  Eq.  (4)  instead  of  Eq.  (5). 
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For  N2,  nonadiabatic  studies  have  not  been  reported  in  the  literature.  However,  a 
comparison  of  the  adiabatic  cross  section  with  experiment  indicates  that  nonadiabatic 
effects,  assumed  to  be  unimportant  for  the  fundamental,  may  well  play  a  nwre  important 
role  in  the  overtone  transition.  Another  possible  source  of  error,  resulting  from  the  use  of 
static-field  (X02  and  Yo2i  is  estimated  to  be  1-2%,  based  on  the  corresponding  frequency 
dependence  found  in  H2. 

Finally,  we  can  compare  the  theoretical  results  for  the  v=:l  polarizability  an  to 
measurements^<^3  of  the  change  in  the  index  of  refraction  between  the  ground  and  first 
vibrational  state  at  693  nm.  The  index  of  refraction  for  the  v=l  level  is  particularly 
important  for  calculating  phase  perturbations  (and  hence  the  beam  quality)  that  occur  in 
large  scale  Raman  amplifiers.24  From  Tables  n  and  HI,  we  obtain  the  value  for  Aa  =  ai  1- 
Ooo  =  7.08  X  10‘26  cm^.  The  theory  predicts  that  the  index  of  refraction  of  the  v  =  1 
vibrational  level  is  8.6%  larger  than  that  of  the  v  s  0  level.  This  can  be  compared  with  the 
previously  measured  value23  of  Aa  =  1.53  ±  0.22  x  10-25  cm3  (19%  larger  than  the  v  =  0 
level)  and  a  semi-empirical  calculation25  of  1.40  x  10-25  cm3,  nie  discrqrancy  is  larger  by 
more  than  the  estimated  error  of  <  5%  in  the  present  tiieoretical  calculation.  Since  the 
present  theory  predicts  Ooo  of  v  =  0  to  within  the  error  bars  of  the  index  of  refraction  data^® 
(<1  %),  otoi  for  the  v=0— >1  transition  to  within  the  OTor  bars  of  Raman  scattering  data,2 
(<5%),  and  092  for  the  v=0->2  transition  to  within  the  experimental  error  bar  (20%),  it  is 
hand  to  rationalize  such  a  large  discrepancy  betweoi  the  experiment  and  our  theory.  We 
therefore  conclude  that  additional  experimental  and  theoretical  research  is  required  to 
resolve  the  differences. 
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CONCLUSION 


In  summaiy,  experimental  and  theoretical  results  are  presented  for  the  fundamental 
and  overtone  Raman  cross  sections  in  H2,  D2,  and  N2.  The  H2  fundamental  measurement 
is  in  better  agreement  with  the  partially  nonadiabatic  calculation  than  with  the  completely 
adiabatic  calculation.  In  addition,  the  H2  overtone  measurement  is  in  satisfactory 
agreement  with  the  partial  nonadiabatic  calculation.  Fot  N2,  although  the  fundamental 
adiabatic  result  is  in  good  agreement  with  expoiment,  the  overtone  adiabatic  result  is  a 
factor  of  2  smaller  than  the  experiment.  These  results  suggest  that  the  partially  nonadiabatic 
treatment  describes  the  Raman  process  in  gas  phase  diatomic  naolecules  better  than  a  purely 
adiabatic  treatment.  To  further  test  this  conclusion,  purely  nonadiabatic  calculations  and 
more  accurate  experimental  measurements  would  both  be  helpful. 
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Uncertainty  estimated  at  a  factor  of  2-3. 
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*  1  a.u.  =  1.481845  X  10-25  cm: 


oo 

O 

rvi 

ON 

«n 

c*^ 

fn 

«M 

r* 

Ov 

oo 

VO 

o 

O 

o 

fO 

r* 

00 

r~; 

ri 

o 

o 

o 

O 

o 

o 

VO 

o 

oo 

•n 

o 

•n 

§ 

cs  CM  <n 

8  8  8  8  8 

*n  oo 
8  8 

r< 

O 

1 

1 

1 

1 

1 

1 

r 

1 

1 

1 

QO 

r- 

VO 

m 

cs 

00 

Ov 

VO 

•r» 

*n 

m 

<s 

o 

n 

<s 

<N 

(M 

r4 

«N 

O 

o 

o 

o 

o 

o 

o 

O 

o 

O 

ZJ 

VO 

rvi 

OO 

Ov 

pH 

«*!• 

• 

o 

oo 

«rj 

p 

o 

CM 

Ov 

PH 

VO 

o 

o 

O 

pH 

VO 

CM 

pH 

CM 

r* 

VO 

CM 

oo 

VO 

8 

3 

o 

3 

3 

cn 

o 

cn 

o 

O 

CM 

cn 

1* 

r 

r 

1* 

r 

r 

• 

» 

r 

r 

p 

v> 

oo 

oo 

oo 

oo 

ov 

Ov 

Ov 

oo 

tn 

cn 

<n 

fO 

cn 

cn 

CM 

pH 

pH 

pH 

1* 

r 

r 

r 

r 

I* 

1* 

m 

r 

o 

r 

o 

Ov 

•n 

o 

O 

• 

vri 

•>4 

•Vt 

0*0 

q 

* 

00 

• 

ov 

• 

><*■ 

p 

ov 

• 

00 

hH 

§ 

S 

PH 

^H 

•p- 

■P” 

r>* 

CM 

3 

Ov 

VO 

o 

cn 

O 

oo 

ov 

p 

O 

pH 

q 

p 

Ov 

• 

p 

as 

p 

^H 

• 

pH 

* 

^H 

^H 

pH 

p*- 

oo 

Ov 

PH 

to 

cn 

o 

p 

p 

p 

p 

PH 

pH 

CM 

tn 

p 

^H 

PH 

pH 

pH 

^4 

pH 

pH 

pH 

PH 

•o 

p 

oo 

o 

oo 

CM 

Ov 

CM 

so 

O 

• 

Ov 

• 

o 

CM 

CM 

CM 

• 

od 

^H 

cn 

CM 

O 

o 

ri 

cn 

«d 

Ov 

1-* 

cn 

VO 

VO 

cn 

VO 

OV 

VO 

cn 

Ov 

Ov 

CM 

pH 

Ov 

cn 

Ov 

OO 

H 

»n 

•O 

p 

<o 

c^ 

r; 

ov 

pH 

o 

CM 

CM 

CM 

CM* 

CM* 

CM* 

CM* 

CM 

cn 

cn 

cn 

CM 

OV 

•o 

•n 

OV 

H 

Ov 

•o 

VO 

CM 

iC 

oo 

r4 

r- 

•n 

CM 

cn 

p 

o 

O 

pH 

hH 

CM 

cn 

*n 

oe 

•o 

*n 

vd 

vd 

vd 

vd 

vd 

vd 

vd 

vd 

1-21 


*  1 JLU.  =  1.481845  X  10-25  cm^ 


FIGURE  CAPTIONS 


Figure  1 .  Relative  frequency  dependence  of  the  square  of  the  average  polarizabilities, 

l«(v)/oe(v  =  0)|2,  for  Rayleigh  scattering  with  Av  =  0,  (solid),  fundamental 
Raman  with  Av  =  1  (medium  dash).  1st  overtone  Raman  vdth  Av  =  2  (short 
dash),  and  2nd  overtone  Raman  Av  =  3  (long  dash). 
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OPTICAL  STARK  SHIFT  SPECTROSCOPY: 
MEASUREMENT  OF  THE  V  =  1  POLARIZABILITY  OF  H2 
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Optical  Stark  shift  spectroscopy:  Measurement  of  the  u  =  1  polarizability  in  H2 

Mark  J.  Dyer  and  William  K.  Bischel* 

Molecular  Physics  Laboratory,  SRI  International,  Menlo  Park,  California  94025 
(Received  12  November  1990) 

We  report  a  quantitative  determination  of  the  u  =  1  polarizability  in  H2  by  measuring  the  optical  Stark 
shift  and  splitting  of  the  ^(0)  and  Qd)  vibrational  Raman  transitions.  For  the  Q{1)  transition,  the 
A/y  =  ±l,0  components  are  resolved,  thereby  allowing  a  determination  of  the  polarizability  anisotropy 
for  the  t;  =  l  state.  The  measured  optical  Stark  shift  of  the  QiO)  transition  at  1060  nm  is  23 
MHzGW"'  cm*,  giving  a  polarizability  difference  between  the  u  =1  and  0  states  of  Aa=0.49±0.03  a.u., 
in  good  agreement  with  ab  initio  theory. 


INTRODUCTION 

The  optical  Stark  effect  will  cause  the  rotational- 
vibrational  energy  levels  in  small  molecules  to  shift  and 
split  due  to  the  presence  of  a  high-intensity  nonresonant 
laser  held.  This  effect  has  recently  been  observed  for  both 
vibrational  [1]  and  rotational  [2]  Raman  transitions  in 
high-resolution  Raman  gain  experiments  using  an  opti¬ 
cally  applied  Stark  field.  Because  shifting  of  the  energy 
levels  effectively  broadens  the  observed  linewidths,  it  is 
one  of  the  important  mechanisms  that  limit  the  spectral 
resolution  of  coherent  Raman  experiments.  Hence,  the 
optical  Stark  effect  has  recently  been  included  in  the 
theory  [3]  of  saturation  processes  in  coherent  anti-Stokes 
Raman  spectroscopy  (CARS).  In  other  experiments,  the 
broadening  and  shift  of  the  anti-Stokes  radiation  generat¬ 
ed  by  multiwave  Raman  mixing  has  been  interpreted  to 
be  due  to  the  optical  Stark  effect  [4,5].  While  all  these  ex¬ 
periments  have  tried  to  interpret  their  results  based  on 
order-of-magnitude  estimates  of  the  optical  Stark  effect, 
there  have  been  no  quantitative  experiments  that  have 
used  the  nonresonant  optical  Stark  effect  to  derive  funda¬ 
mental  constants  for  the  excited  vibrational  states  of 
small  molecules. 

We  report  an  experiment  to  use  optical  Stark  spectros¬ 
copy  for  quantitative  measurement  of  the  polarizability 
and  the  polarizability  anisotropy  of  the  i;=l  vibrational 
state  in  H2.  Hydrogen  is  a  particularly  interesting  first 
molecule  to  study  for  a  number  of  reasons.  First,  there 
are  accurate  ab  initio  calculations  [6]  of  polarizabilities  in 
H2  that  can  be  used  to  make  a  quantitative  comparison 
between  theory  and  experiment. 

Second,  accurate  knowledge  of  the  optical  Stark  shift 
in  H2  is  of  great  practical  importance  in  the  design  of  Ra¬ 
man  amplifiers  and  beam  combining  techniques  [7]  since 
H2  is  the  gas  most  often  used  for  stimulated  Raman  fre¬ 
quency  conversion  applications.  For  this  application,  it 
is  quite  important  to  know  the  polarizability  of  the  u=l 
in  H2  since  the  change  in  the  polarizability  between  the 
v=0  and  the  v=l  vibrational  states  (Aa)  will  determine 
the  phase  perturbations  (and  hence  the  beam  quality)  that 
occur  across  the  spatial  aperture  of  a  large  scale  Raman 
amplifier  [7].  The  exact  value  of  Aa  has  been  an  issue  for 
a  number  of  years  due  to  the  fact  that  there  is  a  factor- 


of-2  discrepancy  between  current  measurements  [8]  and 
ab  initio  theory  [6].  An  accurate  measurement  of  the  op¬ 
tical  Stark  effect  for  the  vibrational  Raman  transition  in 
H2  would  resolve  this  controversy. 

The  goal  of  our  experiments  was  to  measure  the  optical 
Stark  effect  for  the  Q(0)  and  Q(l)  vibrational  Raman  tran¬ 
sitions  in  H2  and  compare  the  measurements  to  ab  initio 
theory.  The  energy  levels  and  splittings  in  an  optical  field 
for  the  Q(l)  transition  are  illustrated  in  Fig.  1,  where  the 
magnitude  of  the  shifts  in  units  of  MHz/(GW  cm “^)  are 
calculated  from  the  data  presented  in  this  paper  and  Ref. 
[6].  Without  the  high-intensity  nonresonant  optical  field, 
the  zero-pressure  Raman  frequency  shift  for  the  01) 
transition  is  4155.247±0.007  cm“'  as  measured  by 
Fourier-transform  Raman  spectroscopy  [9].  When  the 
nonresonant  field  is  applied,  there  is  a  shift  of  the  v=0 
and  1  vibrational  levels  from  their  zero  field  values,  and 


Zero  FMd  Raman  WHh  Nonraaonam 
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FIG.  1.  Schematic  illustrating  level  shifts  of  the  01)  branch 
in  H2  under  the  influence  of  a  strong  electric  field. 
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the  Mj  degeneracy  of  the  J=  1  rotational  level  is  split. 
(The  y=0  levels  will  only  show  a  shift  since  they  are  non¬ 
degenerate.)  As  illustrated  in  Fig.  1,  the  magnitude  of 
shift  in  the  vibrational  levels  (which  is  proportional  to  the 
average  polarizability;  see  theory  section)  is  approximate¬ 
ly  six  times  larger  than  the  magnitude  of  the  Mj  splitting 
(which  is  proportional  to  the  polarizability  anisotropy). 
We  can  indirectly  observe  these  energy-level  shifts  by 
carefully  measuring  the  shift  of  the  Raman  transition  fre¬ 
quency  (illustrated  by  the  arrows  in  Fig.  1)  with  and 
without  the  nonresonant  held.  The  shift  of  the  Raman 
frequency  is  due  to  the  difference  in  the  optical  Stark 
shifts  of  the  t;=0  and  1  levels,  and  hence  the  Raman  fre¬ 
quency  would  not  change  if  the  polarizabilities  for  these 
levels  were  the  same  (i.e.,  if  Aa=0). 

In  our  experiments  we  use  linear  polarization  for  both 
the  Raman  pump-probe  field  and  the  nonresonant  field. 
Therefore  we  have  a  selection  rule  of  AM=0  for  the  Ra¬ 
man  transitions.  As  illustrated  in  Fig.  1,  we  expect  to  ob¬ 
serve  two  Raman  transitions  for  Mj=0  and  ±1.  The 
splitting  between  the  Mj=0  and  ±  1  Raman  transitions  is 
only  ~9  MHz/GWcm“^)  making  these  splittings  ex¬ 
tremely  difficult  to  observe  experimentally. 

EXPERIMENT 

The  experimental  apparatus  is  illustrated  in  Fig.  2  and 
consisted  of  a  quasi-cw  Raman  gain  spectrometer  [10]  us¬ 
ing  a  Molectron  MY-34  Nd:YAG  (where  YAG 
represents  yttrium  aluminum  garnet)  laser,  modified  and 
injection-locked  to  a  Lightwave  Electronics  S-lOO  injec¬ 
tion  seeding  system  [1 1]  as  the  pulsed  pump  laser  source, 
and  an  argon-ion  pumped  Coherent  CR699-29  Ring  Dye 
Laser  as  the  continuous  tunable  narrow-band  probe  laser. 
The  Nd:YAG  laser  typically  produces  several  hundred 


millijoules  of  transform-limited  Ught  at  1.06  ^m  in  a  20- 
ns  FWHM  (full  width  at  half  maximum)  duration,  which 
when  doubled  to  S32  nm  yields  a  14-ns  pulse  with  a  band¬ 
width  of  22  MHz.  The  dye  laser  was  operated  at  683  nm 
to  probe  both  the  ^0)  and  ^1)  transitions  in  H2  with  a 
few  hundred  milliwatts  and  a  time-averaged  bandwidth  of 
roughly  8  MHz. 

As  shown  in  Fig.  2,  both  pump  and  probe  beams  were 
spatially  filtered  immediately  after  their  sources,  expand¬ 
ed  and  collimated  with  telescopes,  then  combined  on  a 
683-nm  dichroic  mirror  before  being  focused  together  by 
a  40-cm  focal  length  lens  into  a  cell  containing  hydrogen 
gas.  The  cell  was  constructed  from  stainless  steel  and 
designed  for  high  pressure.  Both  ends  of  the  cylindrical 
cell  were  sealed  by  fused  silica  windows  with  triple  “V” 
antireflection  coatings  at  the  three  wavelengths  used 
here.  The  spot  diameters  at  the  input  lens  of  the  Raman 
cell  were  1.6  and  2.1  cm  for  pump  and  probe,  respective¬ 
ly,  producing  FWHM  beam  waists,  which  were  scanned 
in  profile  by  a  1-^m  pinhole,  of  25  fim  in  the  center  of  the 
cell.  Assuming  a  shift  coefficient  of  25  MHzGW~'  cm^, 
this  relatively  tight  focusing  of  the  pump-probe  beams 
meant  that  the  energy  of  the  532-nm  pump  had  to  be  held 
below  120  fiJ  for  a  14-ns  FWHM  pulsewidth  (an  intensity 
of  1  GW/cm^)  in  order  to  avoid  more  than  a  5%  shift  or 
broadening  contribution  on  a  500-MHz  Raman  linewidth 
by  the  532-nm  pump.  Intensities  at  or  below  this  level, 
however,  provided  adequate  Raman  gain  signals  at  the 
hydrogen  densities  investigated,  and  the  small  Rayleigh 
ranges  of  the  pump  and  probe  allowed  selection  of  a  tight 
confocal  parameter  for  the  ir  beam  that  would  push  the 
intensity  of  the  applied  optical  field  beyond  100 
GW/cm*,  well  above  the  intensity  expected  to  resolve  the 
shifted  from  unshifted  line  shapes. 

The  high-intensity  nonresonant  optical  field  was  sup 


FIG.  2.  Experimental  layout  for  the  measurement  of  the  optical  Stark  effect  in  Hj. 
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plied  by  the  residual  ir  fundamental  of  the  Nd:YAG 
laser.  The  beam  was  first  passed  through  a  Laser  Metrics 
Electro-Optic  Chopper  (EOC),  then  into  a  vacuum  spatial 
filter,  and  collimated  for  delivery  to  the  Raman  cell.  The 
chopper  served  dual  purposes  as  a  switch  and  a  polarizer 
for  the  ir  beam.  The  spatial  filter  under  vacuum  cleaned 
up  the  transverse  mode  quality  of  the  beam  and 
significantly  reduced  power  fluctuations  of  the  ir  due  to 
air  particulates  and  pinhole  breakdown.  This  ultimately 
proved  vital  in  reducing  broadening  effects  on  the  shifted 
line  shape  associated  with  ir  field  variations.  The  col¬ 
limated  ir  beam  was  then  passed  through  a  SO-cm  lens 
and  counterpropagated  against  the  pump-probe  beams 
after  being  directed  into  the  Raman  cell  by  an  ir  high 
reflector  (HR)  transparent  to  the  pump-probe  combina¬ 
tion.  The  waist  of  the  beam,  measured  to  be  75  (im  in  di¬ 
ameter  at  FWHM,  was  carefully  positioned  to  overlap 
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FIG.  3.  Spatial  and  temporal  overlap  for  optical  Stark-shift 
measurements,  (a)  1/e^  intensity  data  points  for  spatial  cross 
sections  of  the  ir  and  probe.  The  theoretical  curves  for  the  1/e* 
intensity  plots  for  the  ir  probe  and  pump  (dashed)  beams  are 
shown.  Digitized  oscilloscope  traces  (b)  depict  the  time  evolu¬ 
tion  of  the  ir,  pump,  and  the  Raman  gain  signal. 


the  waists  of  the  pump  and  probe  beams;  the  optical 
paths  between  laser  source  and  cell  center  for  the  1064- 
and  S32-nm  beams  were  made  equal  so  that  the  peaks  of 
both  would  overlap  in  time  as  well.  Figure  3(a)  shows  the 
spatial  profiles  and  Fig.  3(b)  shows  the  temporal  overlap 
achieved  in  the  experiment.  Upon  exiting  the  cell,  the  ir 
beam  was  picked  off  by  an  identical  ir  dichroic  reflector 
and  sent  into  a  power  meter. 

After  passage  through  the  Raman  cell,  the  S32-nm 
pump  was  separated  from  the  probe  by  a  dichroic,  mea¬ 
sured  and  integrated  by  a  pyroelectric  energy  detector 
and  boxcar  averager,  and  recorded  by  computer  for  nor¬ 
malization.  The  transmitted  probe  was  down-collimated, 
filtered  to  remove  residual  pump  light,  mildly  focused 
through  an  acousto-optic  modulator  (AOM)  operating  in 
first  order,  and  sent  into  an  EG&G  FND  1(X)  photodiode 
having  a  1-ns  rise-time  response.  The  resulting  gain  sig¬ 
nal  detected  by  the  photodiode  exhibited  the  form  of  the 
14-ns  pump  pulse  shape  superimposed  atop  the  200-^is- 
wide  chopped  CW  probe  background.  This  signal  was 
enhanced  through  a  S-SOO-MHz  bandpass,  2()-dB  gain 
Avantek  amplifier,  where  the  CW  probe  offset  was 
filtered  out,  and  sent  into  two  SRS  boxcar  integrators 
operating  with  2-ns  gate  widths  and  toggled  at  5  Hz,  half 
the  trigger  rate  of  the  pump  laser. 

An  example  of  the  Raman  gain  signal  is  given  in  Fig. 
3(b).  The  single-shot  trace  was  acquired  as  the  ring  laser 
swept  in  frequency  over  the  transition.  The  result  of  the 
optically  induced  Stark  shift  of  the  transition  away  from 
the  probe  frequency  is  manifested  here  as  a  dip  in  the 
gain  signal  at  peak,  ir  intensity.  All  boxcar  channels  were 
triggered  off  the  Q-switched  Nd;YAG  laser  oscillator 
cavity  emission  in  order  to  reduce  boxcar  gate  jitter  on 
the  signals  to  below  0.2  ns.  The  averaged  output  of  the 
integrators  were  then  recorded  by  computer.  This  entire 
arrangement  made  possible  the  controlled  application  of 
the  ir  field  to  the  Raman  interaction  region  on  altering 
pulses  of  the  pump  laser  as  the  frequency  of  the  probe 
laser  was  tuned  through  the  vibrational  Raman  transi¬ 
tion,  thus  facilitating  collection  of  shifted  and  unshifted 
spectral  data  on  the  same  scan. 


LINE-SHAPE  ANALYSIS 

Analysis  of  the  data  was  performed  on  a  DEC  VAX 
11/750  using  a  routine  designed  to  remove  the  Gaussian 
laser  linewidth  contribution  and  then  generate  a  fit  to  the 
experimental  data.  The  line  shapes  used  to  fit  the  unshift- 
ed  data  were  the  theoretically  predicted  Lorentzian 
profiles,  and  the  optical  Stark-shifted  data  should,  under 
conditions  of  a  homogeneous,  steady-state  field,  remain 
Lorentzian  with  a  shift  in  frequency  being  the  only  dis¬ 
tinction  from  the  optical  field-free  data.  The  ir  Stark  field 
applied  in  our  experiment,  however,  had  Gaussian  tem¬ 
poral  and  spatial  mode  distributions,  so  its  slightly  varied 
temporal  and  spatial  behavior  over  the  regions  sampled 
in  time  and  space  contributed  to  a  broadening  of  the  ideal 
line  shape.  As  a  rule,  our  fits  were  started  using  the  un- 
shifted  linewidth  and  allowing  the  peak  position  to  be  the 
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only  adjustable  parameter;  then  width  and  height  param¬ 
eters  were  freed  to  obtain  the  best  fit  to  the  side  of  the 
Lorentzian  furthest  from  the  unshifted  line  position. 
This  method  of  analysis  was  based  on  the  fact  that  the 
peak  position  is  theoretically  predicted  to  shift  linearly 
with  peak  field  intensity,  and  the  part  of  the  Lorentzian 
shifted  the  furtherst  corresponds  to  gain  in  the  region 
subjected  to  the  most  intense  part  of  the  ir  field;  hence,  if 
there  are  intensity  gradients  away  from  the  peak  intrud¬ 
ing  upon  the  Raman  gain  region,  the  line  shape  will 
broaden  in  a  direction  favoring  the  lower  intensities  or 
smaller  shifts.  As  an  alternative,  assuming  nothing  about 
the  temporal  and  spatial  distributions  of  the  ir  field,  the 
data  were  also  fit  in  entirety  by  floating  all  parameters, 
and  although  the  differences  in  the  fitting  techniques  were 
observable,  shift  results  agreed  to  within  the  experimental 
error  of  the  measurements. 

RESULTS 

Typical  data  are  given  in  Fig.  4  at  a  density  of  1.81 
amagat  of  Hj  and  a  laser  intensity  of  nearly  100 
GW/cm^.  The  observed  shifts  were  density  independent 
for  both  transitions  in  the  absence  of  four-wave  mixing 
for  counterpropagating  ir  and  pump-probe  geometries. 
Polarizations  of  the  pump  and  probe  were  identical  for  all 
measurements,  and  the  ir  polarization  was  initially  fixed 
parallel  to  the  pump-probe  beams.  Measurements  involv¬ 
ing  different  combinations  of  ir  polarization  and  propaga¬ 
tion  direction  were  also  made  but  the  data  will  not  be 
presented  in  this  paper.  Figures  4  and  S  show  the  raw 
data  and  fit  as  a  function  of  frequency  for  an  ir  laser  in¬ 
tensity  of  100  GW/cm^  and  Fig.  6  graphs  the  dependence 


FIG.  4.  Raman  gain  of  the  unshifted  and  shifted  ^0)  branch 
of  H2  at  1.8 1  amagat  and  300  K. 


FIG.  S.  Unshifted  and  shifted  line  shapes  from  Raman  gain 
on  the  ^1)  branch  of  H2  at  1.81  amagat  and  300  K. 


of  the  optical  Stark  shift  versus  ir  field  intensity  for  the 
Q(0)  and  Q(0)  transitions,  along  with  a  linear  least-squares 
fit  to  each  of  the  Mj  levels  undergoing  the  shift.  For  the 
00)  transition  (Fig.  6),  the  line  shifts  scaled  with  ir  field 
intensity  and  the  data  yielded  a  shift  coefficient  of  23 
MHz  GW”*  cm^.  Shifting  of  the  Q  (1 )  transition  (Fig.  5) 
was  measured  to  be  27  MHz  GW”'  cm^  for  Mj=0  and 
19  MHzGW”'  cm^  for  —  The  Mj  splittings  of  a 
rovibrational  level  have  been  resolved,  thus  allowing  the 
determination  of  the  polarization  anisotropy  of  the  excit¬ 
ed  state. 


FIG.  6.  Plot  of  the  optical  Stark  shifts  of  00)  and  01)  in  H2 
vs  peak  ir  intensity,  and  linear  least-squares  fits  to  individual  Mj 
levels. 
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TABLE  1.  Wavelength  dependence  of  the  optical  Stark  shift  for  the  Q(0)  and  Q{1)  vibrational  Raman 
transitions  in  H2.  1  atomic  unit  (a.u.)=  1.481 845  X  10~^^  cm^. _ 

Optical  Stark  shift 

Wavelength  Polarizabilities  (a.u.)  Auj  (MHz  GW“'cm^) 

A.  -  fid)  fid) 


(nm) 

Ooo 

Yoo 

Oil 

yii 

Aa 

Ay 

fi(0) 

0 

II 

II 

H- 

00 

5.469 

2.013 

5.923 

2.473 

0.454 

0.460 

21 

27 

18 

600 

5.584 

2.085 

6.062 

2.576 

0.478 

0.491 

22 

29 

19 

500 

5.637 

2.118 

6.125 

2.623 

0.488 

0.505 

23 

29 

20 

400 

5.736 

2.181 

6.249 

2.716 

0.513 

0.525 

24 

31 

21 

300 

5.965 

2.331 

6.529 

2.936 

0.564 

0.605 

27 

34 

23 

250 

6.218 

2.501 

6.845 

3.193 

0.627 

0.692 

29 

38 

25 

200 

6.757 

2.881 

7.536 

3.789 

0.780 

0.910 

37 

48 

31 

Uncertainties  in  the  final  fits  were  relatively  small  com¬ 
pared  to  the  experimental  uncertainty,  which  we  conser¬ 
vatively  estimate  to  be  less  than  10%,  arising  primarily 
from  beam  overlap  and  peak  intensity  measurement  er¬ 
rors. 


THEORY  AND  DISCUSSION 


The  optical  Stark  shift  of  a  molecular  rotational- 
vibrational  level  with  quantum  numbers  v,J,M  can  be 
calculated  from  the  following  expression  [2,12]: 


where  a  is  the  isotropic  part  of  the  polarizability 
[a„  =  |{a||+2ai)]  and  y  is  the  anisotropic  part  of  the  po¬ 
larizability  (y  =a||— a^),  Eq  is  the  ir  laser  field,  and 


2  JU  +  l)-3Af^ 

3  {2J+3){2J-l)  ' 


(2) 


Using  the  definition  for  the  intensity  |£ol^=8ir//c, 
where  /  is  the  peak  ir  laser  intensity,  we  then  calculate 
the  optical  Stark  shift  Avs  for  a  Q  branch  transition  to  be 
the  difference  in  shifts  between  the  upper  and  lower 
states. 


Ai;s{MHz)=46.8/ 


GW 


cm* 


[Aa+C(J,M)Ar]  , 


(3) 


where  we  have  introduced  the  notation  Aa=aii~aoo  to 
denote  the  difference  between  the  i;=l  and  the  v—0  po¬ 
larizabilities  (similarly  for  Ay ).  In  Eq.  (3),  Av,  is  in  units 
of  MHz,  the  polarizabilities  are  in  atomic  units  (1 
a.u.  =  1.481 845X 10“^^  cm^),  and  /  is  in  units  of 
GW/cm^.  The  calculated  optical  Stark  shifts  as  a  func¬ 
tion  of  wavelength  are  listed  in  Table  I,  where  the 
wavelength-dependent  polarizabilities  are  taken  from 
Ref.  [6].  We  see  from  Table  I  that  the  optical  Stark  shift 
around  A. =1060  nm  ranges  from  19  to  29 
MHz  GW~ '  cm^  for  the  different  Mj  components  with  an 
average  value  of  23  MHz. 

The  ab  initio  theoretical  values  from  Table  I  are  com¬ 
pared  to  the  experimental  measurements  in  Table  II, 


where  the  agreement  is  surprisingly  good  for  an  experi¬ 
ment  where  the  laser  intensity  had  to  be  measured. 
Analysis  of  the  fid)  shifts  gives  Aa=0.462±0.02  and 
Ay=0.450±0.02.  No  significant  differences  for  Aa  of 
J=1  and  0  were  found.  These  results  indicate  that  the 
u  =  1  excited-state  polarizabilities  can  be  accurately  calcu¬ 
lated  for  the  simplest  diatomic  molecule  (H2),  and  pro¬ 
vide  a  favorable  basis  to  judge  the  accuracy  of  the  calcu¬ 
lations  of  the  u=  1  polarizability  at  other  wavelengths  as 
well  as  that  of  higher-order  polarizabilities  [6].  We  can 
also  compare  our  results  to  those  of  a  previous  experi¬ 
ment  [13]  that  reported  a  shift  of  approximately  —0.02 
cm~‘  for  the  fid)  transition  under  a  field  of  about  30 
GW/cm^  which  gives  a  shift  coefficient  of  —20 
MHzGW'cml 

The  difference  in  the  index  of  refraction  between  i;=  1 
and  0  (proportional  only  to  Aa)  is  also  given  in  Table  II 
to  be  8.4%.  This  is  compared  with  the  only  other  experi¬ 
ments  [8]  that  have  measured  this  quantity  to  be  19%. 
We  conclude  that  the  previous  experiments  published  a 
An  that  was  too  large  by  over  a  factor  of  2.  The  practical 


FIG.  7.  Index  of  refraction  of  Hj  (w=0)  and  (i>=  1)  vs  wave¬ 
length.  Experimental  data  for  u=0  taken  from  Ref.  [14]. 
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TABLE  II.  Index  of  refraction  measurement  of  Hj  (i>=  1)  at  k=  1060  nm  using  the  Stark  effect.  Er- 
ror  bars  are  one  standard  deviation  for  the  fit  illustrated  in  Fig.  6. _ 


Transition 

Av  (MHz/GWcm~*) 

A(n-l).,  (%) 

Expt. 

Theory 

Expt. 

Theory 

Ref.  [8] 

00) 

23.0±1.2 

22 

01)  Mj=Q 

27.0±0.7 

28 

8.4±0.3 

8.4 

19 

01)  Mj  =  ±\ 

19.0±0.6 

19 

consequence  of  this  measurement  is  that  the  beam  unifor¬ 
mity  of  the  pump  laser  required  for  large  scale  Raman 
amplifiers  [7]  can  be  significantly  relaxed  from  previously 
used  values,  thus  reducing  the  overall  cost  of  the  device. 

Finally,  we  illustrate  in  Fig.  7  the  comparison  between 
theory  [6]  and  experiment  for  the  absolute  index  of  re¬ 
fraction  for  the  i;=0  state  (where  the  data  have  been  tak¬ 
en  from  Ref.  [14])  and  the  t>=l  state  (data  from  the 
present  experiment).  Once  again  we  see  that  the  ground- 
state  polarizabiUties  are  well  known  and  can  be  used  to 
determine  the  absolute  polarizability  of  the  u=l  excited 
state. 

CONCLUSION 

We  report  a  quantitative  measurement  of  the  average 
polarizability  and  the  polarizability  anisotropy  of  the 


»=  1  excited  state  in  H2  from  a  measurement  of  the  non¬ 
resonant  optical  Stark  shift  for  the  00)  and  01)  Raman 
transitions.  In  these  experiments  the  Mj  splitting  of  the 
01)  line  has  been  observed,  thus  allowing  the  determina¬ 
tion  of  the  polarizabiUty  anisotropy.  The  results  are  in 
good  agreement  with  ab  initio  theory,  thus  giving 
confidence  in  the  ability  of  theory  to  predict  higher-order 
polarizabilities  in  H2. 
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The  temperature  and  density  dependence  of  the  rotational  Raman  linewidths  and  line  shifts  for 
the  diatomic  molecules  N;  and  H2  has  been  measured  using  stimulated  Raman  gain  spectroscopy. 
Room-temperature  results  for  the  density-broadening  coefficients,  B,  in  are  compared  with  ab  in¬ 
itio  calculations,  while  the  low-temperature  (80  and  195  K)  results  are  the  first  to  be  reported,  to  our 
knowledge.  Values  of  y  in  the  relation,  B  oc  7’’’,  were  determined  from  fits  to  the  data  for  self¬ 
broadening  in  N;  yielding  0.25  O'  <0.39.  Foreign-gas  (Oi)  broadening  coefficients  for  the  N^  tran¬ 
sitions  were  also  measured  and  found  to  be  10—15%  smaller  than  the  self-broadening  coefficients. 
Our  linewidth  and  line-shift  results  for  H2  are  in  general  agreement  with  previous  measurements. 


I.  INTRODUCTION 

Modeling  the  rotational  Raman  gain  in  Nj  is  of  interest 
since  stimulated  Raman  scattering  will  alter  the  mode 
quality  of  a  high-intensity  beam  propagating  through  air.' 
Accurate  knowledge  of  the  Raman  linewidths  is  neces¬ 
sary^  for  any  useful  model  because  the  steady-state  Ra¬ 
man  gain  is  inversely  proportional  to  the  linewidth.  Re¬ 
sults  of  previous  studies^"’  of  the  room-temperature  rota¬ 
tional  R^an  linewidths  in  N2  typically  differ  by  30%, 
while  results  at  other  temperatures  have  not  been  report¬ 
ed.  The  main  purpose  of  this  work  is  to  measure  the  tem¬ 
perature  dependence  of  the  linewidths  and  line  shifts  for 
the  rotational  Raman  transitions  (5  branch)  in  N2. 

This  work  differs  from  previous  measurements  of  rota¬ 
tional  Raman  linewidths^~^  in  three  aspects.  First,  we 
have  measured  the  temperature  dependence  of  the  line 
broadening  over  the  range  80—300  K.  Previous  work  in 
N2  has  been  only  for  room  temperature.  Second,  the  den¬ 
sities,  used  in  this  work  (0.01—2.0  amagat)  are  2  orders  of 
magnitude  lower  than  those  (1.0—100  amagat)  used  in 
previous  measurements.  Complications  due  to  overlap¬ 
ping  of  adjacent  lines  are  eliminated  by  working  at  lower 
densities  where  the  Raman  lines  are  typically  separated  by 
100  linewidths.  Finally,  all  previous  measurements  have 
used  spontaneous  Raman  'scattering,  whereas  we  have 
used  stimulated  Raman  scattering,  which  has  a  frequency 
resolution  several  orders  of  magnitude  higher  than  its 
spontaneous  counterpart. 

This  paper  reports  experimental  results  for  line¬ 
broadening  and  line-shift  coefficients  for  Stokes  rotational 
Raman  transitions  in  N2  and  H2.  Results  are  presented 
for  the  even-7  transitions,  5(2)  through  5(16),  at  tempera¬ 
tures  of  80,  19S,  and  295  K  in  N2  and  for  5(0)  and  5(1)  at 
temperatures  of  80  and  295  K  in  H2.  The  foreign-gas  (O2) 
broadening  and  shifts  of  the  four  strongest  N2  transitions 
have  also  been  measured  at  195  and  295  K.  The  room- 
temperature  broadening  coefficients  in  N2  are  compared 
with  ab  initio  calculations,  whereas  the  low-temperature 
N2  self-broadening  and  the  O2  foreign-gas  results  are  the 
first  to  be  reported  to  our  knowledge.  The  present  tem¬ 


perature  dependence  of  N2  is  compared  to  that  found  in 
studies  of  CO  and  CO2.  Finally  the  present  H2  results  are 
compared  with  previous  measurements. 

II.  EXPERIMENTAL  PROCEDURE 
A.  Apparatus 

An  overview  of  the  experimental  setup  is  shown  in  Fig. 
1.  This  quasi-cw  stimulated  Raman  spectrometer  is  simi¬ 
lar  to  that  developed  by  Esherick  and  Owyoung*  and  con¬ 
sists  of  a  cw  probe  laser,  a  tunable,  pulstri  pump  laser,  a 
pair  of  gas  (sample  and  reference)  c^s,  and  fast  photo¬ 
diodes  to  detect  the  induced  Raman  gain  on  the  probe 
beams. 

The  probe  laser  is  a  Kr'''-ion  laser,  operating  at  568  nm 
and  forced  into  single-mode  operation  with  the  insertion 
of  a  temperature-stabilized  etalon  into  the  cavity.  This 
laser  has  a  tirne-.iveraged  linewidth  less  than  30  MHz. 
The  probe  beam  «'as  split  into  two  equal  intensity  beams 
so  that  Raman  signals  could  be  simultaneously  obtained 
from  two  separate  gas  samples.  After  passing  through 
their  respective  sample  cells,  both  of  the  cw  probe  beams 
were  chopped  at  10  Hz  with  a  mechanical  chopper  (not 
shown  in  Fig.  1).  This  resulted  in  170-/;'"c  square-wave 
pulses  that  were  incident  on  the  photodi  v  es.  Chopping 
of  the  probe  increases  the  saturation  intensity  for  the  pho¬ 
todiodes.  The  average  power  of  the  probe  beams  was  150 
mW  in  both  the  sample  and  reference  cells. 

The  pump  laser  was  a  tunable,  single-mode  (1-MHz 
linewidth)  cw  dye  laser  that  was  pulse  amplified  using  a 
Quanta-Ray  PDA-1.  The  cw  dye  laser  input  to  the 
PDA-1  was  typically  400  mW.  The  PDA-1  was  pumped 
with  the  532-nm  output  (130  mJ  per  pulse)  from  a 
frequency-doubled  YAG  laser  (where  YAG  denotes  yttri¬ 
um  aluminum  garnet).  At  565  nm,  the  tunable  output  of 
the  PDA-1  consisted  of  2-mJ,  10-nsec  [full  width  at  half 
maximum  (FWHM)]  pulses  at  10  Hz.  The  linewidth  of 
the  pump  laser  was  100  MHz,  about  a  factor  of  2  larger 
than  the  Fourier  transform  for  a  10-nsec  pulse.  Peak 
pump  intensities  at  the  focal  planes  were  estimated  to  be  5 
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GW/cm*  in  the  sample  cell  and  70%  of  that  in  the  refer¬ 
ence  cell. 

The  pump  probe  beams  were  mode-matched  and  then 
were  focused  and  crossed  in  the  gas  cells  using  40-cm  fo¬ 
cal  length  lenses.  The  crossing  angle  was  1.0°  in  both  of 
the  cells.  A  series  of  diaphrams  placed  around  the  probe 
beams  was  used  to  eliminate  scattered  pump  light.  Both 
probe  beams  were  monitored  with  reverse-biased  photo¬ 
diodes  that  had  rise  times  of  less  than  1.0  nsec,  and  were 
terminated  with  50-11  resistors.  The  total  probe  power  in¬ 
cident  on  the  photodiodes  during  the  170-/isec  probe-laser 
pulses  was  10  mW.  Capacitors  (300  pF)  were  used  to 
block  the  170-/xsec  pulses  and  to  pass  the  10-nsec  Raman 
signals.  These  signals  were  amplified  by  a  factor  of  100 
and  then  averaged  with  gated  integrators.  Output  time 
constants  for  the  integrators  and  total  scan  dmes  were 
usually  O.S  sec  and  4  min,  respectively.  Occasionally,  for 
weak  transitions,  time  constants  and  scan  times  were  in¬ 
creased  to  2.5  sec  and  10  min,  respectively.  After  the 
pump  beam  exited  the  reference  cell,  it  was  also  moni¬ 
tored  with  a  third  photodiode.  This  signal  was  used  to 
normalize  the  two  Raman  signals  for  slow-pump  power 
variations  during  the  scans. 

The  sample  cell  was  constructed  by  attaching  a  pair  of 
3.0-cm-diam  Brewster  window  extensions  to  a  drilled-out 
10  cm^  cube  of  aluminum.  One  face  of  the  aluminum 
cube  is  in  contact  with  a  second  reservoir  that  can  be 
filled  with  either  liquid  N2  or  an  acetone  dry-ice  bath. 
The  cooled  sample  cell  was  insulated  with  a  vacuum,  and 
the  temperature  was  monitored  with  a  thermocouple  that 
was  imbedded  in  the  aluminum  cube.  The  reference  gas 
cell  was  kept  at  room  temperature  and  a  pressure  of  10 
Torr  for  N2  scans  and  200  Totr  for  most  of  the  K2  scans. 

The  temperature  dependence  of  the  linewidths  and  line 
shifts  was  measured  by  taking  measurements  at  three  ton- 
peratures:  80,  195,  and  295  K.  At  each  temperature  the 
Une  shape  was  recorded  for  several  densities  over  the 
range  0.01—2.0  amagat.  For  the  low- temperature  runs, 
there  was  a  time  period  ranging  from  a  few  to  15  minutes 
between  scans  with  different  densities,  so  that  any  gas 
added  between  scans  always  had  sufficient  time  to  come 
to  thermal  equilibrium  with  the  cell.  For  each  set  of  runs 
at  constant  temperature,  the  density  was  randomly 
cnanged  rather  than  constantly  adding  or  eliminating  gas 
from  the  sample  cell  for  each  succeeding  run. 

The  pressure  in  the  experimental  cells  was  measured  us¬ 
ing  a  lOOOO-Torr  MKS  Baratron  gauge,  and  the  tempera¬ 
ture  was  measured  using  a  Chromel-Alumel  thermocou- 
pled  gauge  referenced  to  273  K.  The  density  for  both  N2 
and  H2  was  calculated  using  the  perfect  gas  law  since  the 
density  calculated  using  the  second  virial  coefficient’  devi¬ 
ates  from  this  calculation  by  less  than  1%  at  the  densities 
and  temperatures  used  in  these  experiments. 


B.  Data  analysis 

All  of  the  experimental  line  shapes  were  fit  to  either 
Voigt  or  Lorentzian  profiles  rather  than  the  more  accu¬ 
rate  Galatry  profile, which  accounts  for  Dicke  narrow¬ 
ing.  Using  the  Voigt  profile  is  a  good  approximation  for 
the  N2  rotational  transitions  studied  here  since  the 


FIG.  1.  Schematic  of  experimental  apparatus  for  quasi-cw 
stimulated  Raman  spectroscopy. 


Doppler  broadening  (5  MHz),  and  hence  the  Dicke  nar¬ 
rowing,  is  negligible.  This  is  also  a  good  approximation 
for  the  H2  transitions,  since  the  instrument^  width  typi¬ 
cally  accounts  for  one-half  the  total  Unewidth  and  uncer¬ 
tainties  in  the  instrumental  Unewidth  will  mask  the  small 
difference  between  the  Voigt  and  Galatry  profiles.  Thus 
we  have  used  Voigt  profiles  to  fit  our  data  using  the  pro¬ 
cedure  described  here. 

Since  the  laser  Unewidth  was  comparable  to  the  low- 
density  Raman  Unewidth,  it  was  important  to  account  for 
this  instrumental  width  in  the  data  analysis.  Although 
the  pump-laser  frequency  proBle  was  not  directly  mea¬ 
sured,  it  was  indirectly  determined  by  scanning  the  laser 
over  a  narrow  Raman  Une  [e.g.,  S(0)  in  H2  at  5  Torr]. 
This  transition  has  negUgible  (1  MHz)  coUisional  broaden¬ 
ing  compared  to  100  MHz  of  Doppler  broadening.  The 
proHle  for  this  transition  was  found  to  fit  extremely  well 
to  a  pure  Gaussian.  Since  the  convolution  of  two  Gauss- 
ians  results  in  another  Gaussian,  we  concluded  that  the 
pump-laser  profile  could  be  approximated  as  a  Gaussian. 
Fitting  Voigt  profiles  to  low-density,  narrow  N2  Une  pro¬ 
files  with  known  Doppler  and  Lorentzian  contributions 
confirmed  this  conclusion.  The  probe-laser  Unewidth  was 
assumed  to  be  negUgible  when  convolved  with  the  pump 
Unewidth;  thus  all  the  data  presented  here  were  analyzed 


RELATIVE  LASER  FREQUENCY  DIFFERENCE  (MHil 

FIG.  2.  Example  of  S(10)  line  shape  in  Nj  at  295  K  and  10 
Torr  (narrower  profile)  and  at  195  K  and  400  Torr  (wider  pro¬ 
file).  Individual  dots  are  data  and  the  solid  lines  are  Voigt  fits. 
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FIG.  3.  SdO)  linewidths  in  N%  at  temperatures  of  80  K  (□), 
19S  K  (O ),  and  295  K  (A )■  Solid  lines  are  linear  least-squares 
fits  to  the  data. 

assuming  a  Gaussian  profile  for  the  instrumental  width. 

The  fitting  procedure  used  to  determine  the  Lorentzian 
portion  of  the  line  shapes  was  different  for  the  N2  and  H2 
data.  We  first  descri^  the  fitting  procedure  for  the  N2 
transition.  For  high  densities  in  N2>  Doppler  broadening 
(3  MHz)  and  the  laser  linewidth  (100  MHz)  can  be 
neglected.  This  fact  is  conflrmed  with  the  following  re¬ 
sult;  Fitting  Voigt  profiles  with  freely  varying  Lorentzian 
and  Gaussian  parameters  to  our  higher  density  (greater 
than  0.1  amagat)  data  gave  negligible  Gaussian  com¬ 
ponents.  Thus  sample  cell  line  shapes  for  high  densities 
were  first  fitted  to  pure  Lorentzians  to  get  a  first  approxi¬ 
mation  to  the  Lorentzian  width.  The  low-density  refer¬ 
ence  cell  data  were  then  fitted  to  Voigt  line  shapes  with  a 
fixed  Lorentzian  (determined  from  the  high-density  fits) 
and  variable  Gaussian  components.  This  Gaussian  com¬ 
ponent,  representing  the  pump-laser  profile  for  that  par¬ 
ticular  scan,  was  then  used  in  a  final  Voigt  fit  to  the  sam¬ 
ple  cell  profile.  Any  fluctuations  in  the  pump-laser 
linewidth  from  scan  to  scan  were  accounted  for  with  this 
procedure.  These  Gaussian  comfionents  were  typically 
within  ±20  MHz  of  an  approximate  average  value  of  100 
MHz.  The  change  in  Lorentzian  linewidth  determined 
from  the  first  Lorentzian  fit  and  the  final  Voigt  fit  was 
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FIG.  4.  J  dependence  of  the  N;  rotational  Raman  broadening 
coefficients  for  80  (□),  195  ( O ),  and  295  K  ( A ). 


TABLE  I.  Temperature  dependence  of  the  self-broadening 
coefficients  for  the  rotational  Raman  lines  of  N;. 


SU) 

295  K 

B  (MHz/amagat) 
195  K 

80  K 

0 

4730±700 

3140±470 

2 

4160±400 

3560±350 

2700±210 

4 

3580±60 

3230±320 

2490±250 

6 

3560±30 

3070 ±30 

2520±40 

8 

3270±60 

2860±30 

2120±60 

10 

3060±60 

2660+30 

1940±40 

12 

2870±50 

2340±260 

1690±100 

14 

2660±270 

2150±215 

16 

1840±185 

typically  a  few  percent  and  further  iterations  were  not 
necessary. 

For  the  H2  transitions  studied,  experimental  constraints 
did  not  allow  the  density  in  the  cell  to  be  larger  than  3 
amagat.  Therefore  the  approximate  density-broadening 
coefficient  could  not  be  determined  from  the  high-density 
data.  Hence,  the  above  procedure  was  not  applicable.  In¬ 
stead,  all  of  the  H2  line  shapes  were  fitted  to  Voigt  pro¬ 
files  with  freely  adjustable  Lorentzian  and  Gaussian  com¬ 
ponents. 

m.  RESULTS 
A.  Line  broadening  in  N2 

An  example  of  the  5(10)  Raman  profile  in  N2  is  shown 
in  Fig.  2.  The  narrower  profile  was  taken  at  10  Torr  and 
295  K  and  the  wider  profile  was  taken  at  400  Torr  and 
193  K.  Individual  dots  represent  the  data,  and  the  solid 
line  represents  the  fit  for  a  Voigt  line  shape,  which  is 
dominated  by  its  Lorentzian  component  for  the  400-Torr 
data.  A  slight  asymmetry  is  observable  in  the  data  of  Fig. 
2  and  was  characteristic  of  most  of  the  N2  data  obtained. 
An  asymmetry  of  this  magnitude  will  not  have  an  appre¬ 
ciable  effect  on  the  linewidth  determinations,  but  could 
produce  significant  systematic  errors  in  the  line-shift 
determinations  -because  the  line  shifts  are  so  small.  This 
asymmetry  is  discussed  later  in  Sec.  Ill  C. 

For  densities  larger  than  0.(X)4  amagat  ( —  3  Torr  at  298 
K),  the  rotational  Raman  lines  of  N2  are  collision 
broadened  giving  a  Lorentzian  line-shape  function. 
Hence,  the  linewidth  will  be  a  linear  function  of  density. 
We  have  therefore  taken  the  Lorentzian  component  result¬ 
ing  from  the  data-fming  procedure  described  in  Sec.  II B 

TABLE  II.  Comparison  of  different  measurements  of  the 
room-temperature  self-broadening  in  N2. 

B  (MHz/amagat) 


S(/) 

Present  work 

Ref.  3 

Ref  6 

Ref  7 

2 

4160±400 

3370 

4 

3580±60 

3050 

2590 

3280±300 

6 

3570±30 

3050 

8 

3270±60 

2660 

2300 

31501300 

10 

3070±60 

2590 

12 

2870±50 

2400 

1980 

28201300 

14 

2660±270 

2205 

16 

2080 
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LOWER  ROTATIONAL  QUANTUM  NUMBER  (Jl 

FIG.  S.  Comparison  of  theoretical  calculations  with  experi¬ 
mental  data  for  collisional  broadening  of  the  rotational  Raman 
lines  in  N}  at  room  temperature.  The  calculations  are  from  Ref. 
11  (long  dash).  Ref.  12  (short  dash),  and  Ref.  13  (solid  line). 
The  experimental  data  are  from  Ref.  3(D)  and  this  work  (A). 

and  determined  the  density-broadening  coefficients  B  for 
each  temperature  studied  by  fitting  the  linewidths  to  a 
straight  line.  The  data  and  the  linear  fit  (solid  line)  for 
the  5(  10)  transition  are  given  in  Fig.  3  for  temperatures  of 
298,  19S,  and  80  K.  Individual  points  indicate  the  data, 
and  the  solid  lines  represent  the  least-squares  fits.  The 
statistical  uncertainties  (one  standard  deviation)  for  the 
broadening  coefficients  determined  from  the  fits  in  Fig.  3 
are  1—2  %.  Systematic  errors,  due  to  uncertainty  in  the 
dye-laser  scan  width,  are  about  1%  or  less. 

Table  I  summarizes  aU  of  the  Nj  self-broadening  coeffi¬ 
cients  measured  in  this  study.  These  values  are  graphical¬ 
ly  illustrated  in  Fig.  4  as  a  function  of  rotational  quantum 
number.  Odd-/  transitions  in  N2  were  not  studied.  For 
5(4)  through  5(12),  the  specified  uncertainties  are  one 
standard  deviation  for  the  least-squares  fit.  For  5(0),  5(2), 
5(14),  and  5(16),  data  were  obtained  at  only  one  density, 
and  the  uncertainties  shown  for  these  transitions  represent 
the  error  in  the  linewidth  determination  of  a  single-line 
profile. 

Earlier  room-temperature  measurements  by  different 
authors^’*’^  have  shown  considerable  disagreement  (30%). 
Table  II  compares  these  earlier  results  with  our  present  re¬ 
sults.  It  can  be  seen  that  the  three  earlier  studies  show 
significant  differences,  while  our  current  measurement  is 
consistent  with  that  of  Ref.  7.  The  results  of  the  three 
previous  studies  that  are  shown  in  Table  II  were  obtained 
at  densities  where  adjacent  lines  appreciably  overlapped 
with  one  another,  whereas  the  present  results  were  ob¬ 
tained  at  densities  where  there  was  no  overlapping.  We 
postulate  that  the  differences  in  the  experimental  results 
of  the  previous  studies  given  in  Table  II  are  due  to  sys¬ 
tematic  errors  arising  from  the  overlapping  of  adjacent 
lines. 

Theoretical  calculations  of  the  collision-broadened  Ra¬ 
man  linewidth  have  been  perfomed  by  Grey  and  Van 
Kranendonk,"  Srivastava  and  Zaidi,'^  and  Robert  and 
Bonamy.'^  The  difference  in  these  three  approaches  is  in 
the  treatment  of  the  short-range  intermolecular  inter¬ 
actions.  The  results  of  these  different  approaches  are 


FIG.  6.  Lorentzian  (□)  and  Gaussian  (A)  components  from 
Voigt  fits  to  data  as  a  function  of  density  for  5(1)  in  H:.  Solid 
line  is  a  linear  one-parameter  fit  to  Eq.  (1). 

compared  with  our  measurements  in  Fig.  S,  where  the 
room-temperature  broadening  coefficients  are  shown  as  a 
function  of  rotational  quantum  number  /.  We  also  show 
the  experimental  values  of  Jammu  et  al?  for  comparison. 
It  should  be  noted  that  Srivastava  and  Zaidi'^  have  used 
an  adjustable  parameter  to  obtain  their  fit  to  the  data  of 
Jammu  et  al.,^  whereas  the  calculations  of  Grey  and  Van 
Kranendonk"  and  Robert  and  Bonamy"  do  not  use  any 
adjustable  parameters.  Figure  S  shows  the  ab  initio  calcu¬ 
lations  of  R(d)ert  and  Bonamy"  to  be  in  the  best  agree¬ 
ment  with  our  measurements. 

To  compare  the  pure  rotational  linewidths  to  the  pure 
vibrational  linewidths,  we  first  note  that  the  present  rota¬ 
tional  broadening  coefficients  are  10—15%  larger  than 
previously  measured**"**  Q-branch  vibrational  broaden¬ 
ing  coefficients  for  low  densities  (0.1  amagat).  Because 
negligible  Doppler  broadening  of  the  rotational  transitions 
means  that  Dicke  narrowing  is  also  negligible,  the  rota¬ 
tional  line  shapes  are  accurately  described  by  Voigt  pro¬ 
files,  while  the  vibrational  transitions  (with  20  times  more 
Doppler  broadening)  require  Galatry  profiles  for  the  most 
accurate  fits.  Also,  because  the  vibrational  transitions  (0.3 
cm"*)  are  more  closely  spaced  than  the  rotational  (8 
cm"*)  transitions,  coherent  line-mixing  effects  must  be  in¬ 
cluded  to  precisely  describe  the  vibrational  line  shapes  at 
densities  of  1  amagat  or  larger. 

Previous  to  this  work,  only  atomic  foreign-gas  (He  and 
Ar)  broadening  of  the  N2  rotational  transitions  had  been 
investigated.^  We  have  also  measured  diatomic  foreign- 
gas  (O2)  broadening  coefficients  at  195  and  295  K  for  the 
four  strongest  N2  transitions.  Our  O2  results  in  Table  III 


TABLE  III.  Temperature  dependence  of  the  foreign-gas  (O2) 
broadening  coefficients  of  N2  rotational  Raman  lines. 


SU) 

B  (MHz/amagat) 

295  K 

195  K 

6 

2970±130 

2690±120 

8 

2770±80 

2600±100 

10 

2750±50 

2410±100 

12 

2450±70 

I980±100 
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TABLE  IV.  Hi  rotational  Raman  parameters  for  an  ortho:para  ratio  of  3:1. 


T 

vr 

Do‘ 

A 

Pc 

(K) 

SU) 

(cm~') 

(MHz) 

(cm*  amagat  sec“  ‘) 

(MHz  amagat) 

(amagat) 

80 

0 

354 

48 

0.871 

1.37 

0.095 

1 

587 

80 

0.493 

2.10 

0.088 

295 

0 

354 

92 

1.19 

1.87 

0.068 

1 

587 

153 

1.42 

6.15 

0.134 

•Reference  18. 


show  that  the  foreign-gas  broadening  by  O;  is  typically 
85—90%  of  the  self-broadening  at  both  temperatures. 
Thus  the  foreign-gas  broadening  by  O2  is  approximately 
30%  larger  than  the  foreign-gas  broadening  by  the  atomic 
gases.^ 

B.  Line  broadening  in  Hi 

The  Raman  linewidth  for  H2  has  a  much  more  complex 
density  dependence  than  that  for  N2.  This  results  mainly 
from  the  fact  that  the  H2  density-broadening  coefficient  is 
approximately  a  factor  of  30  smaller  than  the  N2  broaden¬ 
ing  coefficient,  and  thus  Dicke'^  or  collisional  narrowing 
of  the  Raman  linewidth  can  be  observed**’*®  in  the  low- 
density  regime  if  the  resolution  is  high  enough.  A  sum¬ 
mary  of  the  relevant  collisional  processes  that  contribute 
to  the  linewidth  at  various  densities  is  given  in  Ref,  18. 

For  densities  above  a  certain  cutoff  density  (p,  defined 
below),  the  Raman  line-shape  function  is  Lorentzian  with 
a  FWHM  linewidth  Av  that  can  be  expressed  as 

Av=A/p-fBp ,  (1) 

where  A  (MHz  amagat)  is  a  coefficient  proportional  to 
the  self-diffusion  coefficient  Dq  (cm^  amagat  see**).  For 
forward  scattering,**  A  is  equal  to 

A  =4i7^i^Z)o  ^2) 

where  vj;  is  the  Raman  transition  frequency  in  cm~'. 

This  model  is  known  as  the  diffusion  model  of  the  Ra¬ 
man  Unewidth  and  it  diverges  as  the  density  approaches 
zero.  There  is  a  cutoff  density,  pc  at  which  this  model 
predicts  a  linewidth  that  is  10%  larger  than  that  predicted 
by  the  “hard  collision”  line-shape  theory.**  We  have  pre¬ 
viously  determined  that 

Pc  =  3.33i4/Av2)  ,  (3) 

where  Av/jCMHz)  is  the  FWHM  Doppler  width  of  the  Ra¬ 
man  transition  given  by  (forward  scattering) 

Ava=7. 15  X  10-‘'v/{ ( r/m  )*'^  ,  (4) 


where  vjj  is  in  MHz,  T  is  the  temperature  in  K,  and  m  is 
the  mass  in  amu.  For  densities  larger  than  p^,  the  line 
shape  is  Lorentzian  and  the  widths  can  be  effectively 
modeled  using  £q.  (1). 

For  example,  the  experimental  data  set  for  the  S(l) 
transition  at  298  K  is  given  in  Fig.  6.  From  the  Voigt  fit, 
we  determine  both  the  Lorentzian  and  Gaussian  contribu¬ 
tions  to  the  experimental  linewidth.  Since  the  Raman  line 
shape  is  Lorentzian  for  all  densities  shown  in  Fig.  6,  the 
Gaussian  contribution  is  due  to  the  spectral  width  of  the 
pulsed  dye  laser.  The  Lorentzian  and  Gaussian  contribu¬ 
tions  to  the  linewidth  are  plotted  in  Fig.  6,  where  we  see 
that  the  Gaussian  contribution  is  approximately  constant 
while  the  Lorentzian  contribution  varies  approximately 
linear  in  density. 

Unfortunately,  the  100-MHz  resolution  of  our  current 
laser  system  does  not  allow  the  Dicke  narrowing  of  the  H2 
rotational  Raman  linewidth  to  be  resolved  with  any  pre¬ 
cision.  We  have  therefore  chosen  a  fitting  procedure  that 
uses  Z>o  determined  from  our  previous  study**  to  calculate 
A  from  £q.  (2).  The  Lorentzian  component  of  the  experi¬ 
mental  linewidth  for  densities  larger  than  p,  is  then  fitted 
to  Eq.  (1)  with  one  free  parameter  S.  The  parameters 
used  in  this  fltting  procedure  are  given  in  Table  IV. 

The  density-broadening  coefficients  B  determined  using 
this  procedure  are  given  in  Table  V,  where  the  uncertain¬ 
ties  shown  are  one  standard  deviation  for  the  one- 
parameter  linear  least-squares  fit.  This  fit  is  illustrated  by 
the  solid  line  in  Fig.  6  for  the  5(1)  transition.  The  max¬ 
imum  correction  to  the  linewidth  due  to  including  the  col- 
lisional  narrowing  term  in  Eq.  (1)  is  on  the  order  of 
5-10%. 

We  have  also  included  in  Table  V  data  from  the  previ¬ 
ous  work  of  Van  Den  Hout  et  al}^  and  Cooper,  May,  and 
Gupta.^'  There  are  small  discrepancies  with  the  results  of 
Refs.  20  and  21.  Some  of  these  differences  could  be  due 
to  small  errors  in  our  estimate  of  A  for  these  rotational 
transitions.  We  emphasize  that  the  N2  results  of  Sec. 
HI  A  are  free  of  this  problem  since  Dicke  narrowing  is 
negligible. 


TABLE  V.  Self-broadening  coefficients  for  the  rotational  Raman  lines  of  Hj. 


SU) 

This  work 

293  K 

Ref.  20 

B  (MHz/amagat) 

Ref.  21 

80  K 

This  work 

Ref.  20 

0 

77±2 

84±2 

84±2 

67±2 

63±1 

1 

114±5 

104±2 

105±4 

110±3 

99±1 

K-5 
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C.  Line  shifts  in  N2  and  H2 

In  all  of  the  above  cases  for  which  density-broadening 
measurements  were  reported,  we  have  also  made  density 
line-shift  measurements.  In  all  cases  the  line  shifts  were 
only  a  few  percent  of  the  linewidth,  and  hence  any  asym¬ 
metries  in  the  line  shapes  need  to  be  carefully  analyzed. 

The  asymmetry  in  Fig.  2  is  typical  of  that  observed  in 
our  data  for  the  higher-density  Ni  data.  The  lower- 
density  N2  data  contained  an  asymmetry  that  occurred 
only  further  out  in  the  wing  of  the  spectral  profiles. 
About  80%  of  the  total  number  of  line  profiles  at  all  den¬ 
sities  and  temperatures  exhibited  this  small  asymmetry, 
where  the  data  were  slightly  larger  than  the  symmetric  fit 
on  the  high-frequency  side  of  the  line.  Roughly  20%  of 
the  line  profiles  were  symmetric  and  only  rarely  did  the 
asymmetry  appear  such  that  the  data  were  smaller  than 
the  fit  on  the  high-frequency  side.  It  was  determined  that 
this  asymmetry  was  not  due  to  the  saturation  of  the  detec¬ 
tion  system  by  either  the  probe  beam  or  the  Raman  signal. 
By  scanning  the  dye  laser  backwards,  it  was  determined 
that  the  asymmetry  always  remains  such  that  the  data  are 
larger  than  the  fit  on  the  high-frequency  side  of  the  line, 
where  the  high-frequency  side  corresponds  to  a  larger  Ra¬ 
man  shift.  For  densities  above  0.5  amagat  and  for  the 
pump-laser  intensities  used  here,  this  asymmetry  is  not 
consistent  with  the  measured  magnitudes  of  Stark  effects 
in  N2.^^  At  this  time  we  are  unable  to  determine  the 
cause  of  these  asymmetries. 

In  N2,  both  the  self-density  shifts  and  the  O2  foreign- 
gas  density  shifts  were  found  to  lie  in  the  range  50—150 
MHz/amagat  for  the  transitions  and  temperatures  dis¬ 
cussed  above.  There  was  no  uniform  dependence  with  ei¬ 
ther  temperature  or  quantum  number  J.  More  important¬ 
ly,  there  were  inconsistencies  observed  in  both  the  data  ac¬ 
quired  during  the  course  of  a  single  day  and  in  the  data 
collected  from  day  to  day.  These  inconsistencies  have 
forced  us  to  conclude  that  the  line-shift  measurements  in 
N2  are  substantially  disturbed  by  the  asymmetric  line 
shapes.  The  data  can,  however,  be  used  to  establish  an 
upper  limit  to  the  magnitude  of  the  line  shifts.  For  all 
transitions  and  temperatures  indicated  for  N2  in  Tables  I 
and  III,  the  density  line  shifts  are  less  than  ISO 
MHz/amagat.  Future  elimination  of  the  asymmetry 
would  enable  us  to  detect  line  shifts  as  small  as  10 
MHz/amagat  with  our  current  experimental  apparatus. 

The  densi^'y  shifts  observed  in  H2,  unlike  those  in  N2, 
were  free  the  inconsistencies  mentioned  above.  Data 
for  5(1)  in  H2  at  295  K  are  illustrated  in  Fig.  7,  where  the 
density  shift  is  defined  as 

8v=V2j(p)— v«(0)  ,  (5) 

where  is  the  Raman  transition  frequency  at  the  den¬ 
sity  p.  Hence  a  negative  shift  results  in  a  smaller  Raman 
frequency.  Density  shifts  were  plotted  as  a  function  of 
density  difference  between  the  sample  and  reference  cells 
rather  than  density  because  data  were  collected  for  refer¬ 
ence  cell  pressures  of  50,  1(X),  and  2(X)  Torr.  The  tem¬ 
perature  is  295  K  for  the  data  of  Fig.  7.  A  summary  of 
our  H2  density-shift  results  is  given  in  Table  VI,  along 
with  the  results  from  Refs.  20  and  21.  The  agreement  be- 


0  0.5  1.0  1.5  2.0  2.5 

DENSITY  (amagat) 

FIG.  7.  Example  of  data  for  line-shift  determination  for  5(1) 
in  H2  at  295  K.  The  solid  line  is  the  least-squares  fit. 

tween  our  results  and  those  of  Ref.  21  is  consistent  with 
the  combined  uncertainties  for  the  two  experiments. 

It  should  be  noted  that  the  density  shift  for  the  rota¬ 
tional  Raman  line  in  H2  is  substantially  different  from 
that  observed  for  the  vibrational  Raman  lines.’*  For  ex¬ 
ample,  at  80  K,  the  shift  of  the  Q(I)  vibrational  transition 
is  13  times  larger  than  the  5(1)  rotational  Raman  line. 
This  difference  results  from  the  perturbation  of  the  vibra¬ 
tional  frequency  of  the  H2  molecule  during  an  elastic  col¬ 
lision,  an  ^fect  not  present  for  the  rotational  transition. 

IV.  DISCUSSION 

A  simple  qualitative  understanding  of  this  temperature 
dependence  for  the  N2  rotational  Raman  broadening  data 
given  in  Table  I  can  be  obtained  from  a  consideration  of 
Anderson’s  impact  theory  of  the  coUisional  broadening.^ 
The  linewidth  can  be  expressed  as 

Av=n<uff(i;)>~n<i;><a(u)>  ,  (6) 

where  n  is  the  molecular  number  density  (cm~^),  u  is  the 
relative  velocity  of  the  collision  partners,  <7  is  the 
collisional-broadening  cross  section,  and  the  angular 
brackets  denote  an  average  over  the  velocity  distribution 
of  the  gas.  As  can  be  seen  from  Eq.  (6),  the  temperature 
dependence  of  the  density-broadening  coefficient  is  impli¬ 
cit  in  the  velocity  dependence  of  the  coUisional  cross  sec¬ 
tion.  For  example,  if  we  assume  a  “hard-sphere”  model 
for  the  cross  section  (velocity  independent),  the  density¬ 
broadening  coefficient  (Av/n)  measured  in  our  experi¬ 
ment  would  be  proportional  to  7’*’  *. 

The  broadening  cross  section  is  composed  of  contribu¬ 
tions  from  several  different  types  of  collisions.  Both  das- 
tic  phase-perturbing  collisions,  and  rotationally  inelastic 


TABLE  VI.  Temperature  dependence  of  the  density  line 
shifts  for  H2  rotational  Raman  lines. 


6v=Vj,(p)  — v*(0) 

(MHz/amagat) 

295  K 

80  K 

5(7) 

This  work  Ref.  21 

This  work 

Ref.  20 

0 

6.5±3  3.111.5 

-22.310.6 

-22.3 

1 

5.911  4.011.5 

-23.O1O.6 

-17.8 
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collisions  contribute  to  the  linewidth.  However,  there  are 
several  cases  where  the  contribution  from  elastic  collisions 
is  small.  For  example,  the  linewidth  for  isotropic  vibra¬ 
tional  2'branch  Raman  scattering  have  no  contribution 
from  elastic  collisions*'*  if  collisions  perturb  the  ground 
and  final  vibrational  level  the  same  amount.  This  appears 
to  be  a  good  approximation  for  all  cases  except  for  hydro¬ 
gen.'* 

Rotational  Raman  linewidths,  on  the  other  hand,  have 
contributions  from  phase-perturbing  elastic  collisions. 
This  is  the  primary  reason  why  rotational  Raman  transi¬ 
tions  usually  have  larger  linewidths  than  vibrational  Ra¬ 
man  transitions.  However,  for  the  case  of  N,,  CO,  and 
CO2,  it  has  been  calculated' '  that  elastic  collisions  contri¬ 
bute  no  more  than  15%  of  the  cross  section.  The  present 
N2  rotational  linewidth  measurements  are  15%  larger 
than  the  previous  vibrational  measurements,'^"'*  in 
good  agreement  with  the  conclusion  that  elastic  collisions 
account  for  15%  of  the  rotational  linewidths.  Therefore, 
in  most  cases,  the  primary  problem  is  to  calculate  the 
velocity  dependence  of  the  cross  section  for  rotationally 
inelastic  collisions. 

The  calculation  of  the  collisional  cross  section  has  been 
the  subject  of  considerable  research.  A  good  review  of  the 
various  approaches  to  this  calculation  for  Raman  transi¬ 
tions  has  been  given  by  Srivastava  and  Zaidi.^*  One  ap¬ 
proach  to  this  calculation  is  to  expand  the  long-range  part 
of  the  intermolecular  potential  in  a  multipole  expan- 
sion,"’^’^^  and  then  calculate  the  contribution  to  the 
cross  section  of  each  term  in  this  expansion.  Using  this 
approach.  Gray  and  Van  Kranendonk"  have  shown  that 
for  molecules  such  as  N2,  CO,  and  CO2,  the  most  impor¬ 
tant  terms  in  this  expansion  for  self-broadening  are  the 
quadrupole-quadrapole  (QQ)  potential  and  the  dispersion 
potentid. 

It  can  be  shown^^  in  general  that  the  velocity  depen¬ 
dence  of  these  contributions  to  the  cross  section  scales  as 

,  (7) 

where  n  is  the  exponent  of  the  radial  dependence  of  the 
long-range  part  of  the  potential  in  the  form  F'(r)~r"". 
For  the  QQ  potential,  n  =  5,  while  for  the  dispersion  po¬ 
tential  n=6.  Substituting  Eq.  (7)  into  Eq.  (6)  and 
remembering  that  vo:T°-^,  we  find  that  the  density- 
broadening  coefficient  should  scale  with  temperature  as 

.  (8) 

From  this  simple  argument,  we  would  expect  our 
density-broadening  coefficients  to  scale  as  de¬ 

pending  on  the  relative  contribution  to  the  cross  section  of 
the  QQ  and  dispersion  forces. 

To  test  this  simple  model  of  the  temperature  depen¬ 
dence,  we  have  fitted  our  data  in  Table  I  to  the  equation 

B=Bo(T/295V  .  (9) 

where  Bq  is  the  density-broadening  coefficient  at  295  K. 

The  results  are  given  in  Table  VII  as  a  function  of  initial 
rotational  state.  We  find  that  0.26 <y  < 0.39,  in  reason¬ 
able  agreement  with  this  simple  picture.  We  believe  that 
most  of  the  scatter  in  7  as  a  function  of  J  is  due  to  the 


TABLE  VII.  Fit  parameters  for  temperature  dependence  of 
density  self-broadening  coefficients  in  N;. 


Bo 

Co 

SU) 

(MHz  amagat”') 

r 

(MHz  amagat~'K“') 

2 

4160 

0.33+0.02 

6.8±0.4 

4 

3580 

0.28±0.01 

5.13:0.7 

6 

3560 

0.26±0.04 

bo 

If 

0 

8 

3270 

0.33±0.004 

If 

p 

10 

3060 

0.35±0.002 

5.2±0.6 

12 

2870 

0.39±0.03 

5.4±0.1 

fact  that  a  much  larger  temperature  range  needs  to  be  in¬ 
vestigated  to  accurately  determine  this  small  temperature 
dependence. 

However,  the  scaling  given  in  Eq.  (9)  may  not  be  valid 
over  a  wide  temperature  range.  By  theoretically  examin¬ 
ing  temperature  dependence  of  the  rotational  Raman 
linewidth  over  a  wide  range,  Pack**  proposes  a  quadratic 
temperature  scaling  law  for  density-broadening  coeffi¬ 
cients  that  can  be  written  as 

B  =Bo+Co(T -295)+Do{T -295)^  ,  (10) 

where  Bq  is  the  broadening  coefficient  at  298  K  [same  as 
in  Eq.  (9)],  and  Cq  and  Do  are  fitting  coefficients.  We 
have  previously  us^  this  temperature  scaling  formula  in 
modeling  Raman  linewidths  in  H2.'*  The  data  from 
Table  I  have  been  fitted  to  Eq.  (10)  assuming  Do—0,  and 
the  results  are  also  given  in  Table  VII.  We  find  that  this 
expression  is  just  as  useful  as  Eq.  (9)  in  modeling  the 
dependence  of  the  density-broadening  coefficient  over  the 
limited  temperature  range  investigated  in  our  experiment. 
Both  Eqs.  (9)  and  (10)  can  be  used  to  predict  B  to  better 
than  3%  over  the  range  of  80—300  K.  Recent  studies 
show  that  the  temperature  dependence  for  vibrational  N2 
line  shape  is  accurately  described  by  a  modified  exponen¬ 
tial  pp  scaling  law^^  or  a  polynomial  inverse  energy-gap 
law.^*  Thus,  although  Eqs.  (9)  and  (10)  are  useful  for 
predicting  low-temperature  rotational  Rxunan  linewidths 
for  N2,  additional  data  (more  data  points  covering  a  lar;ger 
temperature  region)  are  necessary  to  unambiguously  deter¬ 
mine  the  temperature  dependence. 

There  are  few  molecular  systems  for  which  the  rota¬ 
tional  Raman  linewidths  have  been  determined  as  a  func¬ 
tion  of  temperature,  and  no  other  studies  in  N2.  However, 
there  have  been  extensive  studies  of  the  temperature 
dependence  of  absorption  linewidths  in  CO  (Ref.  29)  and 
C02.^  The  absolute  magnitude  of  the  linewidths  deter¬ 
mined  from  absorption  will  be  different  from  those  deter¬ 
mined  from  Raman  scattering.  This  is  due  to  the  fact 
that  different  states  are  involved  in  the  calculation  of  the 
collision  cross  section.  However,  the  temperature  depen¬ 
dence  will  be  similar  if  the  same  terms  in  the  multipole 
potential  are  responsible  for  the  broadening.  This  is  the 
case"  for  both  CO  and  CO2  where  the  most  important 
terms  in  the  potential  are  the  same  as  for  N2. 

Most  of  the  data  for  linewidths  in  the  literature  are 
given  in  terms  of  pressure-broadening  coefficients.  We 
can  compare  results  by  noting  that  y  in  Eq.  (9)  is  related 
to  7j  by  y  =  1  —  7j,  where  rj  is  the  exponent  of  the  tempera- 
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ture  for  pressure-broadening  coefficients.  The  17  values 
for  the  Ni  rotational  Raman  linewidths  of  this  study  cov¬ 
er  a  similar  range  of  values  as  observed  for  the  N;  vibra¬ 
tional  Raman  linewidths’*  and  the  absorption  linewidths 
in  CO  (Ref.  29)  and  00;.*° 

The  details  of  the  exact  theoretical  temperature  depen¬ 
dence  of  the  rotational  Raman  linewidths  requires  an  ab 
initio  calculation  of  the  velocity-dependent  collisional 
cross  section,  and  is  beyond  the  scope  of  this  study.  How¬ 
ever,  the  precision  of  the  data  presented  here  should  allow 
a  detailed  comparison  with  linewidths  calculated  from  the 
various  theoretical  approaches  thus  allowing  new  theoreti¬ 
cal  advances  to  be  made. 

V.  CONCLUSIONS 

The  density  self-broadening  coefficients  have  been  mea¬ 
sured  at  temperatures  of  80,  195,  and  295  K  for  the  even- 
J  rotational  Raman  transitions  in  N2.  The  temperature 
dependence  of  these  coefficients  was  found  to  fit  both 
linear  and  7^  models  to  within  3%,  however,  additional 
temperature-dependent  data  are  necessary  to  accurately 
determine  the  temperature  dependence  of  these  rotational 


linewidths.  The  room-temperature  results  of  this  work 
are  in  agreement  (less  than  5%  differences  for  J>2)  with 
the  most  recent  ab  initio  calculations.  Foreign-gas 
broadening  of  these  Ni  transitions,  due  to  O2,  was  also 
measured  and  found  to  be  10—15%  smaller  than  the 
self-broadening  coefficients.  We  have  also  determined 
that  the  density  self-shifts  and  O2  foreign-gas  shifts  for 
these  same  N2  lines  are  no  larger  than  150  MHz/amagat. 
In  a  previous  paper,’  we  have  used  the  above  results  in  the 
modeling  of  the  temperature  dependence  of  the  steady- 
state  Raman  gain  coefficient  in  N2. 

The  rotational  Raman  density  broadening  and  shifts  in 
H2  were  measured  at  80  and  295  K.  Our  results  were  in 
general  agreement  with  previous  experimental  work.  This 
agreement  of  our  H2  data  with  other  independent  mea¬ 
surements  provides  confidence  in  the  new  temperature- 
dependent  N2  results  presented  above. 
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The  temperature,  wavelength,  and  J  dependence  of  the  rotational  Raman  gain  coefficient  has  been  determined  for 
the  S-branch  transitions  in  N2.  First,  the  temperature  dependence  (80-300  K)  of  the  density-broadening  coefH- 
cients  was  measured  using  stimulated  Raman  spectroscopy.  Second,  the  wavelength  dependence  (250-600  nm)  of 
the  polarizability  anisotropy  was  empirically  determined  by  fitting  to  the  most  recent  experimental  data.  These 
two  results  were  then  used  to  calculate  the  rotational  Raman  gain  coefficients  with  accuracies  estimated  at  5%.  At 
room  temperature,  the  high-density  limit  of  the  steady-state  gain  coefficient  of  the  strongest  line,  5(10),  is  4.8  X 
10~‘^  cm/W  for  a  Stokes  wavelength  of  568  nm. 


There  are  many  applications,  such  as  laser  fusion,  that 
require  the  propagation  of  high-intensity  lasers 
through  the  atmosphere.  Recently  there  has  been 
renewed  interest  in  the  nonlinear  optical  processes 
that  limit  the  maTimnm  transmitted  laser  intensity 
for  a  given  propagation  distance.  Several  different 
nonlinear  effects  may  be  important,  and  some  of  these 
were  summarized  by  Zuev^  in  1982.  However,  it  has 
only  recently  been  suggested  that  rotational  stimulat¬ 
ed  Raman  scattering  (RSRS)  in  N2  may  ultimately 
limit  the  maximum  intensity  that  can  be  propagated 
through  air  if  the  frequency  and  divergence  of  the  laser 
are  to  remain  unchanged.  This  suggestion  is  support¬ 
ed  by  recent  experiments^  at  Lawrence  Livermore  Na¬ 
tion^  Laboratory,  where  RSRS  in  N2  was  observed 
when  the  NOVA  fusion  laser  was  propagated  through 
100  m  of  air  at  intensities  in  the  range  of  a  few  giga¬ 
watts  per  square  centimeter. 

This  observation  has  reinforced  the  need  for  accu¬ 
rate  modeling  of  stimulated  Raman  scattering  in  air. 
Unfortimately,  the  existing  uncertainties^*^  in  the 
steady-state  gain  coefficient  in  N2  limit  the  usefulness 
of  any  potential  model.  In  particular,  no  data  exist  for 
the  temperature  dependence  of  the  gain  coefflcient,  an 
important  consideration  in  any  model.  In  this  Letter 
we  present  new  experimental  results  for  the  tempera¬ 
ture  dependence  of  the  Raman  gain  coefflcient  for  the 
most  important  S-branch  rotational  transitions  in  N2. 

Direct  measurement  of  the  gain  coefflcient  requires 
an  accurate  measurement  of  the  laser  intensity  in  the 
interaction  region.  For  low-gain  Raman  systems  such 
as  N2,  this  is  extremely  difficult  since  a  focused  geome¬ 
try  is  usually  required  if  adequate  signal  strengths  are 
to  be  obtained.  Alternatively,  it  was  recently  experi¬ 
mentally  demonstrated*  that  the  gain  coefflcient  can 
be  accurately  calculated  if  the  line-shape  function  for 
the  Raman  transition  and  the  polarizability  anisotro¬ 
py  are  known  precisely.  We  have  chosen  tUs  alterna¬ 
tive  method  for  determining  the  rotational  Raman 
gain  coefflcients  in  N2.  Using  stimulated  Raman  gain 
spectroscopy,^  we  have  measmed  the  density  and  tem¬ 
perature  dependence  of  the  Raman  linewidths  for  the 


Stokes-branch  rotational  transitions  in  N2.  In  addi¬ 
tion,  the  relative  scattering  cross  sections  were  mea¬ 
sured  for  these  same  rotational  lines.  Using  the  re¬ 
sults  of  these  measurements  and  previously  published 
values^^  of  the  polarizability  anisotropy,  we  have  cal¬ 
culated  the  high-density  Raman  gain  in  the  tempera¬ 
ture  region  80-300  K  for  four  of  these  rotational  tran¬ 
sitions. 

For  a  Lorentzian  line-shape  function,  the  peak 
plane-wave,  steady-state  Raman  gain  coefflcient  is^** 

*  hv,irAi>  da 

where  v,  is  the  Stokes  frequency  in  hertz,  X,  is  the 
Stokes  wavelength  in  the  m^ium,  AN  is  the  linewidth 
(FWHM)  in  hertz,  and  AN  is  the4Kq}ulation  differ¬ 
ence,  which  is  equ^  to  N{J )  —  [N(«7’)(2«7  +  1)/{2J'  + 
1)].  N(J)  is  the  Boltzmann  statistical  population 
density  in  a  state  with  rotational  quantum  nvunber  J. 
The  rotational  Raman  scattering  cross  section  for  inci¬ 
dent  and  scattered  waves  linearly  polarized  in  the 
same  direction  is  given  by 

^  ^  2  (J+l)(J  +  2)  2 

dfl  c  /  (2J+l)(2J  +  3)'*'  ’  ^ 

where  y  is  the  polarizability  anisotropy.  The 
linewidth  in  Eq.  (1)  varies  with  J,  density,  and  tem¬ 
perature,  while  7  varies  with  pump-laser  frequency. 
Thus  it  is  necessary  to  know  this  parametric  depen¬ 
dence  in  order  to  calculate  the  Ramw  gain  accurately. 

The  experimental  setup  is  a  quasi-cw  stimulate 
Raman  spectrometer  similar  to  that  in  Ref.  5.  The  cw 
probe  laser,  which  provides  150  mW  of  power  in  the 
interaction  region,  is  a  single-mode  Kr*’-ion  laser  oper¬ 
ating  at  568  nm.  The  pump  laser  is  a  Coherent  699-29 
single-mode  cw  ring  dye  laser  that  is  pulse  amplified 
using  a  Quanta-Ray  PDA-1  amplifler.  At  565  nm,  the 
tunable  output  of  the  pump  laser  consists  of  3-mJ,  10- 
nsec  (FWHM)  pulses  at  10  Hz.  The  frequency  resolu¬ 
tion  of  this  system,  limited  by  the  linewidth  of  the 
pulsed  laser,  is  approximately  100  MHz.  A  250-Mm- 
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Fig.  1.  Measured  temperature  dependence  of  the  density¬ 
broadening  coefficient  for  Stokes-branch  rotational  Raman 
transitions  in  N2.  The  solid  line  indicates  the  most  recent 
room-temperature  ab  initio  calculations  given  in  Ref.  12. 


diameter  pinhole  is  used  to  improve  the  spatial  mode 
of  the  pump  beam  and  to  verify  that  the  beam  does  not 
move  when  the  dye  laser  is  tuned  over  wavelength 
intervals  of  10  nm.  The  pump  and  probe  beams  were 
focused  and  crossed  (angle  of  1  deg)  at  the  center  of  a 
gas  sample  cell  using  a  40-cm  focal-length  lens.  After 
the  cell,  stray  pump  light  was  isolated  from  the  probe 
beam  with  a  second  pinhole  while  an  ac-coupled  pho¬ 
todiode  was  used-to  monitor  the  probe  intensity.  The 
Raman  gain  signal,  which  appears  as  a  10-nsec  pulse 
on  the  probe-l^er  intensity,  is  averaged  with  a  gated 
integrator  and  stored  on  a  disk  for  later  analysis. 

Density-broadening  coefficients  were  determined 
for  the  5(6),  S(8),  5(10),  and  5(12)  rotational  transi¬ 
tions  in  N2-  data  were  taken  at  three  tempera¬ 
tures  (295, 195,  and  80  K)  for  densities  spanning  the 
range  0.01-2  amagats.  The  results  of  these  measure¬ 
ments  are  shown  in  Fig.  1  and  listed  in  Table  1.  The 
values  for  the  density-broadening  coefficient  at  room 
temperature  are  in  agreement  with  ab  initio  theory^^ 
(solid  line  in  Fig.  1)  and  are  in  good  agreement  with 
recent  linewidth  measurements^^  made  using  sponta¬ 
neous  scattering.  Simultaneous  recording  of  the  Ra¬ 
man  signals  at  two  different  pressures  also  permitted 
the  observation  of  density  shifts  of  the  Raman  transi¬ 
tion  frequency.  We  have  determined  that  the  density 
shifts  are  less  than  150  MHz/amagat  for  all  rotational 
lines  and  temperatures  studied.  These  broadening 
and  shift  measurements  will  be  described  in  greater 
detail  in  a  future  publication.^^ 

The  ground-state  polarizability  anisotropy  is  de¬ 
fined  as  7  »  aj  -  aj.,  where  ( J. ,  ll)  indicate  the  perpen¬ 
dicular  and  parallel  components,  respectively,  of  the 
polarizability  a.  We  can  model  the  wavelength  de¬ 
pendence  of  7  by  using  an  empirical  formula^ 


where  v  is  the  pump-laser  frequency,  a,x,||  is  the  static 
or  dc  polarizability,  and  i'ix,ii  is  an  effective  intermedi¬ 
ate  state.  The  vix,!!  were  determined  first  by  fitting 
Eq.  (3)  to  the  results  of  recent  ab  initio  calculations^ 
for  a  I  and  ax.  The  a,;i  and  a,x  were  then  obtained  by 


fitting  Eq.  (3),  with  Vix,!!  fixed,  to  recent  experimental 
determinations®  ®  ®  of  ay  and  ax  to  obtain  a,!]  and  ajx- 
The  N2  parameters  derived  using  this  procedure  are  Wiy 
=  1.260  X  10®  cm->,  vix  =  1.323  X  10®  cm"!,  a,||  =  2.200 
X  10“®^  cm®,  and  a,x  =  1.507  X  10“®^  cm®. 

In  Fig.  2,  the  wavelength  dependence  of  7(1'),  from 
Eq.  (3),  is  shown  along  with  several  experimental  val¬ 
ues.  We  have  used  data  from  Refs.  6, 8,  and  9  only  to 
determine  the  a,x,i  because  these  measurements  have 
substantially  smaller  uncertainties  than  those  of  the 
other  references  cited  in  Fig.  2.  Equations  (2)  and  (3) 
yield  cross  sections  that  are  in  agreement  with  direct 
measurements.^®'^® 

In  N2,  the  J  dependence  of  the  polarizability  anisot¬ 
ropy,  7,  is  negligible^’'!®  for  J  ;S  20.  We  have  verified 
this  conclusion  by  determining  the  relative  cross  sec¬ 
tions  from  our  stimulated  Raman  line-shape  measure¬ 
ments.  Thus  we  have  used  Eq.  (3)  to  determine  the 
polarizability  anisotropy  for  all  transitions  studied  in 
this  Letter. 

The  reliability  of  Eq.  (3)  can  be  checked  by  using  the 
parameters  a,x,|  and  i'(X,||  to  calculate  the  index  of 
refraction  n;  the  result  can  then  be  compared  with 
experimental  data.  The  expression  forn  —  1  is 


n^-1 

n^-i-2 


4irN 

9 


where  N  is  the  density.  Over  the  wavelength  range 
250-600  nm,  expression  (4)  agrees  with  experimental 
data!®  better  than  0.5%,  verifying  a  high  accuracy 
for  y(v)  of  Eq.  (3). 

We  have  used  Eqs.  (l)-(3)  and  our  linewidth  data  to 
calculate  the  Raman  gain  coefficient  for  the  rotational 
transitions  as  a  function  of  temperature  in  the  limit 
where  density  broadening  is  the  dominant  contributor 
to  the  linewidth.  For  N2,  this  high-density  limit  cor- 


Table  I.  Temperature  Dependence  of  Rotational 
Raman  Linewidths,  Fractional  Population 
Differences,  and  Gain  Coefficients  for  S  Branch  in  Na" 


5-Brancb 

Transition 

S(J) 

Fractional 

Population 

Difference 

Broadening 

CoefHcient 

(MHz/amagat) 

Gain 

Coefficient 
(10-12  cm/W) 

T-295K 

6 

0.0282 

3570 

3.6 

8 

0.0337 

3280 

4.6 

10 

0.0335 

3070 

4.8 

12 

0.0289 

2870 

4.3 

r»  195  K 

6 

0.0485 

3070 

7.2 

8 

0.0491 

2860 

7.6 

10 

0.0399 

2660 

6.5 

12 

0.0271 

2340 

4.9 

T»80K 

6 

0.0897 

2530 

16.1 

8 

0.0452 

2120 

9.3 

10 

0.0156 

1950 

3.4 

12 

0.00378 

1690 

0.91 

“  (A,  •  568  nm). 
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Fig.  2.  Wavelength  dependence  of  the  N2  polarizability 
anisotropy.  Solid  line  is  given  by  Eq.  (3) ,  which  was  derived 
with  a  nt  to  the  data  of  Refs.  6, 8,  and  9.  Data  of  Refs.  2  and 
13  were  not  used  in  the  Ht  but  are  shown  for  comparison. 


Fig.  3.  Temperature  dependence  of  the  plane-wave, 
steady-state  R^an  gain  in  Nj  calculated  from  Eq.  (1)  for 
the  high-density  limit  (p  >  0.01  amagat).  The  pump  and 
Stokes  polarizations  are  linear  and  parallel. 


responds  to  densities  greater  than  0.01  amagat.  The 
results  are  summarized  in  Table  1,  while  the  tempera¬ 
ture  dependence  is  illustrated  in  Fig.  3.  These  results 
were  obtained  by  assuming  that  the  rotational 
linewidths  are  linear  with  temperature  and  the  polar¬ 
izability  anisotropy  is  constant  with  J. 

The  accuracy  of  the  g6un  coefflcients  reported  here 
is  limited  by  the  uncertainties  in  the  linewidths  and 
the  polarizability  anisotropy.  Our  linewidth  mea¬ 
surements  have  uncertainties  of  2%  (one  standard  de¬ 
viation)  as  given  by  the  linear  fits  used  to  determine 
the  broadening  coefficients.^^  The  imcertainty  in  the 
polarizability  anisotropy  was  estimated  to  be  1.5% 
from  the  differences  in  the  results  of  different  investi¬ 
gators.  Thus  the  accuracy  of  the  gain  coefficients 
reported  here  is  approximately  5%. 

We  can  compare  the  pure  N2  gain  coefficients  re¬ 
ported  here  with  the  air  coefHcients  reported  by  Hene- 
sian  et  al.,^  since  we  have  also  measured^^  the  foreign 
gas  broadening  coefflcients  of  O2.  After  correcting  for 
N2  and  O2  densities,  X„  and  the  wavelength  depen¬ 
dence  of  the  polarizability  anisotropy  [Eq.  (3)],  we 
obtain  a  N2  gain  coefficient  of  2.0  X  cmAv  for  1 


atm  (80/20%  mixtiu-e  of  N2/O2)  at  X*  =  1.07  ^m.  This 
is  25%  smaller  than  the  value  of  Henesian  et  al.^  but 
consistent  with  their  error  bar  of  a  few  tens  of  per¬ 
cents. 

In  conclusion,  we  have  made  linewidth  and  cross- 
section  measurements  that  will  facilitate  the  modeling 
of  stimulated  rotational  Raman  scattering  in  N2. 
Linewidths  are  linear  with  density  for  densities  from 
0.01  to  2.0  amagats  and  are  also  linear  with  tempera¬ 
ture  for  temperatures  from  80  to  300  K.  Within  the 
accuracy  (±5%)  of  our  relative  measurements,  the  po¬ 
larizability  anisotropy  was  found  to  be  independent  of 
J  for  5(2)  through  5(16).  The  new  linewidth  data 
presented  here  can  now  be  used  to  predict  accurately 
(dk5%)  the  rotational  Raman  gain  over  a  wide  range  of 
temperatures  and  densities  in  N2. 

This  research  was  supported  by  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  under  contract 
N00014-84-C-0256  through  the  U.S.  Office  of  Naval 
Research. 
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Model  of  the  rotational  Raman  gain  coefficients  for 
N2  in  the  atmosphere 

G.  C.  Herring  and  William  K.  Bischel 


A  model  for  stimulated  Raman  scattering  in  the  atmosphere  is  described  for  the  pure  rotational  transitions  in 
N2.  This  model  accounts  for  the  wavelength  dependence  of  the  polariaability  anisotropy,  altitude  and 
seasonal  temperature  variations  in  the  atmosphere,  and  the  O2  foreign-gas  density  broadening.  This 
information  is  used  to  calculate  the  steady-state  plane-wave  Raman  gain  profile  over  the  lower  100  km  of  the 
atmosphere.  Over  altitudes  of  0-40  km,  temperature  variations  produce  30%  changes  in  the  gain  coefficients 
of  1  km~'  cm^  MW~^  for  the  strongest  lines  at  Stokes  wavelengths  of  350  nm. 


I.  Mroduction 

Several  nonlinear  effects,  summarized  by  Zuev,^  are 
possible  when  a  high  intensity  laser  beam  propagates 
through  the  atmosphere.  Recent  work^  has  focused 
on  the  possibility  that  stimulated  rotational  Raman 
scattering  in  N}  may  have  the  lowest  threshold  intensi¬ 
ty.  Stimulated  Raman  scattering  then  limits  the  max¬ 
imum  laser  intensity  that  can  he  transmitted  through 
the  atmosphere  if  the  divergence  of  the  beam  is  re¬ 
quired  to  remain  unchanged.  To  establish  this  maxi¬ 
mum  intensity,  it  is  necessary  to  obtain  an  accurate 
descriptHm  of  the  steady-state  Raman  gain  coefficient 
in  the  atmosphere. 

This  paper  dssc^'bss  a  model  of  the  Raman  gain 
coefficient  for  the  *  2  S-branch  rotational  transitions 
in  the  atmosphere.  This  model  is  based  on  our  recent 
study^  of  the  temper  iture  dependence  of  the  rotation¬ 
al  Raman  gain.  In  the  present  atmospheric  work,  we 
use  the  same  temperature  dependence  as  described  in 
Ref.  3,  although  we  also  include  the  effect  of  foreign- 
gas  (O2)  broadening  on  the  N2  Raman  lines.  This 
temperature-dependent  Raman  gain  coefficient  is 
used  with  a  published  atmospheric  modeP  to  calculate 
the  atmospheric  N2  Raman  gain  coefficient.  The  im¬ 
portant  new  contribution  of  this  paper  is  a  tempera¬ 
ture-dependent  description  of  the  rotational  N2  Ra- 
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man  gain  coefficient  as  a  function  of  altitude  (0-100 
km)  in  the  atmosphere. 

The  modeling  of  the  propagation  of  a  high  intensity 
laser  pulse  through  a  Raman  medium,  such  as  the 
atmosphere,  requires  a  propagation  code  that  includes 
transient  effects,  off-axis  propagation  for  the  generat¬ 
ed  beams,  the  inclusion  of  all  possible  Stokes  and  anti- 
Stokes  orders  that  can  be  generated,  and  the  provision 
for  a  variable  pump  laser  bandwidth.  Codes  including 
some  of  these  effects  are  currently  being  developed  in 
several  laboratories  including  our  own.^  The  m^el  of 
the  steady-state  gain  coefficient  presented  here  is  an 
important  input  to  these  codes. 

1.  Model  Description 

A.  Atmospheric  Model 

In  our  atmospheric  gain  model,  the  Raman  gain 
calculation  uses  the  altitude-dependent  temperature, 
N2  density,  and  O2  density  as  input  data.  We  have 
taken  these  data  from  actual  atmospheric  measure¬ 
ments  at  Wallops  Island,  VA  (38°N),  as  published  by 
Banks  rmd  Kockarts.^  Slightly  different  results  are 
expected  if  the  standard  atmospheric  moder  is  used. 
Temperatures  and  densities  for  N2  and  62  were  taken 
from  Table  3.1  of  Ref.  5.  Densities  were  extrapolated 
for  the  first  15  km  (omitted  from  Table  3.1  of  Ref.  5)  by 
assuming  an  exponential  decay  with  increasing  alti¬ 
tude,  an  e~^  altitude  of  8  km,  and  relative  densities  of 
78  and  21%  for  N2  and  O2,  respectively.  The  average 
atmospheric  temperature  profiles  for  this  model  are 
shown  in  Fig.  1.  Maximum  variations  are  indicated  by 
horizontal  bars  for  summer  and  by  shading  for  winter. 
Temperatures  for  the  first  15  km  (also  omitted  from 
Table  3.1  of  Ref.  5)  were  determined  from  the  average 
of  the  summer  and  winter  values  in  Fig.  1. 
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Fig.  1.  Temperature  profiles  used  in  the  Raman  gain  calculations. 
Maximum  variations  are  indicated  by  horizontal  bars  for  summer 
and  by  shading  for  winter  (from  Ref.  5). 


B.  Raman  Gain  Mcxlel 

The  temperature  dependence  of  the  rotational  Ra> 
man  gain  is  described  in  Ref.  3  and  summarized  below. 
The  steady-state  plane-wave  Raman  gain  coefficient  is 
given  by^ 


g 


iff 

hwg  SQ 


(1) 


where  fidv)  is  the  area-normalized  Raman  line  shape 
and  5p  ^  Vp  —  vs  ~  *^i  where  Vp^  are  the  pump/Stokes 
laser  firequencies,  is  the  Raman  transition  frequen¬ 
cy,  and  As  is  the  Stokes  wavelength  in  the  gain  medium. 


The  spontaneous  scattering  cross  section  per  molecule, 
dff/dQ,  is  temperature  independent,  whereas  the  popu¬ 
lation  difference  AN  and  the  line  shape  f{5v)  are  the 
temperature-dependent  factors.  At  the  peak  (8v  =  0) 
of  a  Lorentzian-shaped  transition,  Eq.  (1)  reduces  to 
Eq.  (1)  of  Ref.  3.  Table  I  briefly  summarizes  the 
parameters  used  in  the  present  Raman  gain  calcula¬ 
tion. 

The  cross  section  is 

da  _  2  (J+1)  J+2)  j 

ail~  15\  c  /  (2J+1)  (2J  +  3) ' 

for  pump  and  Stokes  polarizations  that  are  linear  and 
parallel,  which  is  the  only  geometry  considered  in  this 
work.  The  polarizabilty  anisotropy  y  is  wavelength- 
dependent,  and  the  details  for  determining  the  wave¬ 
length  dependence  are  given  in  Ref.  3.  AH  the  results 
presented  in  this  study  are  for  a  Stokes  wavelength  As 
of  350  nm  and  a  polarizability  anisotropy  of  7.39  X 
10"2®  cm^.  We  neglect  the  rotational  quantiun  num¬ 
ber  J  dependence  of  y,  an  excellent  approximation  for 
N2,  as  discussed  in  Ref.  3. 

The  population  density  difference  is  given  by 

AN  =  N(J)-(^^)w').  (3) 

N{J)  is  the  number  density  of  molecules  with  rotation¬ 
al  quantum  number  J,  where  primed  and  unprimed  J 
denote  the  upper  and  lower  state,  respectively.  The 
population  difference  AN  was  calculated  from  a  Boltz- 
man  distribution  using  the  N2  rotational  constants  in 
Ref.  8. 

Only  the  high  density  limit  wa  .onsidered  for  the 
gain  calculations  of  Ref.  3,  and  thus  Lorc^tzian  func¬ 
tions  were  used  for  f(Sv).  For  the  present  atmospheric 
model,  we  use  a  more  general  Voigt  function  for  f(dv). 
In  the  most  general  case,  the  Galatry  (or  other  related 
line  shapes^)  line  shape  is  a  more  accurate  line  profile 
for  transitions  that  have  line  shapes  with  significant 
contributions  from  coUisional  narrowing.  However, 
our  use  of  Voigt  line  shapes  is  a  good  approximation  for 
rotational  N2  lines  in  the  atmosphere  for  two  reasons. 
First,  at  low  altitudes  (<20  km),  pressure  broadening 
dominates  over  Doppler  broadening  (5  MHz),  and  thus 
coUisional  narrowing  is  negUgible.  Second,  at  higher 
altitudes  where  Doppler  broadening  is  appreciable,  the 


TaM«  I.  In  «»•  8(8)  Qnin  Cniculallon  tor  A,  =  350  nm 


Altitude 

(km) 

Temperature 

(K) 

£kN 

N 

Voigt  A»twhm 
(MHz) 

f(iv  »  0) 

Gain  coefficient 
(cm^km-'MW-i) 

0 

285 

0.0349 

2960 

0.000215 

0.68 

10 

225 

0.0441 

784 

0.000813 

0.92 

20 

219 

0.0451 

180 

0.00354 

0.95 

30 

235 

0.0424 

39.8 

0.0164 

0.87 

40 

268 

0.0373 

10.5 

0.0595 

0.64 

50 

274 

0.0365 

5.7 

0.121 

0.36 

60 

253 

0.0396 

4.9 

0.166 

0.17 

70 

211 

0.0464 

4.5 

0.209 

0.067 

80 

197 

0.0489 

4.3 

0.216 

0.016 

90 

197 

0.0489 

4.3 

0.218 

0.0029 

100 

209 

0.0468 

4.4 

0.212 

0.00049 
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density  is  too  small  for  signiHcant  coUisional  narrow¬ 
ing.  Thus  pure  Voigt  line  shapes  were  used  for  all  our 
results. 

The  width  of  the  Lorentzian  contribution  to  the 
Voigt  was  obtained  by  adding  the  linewidths  due  to 
self-broadening  and  to  O2  foreign-gas  broadening. 
The  linewidth  measurements  used  here  are  summa¬ 
rized  in  Ref.  10.  'Self-broadening  coefficients  not  giv¬ 
en  in  Table  I  of  Ref.  10  were  determined  by  assuming  a 
linear  J  dependence.  Foreign-gas  broadening  coeffi¬ 
cients  for  temperatures  and  J  terms  not  given  in  Table 
V  of  Ref.  10  were  determined  by  assuming  that  the 
foreign-gas  values  are  85%  of  the  corresponding  self¬ 
broadening  values  in  Table  I  of  Ref.  10.  Finally,  all 
broadening  coefflcients  were  assumed  to  have  the  lin¬ 
ear  temperature  dependence  described  in  Ref.  10. 
These  three  assumptions  are  based  on  the  trends 
shown  by  the  data  of  Ref.  10. 

Oiur  c^culations  assume  that  the  laser  linewidth  is 
much  smaller  than  the  Raman  linewidth.  These  cal¬ 
culations  will  also  be  valid  if  the  laser  linewidth  is 
larger  than  the  Raman  linewidth,  provided  that  the 
Stokes  Held  is  phase  locked  to  the  piunp  laser  neld.^^~^^ 

III.  Results 

A.  Altitude  Dependence 

Figure  2(a)  shows  the  peak  (6v  =  0)  Raman  gain  coeffi¬ 
cient  of  S(8)  as  a  function  of  altitude  for  the  tempera¬ 
ture  proflle  described  in  Sec.  II.A.  The  two  quantities 
AN  and  A»t  that  combine  to  produce  this  altitude 
profile  are  also  shown  in  Figs.  2(b)  and  (c)  for  the  same 
temperature  profile  (Apt  is  the  FWHM  Raman 
linewidth.)  The  Raman  gain  is  about  constant  below 
40  km  and  drops  to  zero  above  40  km.  This  occurs 
because  the  peak  Raman  gain  depends  on  the  ratio 
ANIApt-  Below  40  km.  Figs.  2(b)  and  (c)  show  that 
AN  and  Apt  decrease  with  altitude  at  about  the  same 
rate,  yielding  an  approximate  constant  value  of  gain. 
Above  40  km.  Fig.  2(b)  shows  that  AN  continues  to 
decrease,  whereas  Fig.  2(c)  shows  that  Apt  becomes 
constant  with  altitude;  thus  the  Raman  gain  also  de¬ 
creases  with  altitude.  The  Raman  linewidth  is  con¬ 
stant  at  high  altitudes  because  the  coUisional  broaden- 
ing  becomes  negligible  compared  with  the 
approximate  constant  Doppler  broadening. 

Because  the  Raman  line  shape  is  pr^ominantly 
coUisionaUy  broadened  over  the  0-40-km  region,  a 
pure  Lorentzian  line  shape  should  be  a  good  approxi¬ 
mation  in  the  atmospheric  gain  calculation.  For  com¬ 
parison,  the  dashed  line  in  Fig.  2(a)  shows  the  gain  if 
the  line  shape  function  f{Ap)  in  Eq.  (1)  is  calculated 
from  a  Lorentzian  proflle  with  a  FWHM  linewidth 
given  by 

Ayj-  ■  (4) 

The  quantities  Apq  and  Ai/l  are  the  Doppler  and  coUi- 
sional  FWHM  linewidths,  respectively.  The  approxi¬ 
mation  is  in  good  agreement  with  the  fuU  Voigt  calcu¬ 
lation  except  for  a  20%  discrepancy  in  the  50-70-km 
region.  This  smaU  difference  is  important  for  the  gain 
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ALTITUDE  (km) 

Fig.2.  Raman  gain  vs  altitude.  The  solid  curves  show  (a)  the  peak, 
steady-state  plane-wave  Raman  gain  coefficient;  (h)  the  population 
difference  aN,  and  (c)  the  Raman  linewidth  (Voigt  FWHM)  as  a 
function  of  altitude  for  5(8)  in  Nj.  The  dashed  curve  in  (a)  shows 
the  gain  if  the  approximation,  *  Ai^  +  Av^  is  used  in  place  of 
the  Voigt  calculation  of  (c).  Ai>r.  Arc,  and  Aii.  are  the  total,  Gauss¬ 
ian,  and  Lorentzian  linewidths,  respectively. 


calculation  of  the  5(0)  through  5(4)  transitions.  The 
comparison  of  Fig.  2(a)  shows  that  a  simple  approxi¬ 
mation  of  the  N2  line  shape  is  sufficient  to  predict 
accurately  the  atmospheric  Raman  gain  for  the  most 
important  transitions.  However,  we  note  that  the  rest 
of  the  results  presented  here  were  obtained  with  com¬ 
plete  Voigt  calculations. 

The  peak  Raman  gain  coefficients  for  the  strongest 
Stokes  transitions  are  summarized  as  a  function  of 
altitude  in  Fig.  3.  The  difference  between  Figs.  3(a) 
and  (b)  is  that  the  vertical  scale  of  Fig.  3(b)  is  magni¬ 
fied  by  a  factor  of  3  to  better  illustrate  the  small  gain 
coefficients  of  5(0)  through  5(4).  The  number  indi¬ 
cated  next  to  each  curve  is  the  lower  rotational  quan¬ 
tum  number  J  for  that  particular  transition.  Peaks 
and  valle3rs  that  occur  at  an  altitude  of  15  km  are  due  to 
the  temperature  minimum  (200  K)  that  exists  at  this 
altitude.  Decreasing  temperatures  tend  to  transfer 
molecules  from  the  higher  rotational  levels  to  lower 
levels;  thus  AN  for  5(14)  and  5(16)  is  decreasing  while 
AN  for  5(6)  and  5(8)  is  increasing  at  the  temperature 
minimum.  In  Fig.  3,  5(8)  and  5(10)  have  the  largest 
gains  at  sea  level,  in  agreement  with  the  observation  of 
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Fig.  3.  Peak  Raman  gain  coeflicients  for  (a)  the  strongest  transitions  and  (b)  some  of  the  weaker  transitions  as  a  function  of  altitude. 


Henesian  et  al?  and  our  previous  calculations.^  The 
results  of  integrating  the  gain  coefficients  over  the 
depth  of  the  atmosphere  (from  sea  level  to  1(X)  km)  are 
tabulated  in  Table  II,  where  we  see  that  5(6),  5(8),  and 
5(10)  dominate  the  Stokes  conversion  process  in  the 
atmosphere.  The  odd-numbered  J  transitions  are 
only  half  as  strong  as  their  adjacent  even-numbered 
neighbors  and  are  omitted  from  this  investigation. 

Raman  gain  coefficients  for  frequencies  other  than 
those  at  the  peak  [Sp  «  0  in  Eq.  (1)]  of  the  Raman 
transitions  were  alro  calculated.  Figure  4(a)  shows 
the  5(8)  gain  coefficient  as  a  function  of  altitude  for 
detuning  values  of  0,  10,  100,  and  800  MHz.  The 
integrated  (from  0  to  100  km  in  altitude)  gain  coeffi- 


TaM*  S.  PMk  Raman  Oaki  Coafflcianla  (bilagralnd  (mm  Saa  Laval  to 
_ 180  km)  tof  Stofcaa  BotoMonal  TfanaHlona  at  Ma 


Transition  S(J) 

Gain  coefficient 
(cm*MW->) 

0 

2.16 

2 

13.8 

4 

29.1 

6 

39.3 

8 

43.4 

10 

39.8 

12 

32.6 

14 

22.6 

16 

15.2 
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Fig.  4.  S(8)  Raman  gain  coefficient  (a)  as  a  function  of  altitude  for 
detunings  of  0, 10, 100,  and  800  MHz  from  the  line  center  and  (b) 
integrated  (from  0  to  100  km  in  altitude)  gain  coefficient  as  a  func¬ 
tion  of  detuning. 

cient  for  5(8)  is  plotted  vs  detvming  in  Fig.  4(b).  The 
information  illustrated  in  Fig.  4  is  useful  for  determin¬ 
ing  the  effective  bandwidth  that  experiences  substan¬ 
tial  gain  over  the  entire  lower  100  km  of  the  atmo¬ 
sphere.  This  bandwidth  is  needed  to  calculate  the 
pump  laser  intensity  required  to  reach  threshold.^-^ 
We  estimate  this  effective  bandwidth  from  the  detun¬ 
ing  value,  where  the  generated  Stokes  intensity  is 


down  by  a  factor  of  2.  From  Fig.  4(b),  this  bandwidth 
is  2  MHz  and  is  to  be  compared  to  the  value  of  ~3  MHz 
that  we  estimate  from  Eq.  (16)  of  Ref.  4. 

B.  Temperature  Variatioiis 

All  the  results  of  Figs.  2-4  were  obtained  for  the 
speciHc  temperature  profile  of  Table  3.1  of  Ref.  5. 
The  effect  of  atmospheric  temperature  fluctuations  is 
summarized  in  the  four  curves  of  Fig.  5.  The  change  in 
Raman  gain  of  5(8),  due  to  seasonal  variations  of  tem¬ 
perature,  is  illustrated  with  four  different  gain  profiles. 
The  four  temperature  profiles  used  in  the  calculations 
of  Fig.  5  were  determined  from  Fig.  1.  The  curves 
labeled  maximum  are  the  gains  for  the  highest  tem¬ 
peratures  shown  for  that  season,  whereas  those  labeled 
minimum  are  the  gains  for  the  lowest  temperatures 
shown  for  that  season.  Figure  5  shows  that  irmTitnnm 
variations  in  the  Raman  gain  due  to  temperature 
changes  in  the  atmosphere  are  '^10%  at  all  altitudes. 

iV.  DIseusskxi 

This  is  the  first  study  to  consider  the  effect  of  tem¬ 
perature  variations  on  the  Raman  gain,  and  thus  it  is 
useful  to  compare  the  present  work  with  a  previous 
atmospheric  gain  calculation  that  neglects  ^e  tem¬ 
perature  effects.^  These  two  models  are  compared  in 
Fig.  6,  where  the  5(8)  gain  coefficient  is  plotted  vs 
altitude  for  each  model.  Table  in  lists  the  most  im¬ 
portant  parameters  associated  wtih  the  gain  coeffi¬ 
cient  for  both  models.  The  two  largest  differences  in 
these  two  models  are  the  temperature  distribution  of 
the  atmosphere  and  the  density  broadening  coeffi¬ 
cients  that  were  used.  Both  the  density  broadening 
coefficients  and  the  population  difference  change  by 
20-25%  for  5(8)  when  the  temperature  changes  from 


0  20  40  60  80  100 
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Fig.  5.  5(8)  peak  Raman  gain  variations  for  typical  variations  in  atmospheric  temperatures. 
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Fig.  6.  Comparison  of  the  temperature-independent  model  of  Ref. 
4  and  the  present  temperature-dependent  model  of  the  atmospheric 
Raman  gain  coefficient  for  S(8). 


Table  III.  Atmoapherle  Rotational  Raman  Gain  Coefficient  and  Aaaoclatad  Paramatara  lor  the  S(S)  Una  In  Nj  at  X(  ^  350  nm 


Room  temperature 
self-broadening 
coefficient 
(MHz/amagat) 

Polarizability 
anisotropy 
(10"“  cm"^) 

Sea  level 
temperature 
(K) 

Sea  level 
£lN 

N 

Sea  level 
gain  coefficient 
(cm2  lun-i  MW-‘) 

Threshold 

intensity 

(MW/cm2) 

This  work 

3280 

7.39 

285 

0.0349 

0.68 

0.85 

Ref.  2 

0.85“ 

Ref.  4 

2600 

7.47 

300 

0.035 

0.83“ 

0.86 

“  Gain  coefficients  were  scaled  to  350  nm  using  a  1/Xs  dependence. 


300  to  200  K.  The  combined  effect  of  the  atmospheric 
temperature  variation  on  Raman  gain  coefficient  is 
illustrated  by  the  difference  between  the  solid  and 
dashed  curves  of  Fig.  6.  The  broadening  coefficients 
used  in  the  present  work'®  are  ~20%  larger  than  those'® 
used  in  Ref.  4,  and  this  difference  is  the  major  contrib¬ 
utor  to  the  15%  difference  in  the  two  sea-level  (or  room 
temperature)  gain  coefficients  of  Fig.  6.  To  better 
visualize  the  effects  on  the  gain  due  only  to  temperture 
variations,  this  15%  sea-level  difference  should  be  ig¬ 
nored.  Thus  the  maTiTnnm  difference  between  the 
temperature  dependent  gain  and  the  temperature-in¬ 
dependent  gain  occurs  at  15  km,  where  the  tempera¬ 
ture-dependent  model  would  give  a  gain  coefficient 
30%  larger  than  the  temperature-independent  modeL 
Once  the  altitude  dependence  of  the  gain  coeffi¬ 
cients  is  known,  the  integrated  gain  can  be  computed 
and  used  to  estimate  the  threshold  intensity  for  appre¬ 
ciable  Raman  conversion.  We  calculate  the  threshold 
intensity  with  the  expression 

I,  ■  Is  ezp(-IpSg(z}dz),  (5) 

where  z  is  the  atmospheric  altitude,  g(z)  is  the  Raman 
gain  coefficient.  Ip  is  the  pump  intensity.  Is  is  the 
generated  Stokes  intensity,  and  In  is  the  spontaneous 
noise  intensity  from  the  Hrst  gain  length.^  To  com¬ 
pare  with  Ref.  4,  we  consider  an  atmospheric  path 
perpendicular  to  the  ground  and  define  threshold  to  be 
Is/Ip  *  0.01.  The  threshold  intensities  for  the  two  gain 
proffies  shown  in  Fig.  6  are  listed  in  the  last  column  of 
Table  III.  The  agreement  between  these  two  numbers 
illustrates  the  accuracy  of  the  temperature-indepen¬ 
dent  model  when  averaging  over  the  entire  depth  of  the 
atmosphere. 


Other  smaller  differences  between  our  work  and  that 
of  Ref.  4  include  the  wavelength  dependence  of  the 
polarizability  anisotropy  and  the  O2  foreign  broaden¬ 
ing.  We  have  used  the  wavelength  dependence  of  the 
rotational  transitions,®  and  the  results  of  Ref.  4  are 
based  on  the  wavelength  dependence  of  the  vibrational 
transitions.'®  Another  sm^  difference  is  that  we  have 
accounted  for  the  foreign-gas  broadening  contribution 
due  to  O2.  Each  of  these  effects  alters  the  gain  coeffi¬ 
cient  by  a  few  percent. 

The  uncertainties  in  the  Raman  gain  coefficients 
reported  for  this  study  depend  on  how  well  the  tem¬ 
perature  of  the  gain  medium  is  known.  If  the  tem¬ 
perature  is  accurately  known  (e.g.,  laboratory  condi¬ 
tions),  the  uncertainties  in  the  gain  coefficients  are 
limited  by  the  uncertainties  in  the  linewidths  and  the 
polarizability  anisotropy.  Both  of  these  uncertainties 
are  ~l-2%;  thus  the  gain  coefficient  uncertainties  are 
~5%.  In  the  atmosphere,  where  the  temperature  is 
not  well  known,  the  uncertainties  in  the  gain  coeffi¬ 
cients  are  dominated  by  temperature  uncertainties. 
Therefore,  gain  coefficient  uncertainties  are  ~10%,  as 
illustrated  in  Fig.  5. 

V.  Summary 

The  Raman  gain  of  the  Stokes  rotational  lines  in  N2 
has  been  modeled  as  a  function  of  altitude  for  the  lower 
100  km  of  the  atmosphere.  Calculations  were  com¬ 
pleted  for  the  even  J  transitions,  5(0)  through  5(16). 
The  strongest  lines,  5(6)  through  5(10),  have  maxi¬ 
mum  gain  coefficients  of  1  cm®  MW~'  km~'  at  an 
altitude  of  15-20  km  and  smoothly  decrease  to  zero 
over  the  20-80-km  region.  Effects  on  the  Raman  gain 
are  described  for  altitude  and  seasonal  temperature 
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fluctuations  in  the  atmosphere  and  detuning  from  the 
line  center  of  the  Raman  transition.  The  results  of 
this  work  provide  useful  input  for  simulation  codes 
that  model  the  stimulated  Raman  process  in  the  atmo¬ 
sphere. 

We  thank  David  L.  Huestis  for  stimulating  conver¬ 
sion  during  this  study.  This  work  was  supported  by 
the  Office  of  Naval  Research  under  contract  N00014- 
84-C-0256. 
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CW  Stimulated  Rayleigh-Brillouln  Spectroscopy 
of  Gases.  G.  C.  HERRING.  MARK  J.  DYER,  and  WILLIAM  K. 
BISCHEL,  SRI  International.  Menlo  Park.  CA  94025 --We 
report  the  results  of  CW  stimulated  Raylel^*Brlllioun 
spectroscopy  of  gases  at  568  nm.  These  are  the  first 
measurements  that  fully  resolve  the  details  of  the 
line  shapes  and  line  shifts  for  the  Brlllouln  doxiblet 
In  Xe,  SFg,  and  freons  with  a  resolution  limit  of  less 
than  10  MHz.  The  experiment  consists  of  two  single¬ 
frequency  CW  lasers,  a  Kr'*'  laser  at  568  nm  and  a  tun¬ 
able  ring  dye  laser  operating  at  568  nm,  that  counter- 
propagate  through  a  multipass  SBS  cell  that  refocuses 
the  beams  to  common  foci  for  30-60  passes.  The  lasers 
are  configured  In  a  counter-propagating  geometry  to 
probe  the  backward  SBS  gain  used  In  phase  conjugation 
applications.  The  ring  dye  laser  is  amplitude  modu¬ 
lated  at  0.1-3  MHz.  This  modulation  Is  transfered  to 
the  Kr'*'  probe  laser  by  SBS  gain  and  is  detected  using 
a  lock- In  amplifier.  The  resulting  line  shapes  are 
subsequently  fit  using  line  shape  theories  for 
Rayleigh-Brillouln  scattering.  The  data  presented 
here  is  important  for  applications  involving  phase - 
conjugation  using  stimulated  Brlllouln  scattering. 
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ABSTRACT 


We  have  developed  a  high  resolution  (10  MHz)  coherent  Brillouin  spectrometer  that,  for 
the  first  time,  clearly  resolves  the  Rayleigh  and  Brillouin  components  of  the  low  frequency 
scattering  spectrum  of  gases.  Two  narrow  band  cw  lasers  are  used  with  a  multi-pass  cell  to  obtain 
signal-tonoise  ratios  of  1(X).  We  have  used  this  apparatus  to  directly  measure  the  stimulated 
Brillouin  scattering  (SBS)  gain  coefficents,  linewidths,  and  line  shifts  for  Xe  and  SFe.  Our  results 
improve  the  accuracy  of  previous  measurements  and  illustrate  the  power  of  the  technique  for 
characterizing  other  potential  SBS  media  for  phase  conjugation  applications. 
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In  recent  years,  phase-conjugate  reflectors  have  generated  much  interest  for  a  variety  of 
applications.^  Many  nonlinear  techniques  can  be  used  as  phase-conjugators,  but  stimulated 
Brillouin  scattering  (SBS)  is  one  of  the  techniques  receiving  the  most  attention.  1*3  Optimum 
design  of  phaseconjugate  reflectors  will  require  high  resolution  experiments  to  accurately  measure 
the  lineshape  and  related  properties  of  Brillouin  scattenng.  In  the  past,^  the  frequency  resolution  of 
SBS  experiments  in  gases  has  been  limited  to  100  MHz  by  the  Fourier  limited  bandwiths  of  the 
pulsed  lasers  used  to  drive  the  Brillouin  interaction.  More  recently,  a  resolution  of  4  MHz  has 
been  demonstrated  using  two  cw  single-frequency  dye  lasers, ^  but  the  signalto-noise  ratio  (SNR) 
was  too  small  to  observe  the  Rayleigh  component.  In  the  present  work,  the  combination  of  large 
SNR  and  high  resolution  has  resulted  in  a  clear  separation  of  the  Rayleigh  and  Brillouin 
components. 

In  this  paper  we  report  the  developement  of  a  high  resolution  (10  MHz)  SBS  spectrometer 
to  study  backward  SBS  in  gases.  We  have  demonstrated  our  technique  on  Xe  and  SF6.  two  gases 
that  have  been  previously  studied  because  of  their  high  SBS  conversion  efficiency.  To  achieve  the 
highest  resolution,  two  single-frequency  cw  lasers  are  used,  while  a  multipass  cell  is  used  to 
increase  the  sensitivity.  The  high  resolution  of  our  system  has  allowed  the  direct  measurement  of 
the  Brillouin  linewidths  and  line  shifts,  while  our  high  SNR  and  the  relatively  well-known 
Gaussian  spatial  profiles  of  the  lasers  has  allowed  accurate  absolute  gain  coefficient  measurements 
for  both  Xe  and  SF6. 

The  experimental  set-up  is  shown  in  Fig.  1.  The  probe  is  a  Kr*’  -ion  laser  operating  at  568  nm 
and  uses  a  temperature  stabilized  oven  to  reduce  the  bandwidth  to  about  10  MHz.  The  pump  is  a 
tunable  ring  dye  laser  that  has  a  bandwidth  of  1  MHz.  Both  of  these  laser  beams  are  collimated  to 
the  same  size  (3  mm  diameter)  and  then  directed  into  a  multipass  cell  (30cm  focal  length),  where 
the  beams  are  overlapped  in  a  counterpropagating  geometry.  Both  the  pump  and  probe  beams  are 
passed  through  a  Faraday  isolator  to  prevent  laser  instabilities  due  to  optical  feedback.  The  pump 
beam  is  amplitude  modulated  at  100  kHz  with  an  acousto-optic  modulator  that  gives  85%  depth  of 


0-2 


modulation.  This  modulation  is  transferred  to  the  probe  through  the  SBS  process  when  the  angular 
frequency  difference  of  the  two  lasers  is  equal  to 

COB  =  2  (toi  V  n/c)  sin  (0/2), 

where  0  is  the  crossing  angle  of  the  pump  and  probe  beams,  coj  is  the  probe  laser  angular 
frequency,  v  is  the  speed  of  sound,  and  n  is  the  index  of  refraction.  A  small  percentage  of  the 
probe  is  then  picked  off  with  a  glass  slide  and  monitored  with  an  EG&G  SGD-444  photodiode, 
where  the  100  kHz  SBS  signal  is  detected  with  a  l(x:k-in  amplifier. 

An  example  of  the  SF6  Rayleigh-Brillouin  lineshape  is  given  in  Rg.  2.  The  dots  are  the  data  and 
the  solid  line  is  a  fit  to  the  data  based  on  a  Lorentzian  lineshape  mcxlel.  The  left-hand  and  right- 
hand  peaks  are  the  Brillouin  loss  and  gain  peaks  respectively.  The  central  structure  is  the  Rayleigh 
gain/loss  profile,  seen  in  SBS  as  the  derivitive  of  the  spontaneous  Brillouin  profile.  Our  best  SNR 
is  better  than  twice  that  shown  in  Fig.  2  and  is  about  100.  The  exact  stimulated  Rayleigh-Brillouin 
lineshape^  for  an  atomic  medium  is  actually  more  complicated  than  a  Lorentzian,  particularly  at  low 
pressures.  However  at  the  high  densities  used  in  most  of  this  study,  the  Lorentzian  approximation 
for  the  Brillouin  profiles  will  not  produce  significant  errors.  At  these  high  pressures  the  Rayleigh 
linewidth  is  substantially  narrowed  by  collisions,  thus  the  central  Rayleigh  profile  can  also  be 
approximated  as  the  derivative  of  a  Lorentzian. 

We  have  measured  the  pressure  dependence  of  the  Brillouin  linewidths  for  pressures  high 
enough  that  the  Brillouin  and  Rayleigh  components  remain  distinguishable.  These  measurements 
are  shown  in  Rg.  3  along  with  the  results  of  a  Damzen  et  al.^  who  have  recently  measured  the 
decay  time  of  the  acousdc  intensity  in  pulsed  SBS  experiments.  Hieir  results  yield  the  empirical 
relations. 


tB  (ns) — 0.65  X2  p 


(2) 


for  Xe  (p  =  pressure  in  atm  and  X  =  laser  wavelength  (air)  in  |im)  and 


XB-*  (sec*^)  =  ^  —  +  1.6  X  lO^  p  (3) 

pX 

for  SFg  (p  =  density  in  kg  m-3),  that  give  decay  times  as  a  function  of  pressure  or  density.  In 
reproducing  Eqn.  (3)  in  Fig.  3,  we  have  used  the  second  and  third  viiial  coefficients  for  the 
equation  of  state.8  We  have  converted  the  intensity  decay  times  of  Ref.  7  to  FWHM  linewidths 
using  Ato  =  1/t.  The  comparison  in  Fig.  3  shows  good  agreement  between  the  decay  time 
measurements  and  the  linewidth  measurements.  At  the  highest  pressures,  the  present  linewith 
measurements  are  about  10  MHz  larger  than  the  linewidths  predicted  by  Eq.  3. 

We  have  also  used  our  signal  levels  to  determine  the  absolute  steadystate  plane-wave  gain 
coefficients.  Assuming  that  the  pump  and  probe  beams  have  the  same  waist  and  confocal 
parameter,  the  plane-wave  gain  coefficient,  g,  is  related  to  the  multi-pass  focused  fractional  gain, 
G,  of  the  probe  by9 


G  =  (ngPA) 


(1-RN) 

(1-R) 


+  N  In  (R) 


(4) 


where  P  is  the  total  power  in  the  pump  beam,  N  is  the  number  of  passes  through  the  cell, 
and  R  is  the  mirror  reflectivity.  We  determine  G  by  directly  measuring  the  ratio  of  the 
Brillouin  signal  and  the  probe  power  levels.  This  ratio  is  then  corrected  for  pump  and  probe 
laser  power  changes.  The  results  of  our  measurements  are  given  in  Fig.  4  by  triangles  for 
Xe  and  squares  for  SF6. 

The  steady-state  SBS  gain  coefficient  can  be  calculated  fnxnl^ 
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(5) 


g  = _ _ 

—  n  c  V  p 
Tb 


where  Ye  is  the  electrostrictive  coefficient  given  by 

Ye  =  p  (9e/3p)T  =  j  (n^  +  1)  (0^+  2). 


(6) 


Using  the  literature^  values  of  the  density  dependence  of  the  index  of  refraction  to 
determine  Ye.  we  have  calculated  the  gain  for  Xe  and  SF6.  These  results  are  plotted  by  the 
dashed  lines  in  Fig.  4.  This  figure  shows  the  agreement,  within  10%,  of  absolute  gain 
coefficient  determinations  using  two  different  and  independent  methods.  The  agreement 
between  the  measurements  and  the  calculations  is  remaricable  given  the  uncertainties 
inherent  in  mode-matcing  and  overlapping  two  counter-propagating  laser  beams  for  50 
passes.  Thus  we  believe  the  present  technique  can  be  used  to  obtain  absolute  gain 
coefficients  at  the  10~  uncertainty  level  for  unknown  gases  and  liquids. 

Lastly,  we  have  measured  the  Brillouin  frequency  shifts  (center  of  Rayleigh  line  to  center 
of  Brillouin  line)  as  a  function  of  gas  density.  The  data  are  plotted  in  Fig.  5.  In  SFg,  the 
vertical  spread  in  the  data  at  20.5  atm  shows  typical  variations  for  three  successive  scans 
without  any  changes  in  the  experimental  apparatus,  while  the  vertical  spread  near  the  6  atm 
region  shows  the  variations  for  data  taken  on  different  days  with  die  frequency  scale 
recalibrated  between  the  two  different  days.  These  variations  give  an  indication  of  the  size 
of  the  systematic  errors  associated  with  our  measurements.  The  Brillouin  shift,  (Db.  is 
related  to  the  sound  velocity,  v,  by  Eq.  1  and  has  been  previously  studied^^  in  SFs  with 
spontaneous  Brillouin  experiments.  The  result  of  Ref.  12,  also  shown  in  Fig.  5,  consists 
of  an  empirical  relation  for  the  density  dependence  for  the  sound  velocity  and  is  consistant 
with  the  current  measurement.  In  contrast,  the  current  Xe  measurements  are  larger  than  the 
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standard  sound  velocity  equation  (see  caption)  at  20  atm.  We  believe  that  this  discrepancy 
in  Xe  is  real  since  the  difference  is  slightly  larger  than  the  error  bars  determined  by  the 
scatter  in  the  SF6  data. 

In  conclusion,  we  report  the  developement  of  a  high  resolution,  high  sensitivity 
stimulated  Brillouin  spectrometer.  The  system  has  been  demonstrated  with  gas  phase  Xe 
and  SFg,  in  which  we  have  measured  absolute  gain  coefficients  to  10%  accuracy  and 
lineshapes  (linewidths  and  line  shifts)  to  a  few  percent  accuracy.  These  types  of 
measurements  will  be  important  for  the  modeling  of  the  SBS  process.  In  addition,  with 
modest  improvements  in  the  signal-to-noise  ratio  and  the  frequency  resolution,  it  should  be 
possible  to  obtain  lineshapes  accurate  enough  to  yield  new  information  on  the  collision 
potentials'^  associated  with  atomic  collisions. 
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FIGURE  CAPTIONS 


1 .  Experimental  set-up  for  backward  coherent  Brillouin  scattering.  We  use  50  passes 
in  a  multipass  cell  that  has  a  single-pass  transmission  (window  losses  and  mirror 
reflectivity  combined)  of  R  -  0.985.  Laser  powers  into  the  cell  are  70  and  150  mW 
for  the  pump  (ring-dye)  and  probe  (ion)  respectively. 

2 .  Rayleigh-Brillouin  lineshape  for  5540  torr  of  SFg  at  room  temperature.  Dots  show 
data,  while  the  solid  line  shows  a  fit  that  consists  of  Lorentzians  for  the  Brillouin 
gain  and  loss  peaks  and  a  Lorentzian  derivative  for  the  central  Rayleigh  peak. 

3 .  Pressure  dependence  for  Brillouin  linewidths.  The  solid/dotted  lines  show 
linewidths  derived  from  the  decay  time  measurements  of  Ref.  7,  while  the  direct 
measurements  of  this  work  are  shown  for  Xe  (triangles)  and  SFg  (squares). 

4 .  Pressure  dependence  of  Brillouin  gain  coefficients  at  568  nm.  The  present  results 
for  Xe  (triangles)  and  SFg  (squares)  are  shown  along  with  calculations  of  Eq.  5 
(dashed  lines).  The  H2  Raman  gain  for  the  (3(1)  vibrational  transition  (solid  line)  is 
shown  for  comparison. 

5 .  Pressure  dependence  of  the  Brillouin  shifts.  Current  results(squares  for  SFe  and 
triangles  for  Xe)  are  compared  to  the  results  of  Ref.  12  (solid  curve)  for  SFe  and 
the  sound  velocity  equation  (dotted  curve),  v  (y  dp/dp)0^  S,  for  Xe,  where  y  is  the 
heat  capacity  ratio  C^Cy  -  5/3. 
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ABSTRACT 

The  construction  and  operation  of  a  long  pulse,  single-mode  Nd:YAG  laser  is 
described.  The  laser  is  configured  around  an  injection-seeded,  unstable  oscillator/single 
amplifier  and  produces  80  ns  pulses  at  10  Hz  with  pulse  energies  of  over  300  mJ.  When 
frequency  doubled,  the  infrared  is  converted  to  a  56  ns  pulse  at  532  nm  with  a  Fourier- 
limited  linewldth  of  8  MHz.  This  approach  brings  the  Nd:YAG  source  to  within  the  bandwidth 
domain  normally  dominated  by  visible  CW  lasers,  but  at  a  considerably  higher  power,  thus 
providing  an  attractive  alternative  for  applications  in  high  resolution  spectroscopy  and 
non-linear  studies. 


Li.iNTRQP.yCTION 

The  Importance  of  narrowband  laser  light  to  high  resolution  spectroscopic  studies  is 
well  recognized.  While  pulsed  sources  have  found  some  application,  the  use  of  linewidth- 
narrowed  and  stabilized  CW  ion  and  dye  lasers  has  long  been  the  convention  in  fulfilling 
experimental  requirements  for  narrow  bandwidths  of  better  than  10  MHz.  Unfortunately,  CW 
lasers  rarely  provide  line-narrowed  powers  greater  than  several  watts  without  extraordinary 
attention  to  stabilization  and  cooling  requirements.  For  kilowatt  or  greater  powers,  pulsed 
sources  appear  far  more  practical  and  efficient,  but  only  if  the  coherence  properties  of  CW 
lasers  can  be  emulated.  With  the  recent  development  of  laser-diode-pumped,  single-mode 
Nd:YAG  lasers^ and  the  demonstrations  of  injection  locking**’®  Fourier-limited  operation 
of  high-power  Nd;YAG  lasers  is  now  easily  performed  and  commercially  available.  Though  most 
attention  has  been  directed  toward  highest  power  pulse  durations  of  around  10  ns,  to  produce 
linewldths  of  less  than  10  MHz,  pulses  of  greater  than  44  ns  must  be  generated. 

The  production  of  long  single  mode  pulses  in  Nd:YAG  may  be  approached  in  several 
ways.  Certainly  the  most  versatile  source  would  result  from  chain  amplifying  an  acousto-  or 
electro-optically  shaped  single  mode  CW  laser  emission.  In  order  to  compete  with  the  power 
levels  of  other  available  continuous  lasers,  however,  this  method  can  become  quite  costly 
and  complex,  though  multipass  slab  amplifier  geometries  promise  to  make  this  a  more 
attractive  technique.®  Long  pulses  may  also  be  generated  in  stable  or  unstable  Q-switched 
Nd:YAG  oscillator  configurations.  Lengthening  the  pulse  build-up  time  or  forcing  the 
stimulated  photons  within  the  cavity  to  circulate  longer  before  saturating  the  rod  gain 
results  in  a  longer  ejected  pulse.  Furthermore,  injection  locking  these  longer  pulses  to  a 
single  mode  master  frequency  yields  narrowband  light  robust  enough  for  immediate  power 
amplification . 

We  have  constructed  a  stable  source  of  high  power,  Fourier  transform-limited,  80  ns 
wide  laser  pulses  using  an  injection-seeded  unstable  resonator  and  a  single  stage  of 
amplification,  which  provides  more  chan  300  mJ  of  output  energy  per  pulse.  The  design  is  a 
slightly  modified  version  of  an  earlier  system,^  but  Che  technique  may  certainly  be  applied 
to  virtually  any  Q-swltched  resonator. 

2.  SYSTEM  DESCRIPTION 


The  entire  system.  Illustrated  in  Fig.  1,  consists  of  five  principle  stages:  master 
oscillator,  slave  oscillator,  Isolation,  amplification  and  frequency  doubling.  The  master 
oscillator  is  a  Lightwave  Electronics  S-100  Injection  Seeding  System.  It  is  a  compact 
package  housing  a  well-stabilized,  laser-dlode-pumped ,  single-mode  Nd:YAG  ring  laser, 

Faraday  isolator  and  associated  electronics,  and  provides  almost  2.5  mW  of  CW  power  in  a 
1  mm  diameter  beam  at  its  output .  The  seeder  is  directly  mounted  Co  Che  Invar  resonator 
structure  within  a  Molectron  MY34-10  Nd:YAG  laser.  This  modified  chassis  contains  the  slave 
oscillator,  an  adaption  of  a  self-filtering  unstable  resonator®  design,  isolation  and 
amplification  stages,  and  frequency  doubling  and  harmonics  separation  optics. 
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Figure  1.  Component  schematic  for  single-mode  long  pulse  generation  and  amplification. 


2.1  ^ster/slave  oscillator 

The  master  frequency  from  the  seeder  is  polarization-coupled  into  the  self-filtering 
unstable  resonator  (SPUR),  where  It  effectively  fills  the  entire  volume  of  a  7inm  X  110mm 
Nd:YAG  rod  housed  in  a  Quantel  SF51107  dual  flashlamp  pump  chamber.  When  the  flashlamps  are 
fired,  the  rod  gain  is  swept  by  aproximately  .65  mW  from  the  seed  laser  until  Q-switch 
action  is  initiated.  At  the  peak  of  the  rod  gain,  the  Q-switch  is  opened,  momentarily 
closing  the  cavity  off  to  the  seed  and  permitting  the  master  frequency  to  be  resonantly 
amplified  within  the  slave.  A  piezo  stack  mounted  behind  the  end  reflector  M3  and 
controlled  by  Q-switch  buildup  time  (QSBOT)  reduction  electronics®  in  the  S-100  seeder 
provides  active  stablization  of  the  slave  cavity  length  in  order  to  support  the  master 
frequency.  Quarterwave  plates  on  either  side  of  the  rod  provide  the  necessary  suppression 
of  spatial  hole  burning^®  and  adjustable  outcoupling.  When  operated  at  high  flashlamp 
energies  near  the  limits  of  Q-swltch  hold-off,  this  resonator  produces  more  than  ZOO  mJ  in 
10  ns  on  a  single  axial  mode  and  with  near  diffraction-limited  beam  quality. 


Figure  2.  Plot  of  the  energy  per  pulse  as  a 
function  of  pulse  duration  for  the  injection 
seeded  SPUR  cavity. 


2.1.1  Long  pulse  generation 

Stretching  the  pulses  of  a  Q- 
switched  resonator  was  rather 
straightforward.  It  could  have  been 
accomplished  either  by  increasing  the 
length  of  the  oscillator,  introducing 
higher  intracavity  losses,  or  by 
changing  the  gain  of  the  active 
medium.  Lengthening  the  slave  cavity 
considerably  to  obtain  80  ns  pulse 
widths  was  quickly  dismissed  due  to 
considerations  of  mechanical  stability 
and  space  limitations.  We  avoided 
Introducing  additional  lossy  elements 
into  the  cavity  because  they  would 
also  reduce  the  seed  laser  power 
presented  to  the  rod  and  possibly 
degrade  spatial  performance.  The  two 
alternatives  available  to  us  to 
increase  the  losses,  Q-switch  control 
and  quarterwave  plate  rotation  to 
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diminish  feedback,  met  with  limited  success. 

Long  pulses  were  created,  but  they  proved  more 
difficult  to  control  and  injection  seed  than 
those  produced  simply  by  lowering  the  applied 
flashlamp  energy  and,  therefore,  the  gain  of  the 
rod.  With  some  repositioning  of  end  reflector 
M4  to  compensate  for  differences  in  rod  lensing, 
this  master/slave  combination  yielded  seeded 
pulses  continuously  variable  in  width  from  10  ns 
to  as  much  as  180  ns  at  FWHM,  with  relatively 
high  energies.  Fig.  2  shows  the  variation  of 
output  energy  with  pulse  length  and  Fig.  3 
illustrates  the  pulse  shapes  of  the  long  pulses. 


(b) 


Figure  3.  Oscilloscope  traces  of  the 
injection-locked  long  pulses  at  1 .064um 
for  (a)  80  ns  and  (b)  180  ns  at  FWHM. 


(a) 


Active  stabilization  of  the  slave  cavity 
for  80  ns  pulse  generation  is  only  slightly  more 
difficult  than  for  10  ns  pulse  operation. 

Because  the  80  ns  pulses  emerge  from  the  slave  a 
few  hundred  nanoseconds  past  the  Q-switch 
trigger,  instead  of  less  than  100  ns  for  the  10 
ns  pulse,  a  corresponding  delay  had  to  be 
introduced  into  the  trigger  of  the  seeder 
for  its  QSBUT  reduction  electronics  to  lock 
properly  (this  was  less  invasive  than  changing 
the  dynamic  range  of  the  reduction  circuit) . 
Also,  since  the  output  timing  of  the  slave  is 
processed  by  Che  reduction  electronics  via  a 
photodiode  monitoring  the  weak  back  reflection 
off  a  polarizer  in  the  seeder,  the  photodiode 
signal  had  to  be  boosted  for  the  lower  long 
pulse  intensities.  Taking  precautions  to 
isolate  the  slave  resonator  from  mechanical 
vibrations  was  the  last  concern  before  long  term 
stable  operation  was  possible.  The  final  output 
at  80  ns  was  10  mJ  at  1 .064  um  with  aprox.  *  10 
ns  timing  jitter,  attributed  mostly  to  pulse 
build-up  statistics,  and  less 
than  2  MHz  frequency  dither, 
determined  by  Che  piezo  movement 
and  the  slave  cavity's  free- 
spectral  range,  imposed  on  Che 
bandwidth  from  Che  active 
stabilization  by  the  seeder. 
Figure  4  presents  a  spectrum  of 
Che  frequency  doubled  output  at 
532  nm  taken  with  a  scanning 
confocal  interferometer. 


(b) 


L 


2.2  Isolation 

Of  paramount  importance  to 
the  stability  of  Che  seeded 
resonator  for  the  purpose  of 
power  amplification  is  the 
insertion  of  an  optical  isolator 
between  Che  oscillator  and 
amplifier.  This  decouples  Che 
two  stages,  preventing 
spontaneous  emission  in  Che 
amplifier  from  propagating  back 
into  the  oscillator,  during  Che 
buildup  of  Che  long  pulse,  and 
competing  with  Che  injected 
frequency.  Without  adequate 
isolation,  seeding  is  at  best 
erratic  and  active  locking 
impossible.  The  isolation  stage  consists  of  a  permanent  magnet  Faraday  rotator  manufactured 
by  Electro-Optics  Technologies  and  positioned  between  paired  thin  film  polarizers  and  a 
Glan-laser  polarizer.  With  AR-coatings  and  a  clear  aperture  of  8.5  mm,  Che  entire  isolator 


1.50  GHz 

Figure  4.  (a)  Temporal  profile  of  Che  second 

harmonic  and  (b)  optical  spectrum  of  Che  .532(im 
shows  good  agreement  with  Che  expected  Fourler- 
transform  bandwidth  of  8  Mhz  to  within  Che  resolution 
of  Che  scanning  confocal  interferometer. 
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assembly  furnishes  an  isolation  of  greater  than  40  dB  with  less  than  10%  total  insertion 
loss  . 


Tracing  the  optical  path  in  Fig.  1,  the  linearly  polarized  output  from  the  slave  cavity 
is  rotated  45"  from  vertical  by  the  halfwave  plate  HW1 ,  then  transmitted  through  the  thin 
film  polarizers  P4  and  P5  before  entering  the  Faraday  rotator,  where  the  polarization  is 
rotated  45“  back  to  the  vertical  and  passed  through  the  Glan-laser  polarizer  P6.  The  beam 
is  then  retroref lected  through  a  quarterwave  plate,  used  only  for  beam  attenuation  in  this 
scheme,  and  presented  to  the  amplification  stage  by  polarization  coupling  off  polarizer  P7. 
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Figure  5.  Single  and  double  pass  energy  plots  versus  applied  flashlamp  energy  for  seeded 
and  unseeded  operation. 


Figure  6.  Pulse  broadening  related  to 
sain  saturation  for  (a)  450  mJ  and  (b) 
340  mJ  double-pass  output  at  1 .064  pm. 


2.  _3 ,  FI  CAT  ION 

Amplified  energies  of  several  hundred 
millijoulos  were  initially  achieved  with  single¬ 
pass  amplification  through  a  Nd:YAG  pump  chamber 
identical  to  that  used  in  the  oscillator 
cavity.  With  less  than  10  mJ  into  the  rod, 
however,  flashlamp-pumping  the  rod  hard  enough 
to  reach  these  energies  severely  distorted  the 
spatial  quality  of  the  output  beam  and  left  the 
gain  unsaturated.  A  polarization-multiplexed, 
double-pass  geometry  was  then  adopted.  The 
vertically  polarized  beam  reflection  from 
polarizer  P7  passes  through  the  rod, 
retroreflects  off  mirrcr  M9  through  a 
quarterwave  plate,  where  it  undergoes  a  full 
halfwave  rotation,  and  then  passes  back  through 
the  amplifier.  The  amplified  beam  is  then 
transmitted  by  P7  and  down-collimated  for 
frequency  conversion.  For  this  amplification 
scheme,  the  Faraday  isolator  was  particularly 
valuable  because  it  removed  any  undesirable 
energy  rejected  from  P7  that  might  retrace  the 
optical  path  back  to  the  oscillator.  The 
results  of  the  single-  and  double-pass 
amplification  schemes  are  shown  in  Fig.  5 
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We  confirmed  chat  gain  saturation  was  occuring  above  300  mJ  for  double-pass 
amplification  but  observed  that  the  pulse  broadened  in  time.  As  seen  in  Fig.  6,  the  pulse 
width  increases  by  150%  for  a  flashlamp  energy  of  60  Joules  and  a  total  pulse  energy  of  540 
mJ . 

Fre^uencj  doublings 

Frequency  doubling  efficiencies  of  the  high  power  pulses  are  under  investigation  at 
this  time.  With  50%  down-coil imation  after  Che  amplifier,  we  expect  at  least  10%  conversion 
efficiency  for  the  unbroadened  80  ns  output.  Spectral  analysis  of  the  doubled  high  power 
will  determine  if  additional  frequency  components  are  created  by  the  amplification  and 
doubling  stages,  and  significantly  influence  the  bandwidth  integrity  of  the  unamplified 
emission. 

_3i_piscos^iqN 

We  found  that  pushing  this  technique  further  to  longer  pulses  appears  possible,  but 
owing  to  mechanical  and  thermal  instabilities  in  our  slave  resonator,  the  extremely  long 
pulse  build-up  times  and  subsequent  timing  jitter  of  Che  pulse  evolution,  the  seeder  could 

not  maintain  adequate  lock  for  more  than  a  few 
seconds  for  pulses  greater  than  100  ns.  By  modifying 
Che  dynamic  range  and  time  constant  of  the  time- 
reduction  circuit  within  the  seeder,  and  improving 
the  mechanical  and  thermal  stability  of  the  slave 
cavity,  we  hope  to  push  this  system  as  far  as  the 
master  oscillator  will  successfully  seed  the  SFUR 
cavity.  With  .65  mW  of  injected  power,  this  limit 
appears  to  be  just  under  200  ns  before  mode-beating 
oscillations  appear  on  the  pulse  waveform.  This 
restriction,  along  with  Che  relative  pulse  buildup 
times  of  different  pulse  lengths,  is  illustrated  by 
Fig.  7. 


Figure  7.  Oscilloscope  traces  superimposed  to  demonstrate  the  limits  of  this  technique  for 
.65  mW  of  injected  seed  power.  Long  pulses  may  be  sueded  from  10  ns  to  just  under  200  ns 
before  oscillations  appear  on  Che  waveform. 

4.  CONCLUSION 

We  have  demonstrated  that  long  pulse,  narrowband  operation  of  injection-seeded  Nd:YAG 
lasers  is  possible  by  varying  the  applied  flashlamp  energy  to  the  slave  oscillator.  With 
minor  modifications  and  Che  addition  of  optical  isolation  between  oscillator  and  amplifier, 
the  output  is  stable  and  may  be  power  amplified  directly  to  produce  several  hundred 
millijoules  of  1.064  ym  radiation  with  a  bandwidth  of  less  than  10  MHz. 
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VM34  LASER  FREQUENCY  CONVERSION  USING  TWO-PHOTON 
EMISSION  IN  BARIUM  VAPOR 


James  J.  Silva,  C.  C.  Liu,  Kenneth  Y.  Tang,  Western  Research  Corporation; 
David  L.  Huestis,  SRI  International 


ABSTRACT 

We  calculate  the  nonlinear  susceptibility  for  conversion  of  ultraviolet  krypton  fluoride 
excimer  laser  radiation  to  the  visible  using  stimulated  two-photon  emission  in  barium  vapor  and 
discuss  results  of  an  experiment  to  measure  this  conversion. 


Q-1 


FREQUENCY  HALVING 


—  USES  NONLINEAR  MEDIUM  (ATOMIC  BARIUM 
VAPOR)  TO  GENERATE  TWO  VISIBLE  PHOTONS  FROM 
ONE  UV  PHOTON. 


(0  -  CO  =  2co  +  As 

KrF  CO2  F.H.  Metastable 


—  SCHEME  IS  A  HYPER-RAMAN  PROCESS.  NO  PHASE 
MATCHING  IS  REQUIRED,  AND  HIGH  EFFICIENCY  AND 
BEAM  CLEANUP  CAN  BE  EXPECTED. 


—  LARGE  FREQUENCY  DOWNSHIFT  FROM  KrF 
(248nm)  TO  VISIBLE  (660nm)  DOES  NOT  RESULT  IN 
LARGE  ENERGY  WASTE  AS  IT  WOULD  IN  RAMAN 
SHIFTING  BECAUSE  TWO  VISIBLE  PHOTONS  COME 
OUT  FOR  ONE  UV  PHOTON  IN. 


—  HIGH  EFFICIENCY  VISIBLE  LASER  POSSIBLE 
STARTING  WITH  EFFICIENT  KrF  EXCIMER  LASER 


~  OUTPUT  AT  660nm  IS  IDEAL  BECAUSE  OF  HIGH 
OPTICS  DAMAGE  THRESHOLDS  AND  LOW 
ATMOSPHERIC  SCATTERING 
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FREQUENCY  HALVING:  BARIUM  ENERGY  LEVEL  SCHEME 


C87-6559 


FREQUENCY  HALVING  GAIN 

EQUATION 


—  FREQUENCY  HALVING  PROCESS  INVOLVES  A 
SEVENTH  ORDER  NONLINEAR  SUSCEPTIBILITY 


—  HIGH  SUSCEPTIBILITY  REQUIRES  NEARLY 
RESONANT  INTERMEDIATE  LEVELS 


—  KrF  AND  COg  LASERS  CAN  BE  TUNED  FOR 

STIMULATED  RAMAN  SCATTERING  TO*  BARIUM  LEVEL 
AT  39335  cm“i 


—  SUBSEQUENT  STIMULATED  TWO-PHOTON 
EMISSION  TO  LEVEL  AT  9034  cm'^  IS  DETUNED  BY 
ONLY  7cm-^  FROM  INTERMEDIATE  LEVEL  AT  24192 
cm"^ 


CM 


—  SUSCEPTIBILITY  IS  A  FUNCTION  OF  LASER 
POLARIZATION.  FOR  LINEAR  POLARIZATIONS,  KrF 
AND  COg  POLARIZATIONS  MUST  BE  PERPENDICULAR 

TO  EACH  OTHER. 


—  CALCULATION  OF  SUSCEPTIBILITY  INVOLVES 
SUMS  OVER  INTERMEDIATE  LEVELS  OF  PRODUCTS 
OF  EIGHT  TRANSITION  DIPOLE  MOMENTS  DIVIDED  BY 
PRODUCTS  OF  SEVEN  DETUNINGS  OR  LINEWIDTHS 

(7) 

^  ^  *^KrF  '  ^C02  *  ^F.H.^F.H.^C02  ‘  ^KrF^^F.H.^ 

=  (N/irA^)(2x10‘^^)  cm^erg^ 

WHERE  N=ATOMIC  NUMBER  DENSITY  (cm*^), 
A=DETUNING  FROM  S.R.S.  RESONANCE  (cm“^)  AND 
r=LASER  LINEWIDTH  (cm’^) 

—  GAIN  EQUATION  FOR  VISIBLE  SEED 

dlp  H/dz  =  (6x10-42 )  (NA’a2)I^Ico^(Ip  „  )2  (cm2/Watt)3{1/cm) 
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FREQUENCY  HALVING:  EXPERIMENTAL  SET-UP 


Q-6 


C87-6546 


FREQUENCY  HALVING 
EXPERIMENTAL  PARAMETERS 


—  TO  AVOID  OPTICAL  DAMAGE  TO  COATINGS.  FINAL 
EXPERIMENT  INCORPORATED  COLLIMATED  BEAMS 
CROSSING  AT  VERY  SMALL  ANGLES  IN  BARIUM 
VAPOR 

—  BARIUM  VAPOR  CONTAINED  IN  ALUMINA  HEAT 
PIPE  (TOP  TEMPERATURE  ABOUT  1200  CELSIUS) 
WITH  ARGON  BUFFER  GAS  (ABOUT  20  TORR)  AND 
BaFg  WINDOWS  TO  TRANSMIT  BOTH  KrF  (248nm)  AND 

CO2  (10.6]Li)  LASER  RADIATION 

—  CO2  LASER  AND  VISIBLE  DYE  LASER  POLARIZED 

VERTICALLY  AND  INCIDENT  FROM  ONE  END,  KrF 
LASER  POLARIZED  HORIZONTALLY  AND  INCIDENT 
FROM  THE  OTHER  END 

—  LASER  PARAMETERS  FOR  CENTRAL 
OVERLAPPING  50  cm  IN  BARIUM  VAPOR  (+/-  20%) 

CO2:  1000  mjoule  on  0.5  cm^  in  about  1|iS  with  a  peak 

intensity  in  first  100  ns  of  3.3  MW/cm^.  Single  line 
bandwidth  of  0.07  cm"^ 

KrF:  55  mjoule  on  0.2  cm^  in  16.2  ns,  with  75% 
narrowband  in  correct  polarization  for 
12.7  MW/cm^.  Bandwidth  0.05  cm■^  tunable  in 
mode  hopping  increments  of  about  0.2  cm’^ 

Visible:  0.52  mjoule  on  0.005  cm^  in  8.7  ns  for 

11.5  MW/cm^.  Bandwidth  0.1  cm■^  tunable  In 
0.05  cm"^  Increments 
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EXPERIMENTAL  PROCEDURE 


—  BEAM  OVERLAP  ASSURED  BY  CHECKING  BEAM 
POSITIONS  AT  BOTH  ENDS  OF  BARIUM  CELL 


—  KrF  LASER  TUNED  TO  GIVE  MAXIMUM  COg 

STIMULATED  RAMAN  SCATTERING  TO  LEVEL  AT 
39335  cm-^  DETERMINED  BY  MAXIMIZING 
SUBSEQUENT  470  nm  AMPLIFIED  SPONTANEOUS 
EMISSION  TO  LEVEL  AT  18060  cm-^ 


—  RELATIVE  TIMING  OF  LASER  PULSES  SET  TO 
MAXIMIZE  A.S.E. 


—  TIMING  OF  VISIBLE  L^SER  SET  SO  THAT 
BEGINNING  OF  VISIBLE  LASER  PULSE  ENTERS  CELL 
WHEN  BEGINNING  OF  KrF  LASER  PULES  EXITS  CELL 
(RELATIVE  JITTER  ABOUT  1  ns) 


—  KrF  LASER  DETUNED  SOME  AMOUNT  FROM  PEAK 
OF  S.R.S.  TRANSITION 


—  VISIBLE  LASER  SCANNED  IN  WAVELENGTH  WHILE 
•out/*in  monitored  on  BOXCAR  INTEGRATOR  AND 
PLOTTED  AS  A  FUNCTION  OF  TIME 
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EXPECTED  FREQUENCY  HALVING 

GAIN 


—  FOR  TEMPERATURE  OF  1150  CELSIUS,  BARIUM 
VAPOR  PARTIAL  PRESSURE  IS  ABOUT  28  TORR, 
CORRESPONDING  TO  N=2x10'Vcm3 


—  =  1 2.7  MW/cm2,  Iqq^  =  3.3  MW/cm2, 

VISIBLE  SEED  LASER  Ip  =  1 1 .5  MW/cm^ 


—  WIDTH  OF  KrF  -  COg  STIMULATED  RAMAN 

TRANSITION  WAS  MEASURED  TO  BE  ABOUT  1 .5  cm'' 
(AT  THIS  TEMPERATURE  WITH  THESE  LASER 
INTENSITIES).  A  DETUNING  OF  1  cm'’  LEADS  TO  A 

QUADRATURE  A  OF  1.8  cm'’.  LASER  LINEWIDTH  r=0.1 
crtrL 
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—  GAIN  EQUATION  FOR  VISIBLE  SEED 

dip  H/dz  =  {exIO-'ts )  {NA'A2)l^plco^(lp  „  )2  (cm2/Watt)3(1/cm) 

LINEAR  GAIN  (dip  ^/dzl/lp  „  =  .18%/cm  GIVING 

MEASURABLE  GAIN  OF  9%  FOR  50  cm  INTERACTION 
LENGTH. 


—  EQUATION  dl/dz  =  a|2  HAS  SOLUTION 
l(z)  =  (1/l(z=0)  -  az)-^  LEADING  TO  COMPLETE 
CONVERISON  AT  z  =  560  cm  (IGNORING  DEPLETION 
OF  KrF  PUMP  BEAM) 


—  NO  EXPERIMENTS  (2%  SIGNAL  TO  NOISE)  SHOWED 
GAIN 
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TYPICAL  FREQUENCY  HALVING  GAIN  MEASUREMENT  SPECTRUM 
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DISCUSSION  OF  EXPERIMENT 


—  BARIUM  CHOSEN  AS  FREQUENCY  HALVING 
CANDIDATE  ATOM  BECAUSE  OF  NEAR 
INTERMEDIATE  RESONANCES  AND  ADEQUATE 
VAPOR  PRESSURE  AT  ACCESSIBLE  TEMPERATURES 


—  FAILURE  TO  MEASURE  GAIN  ON  VISIBLE  SEED 
LASER  PROBABLY  DUE  TO  TWO  FACTORS:  LARGER 
THAN  EXPECTED  TRANSITION  LINEWIDTHS  AND 
LOWER  THAN  EXPECTED  DIPOLE  TRANSITION 
PROBABILITIES 


—  MEASURED  STIMULATED  RAMAN  SCATTERING 
LINEWIDTH  OF  1.5  cm-*  WAS  A  FACTOR  OF  THREE 
LARGER  THAN  EXPECTED.  LARGE  WIDTH  OF 
RYDBERG  TRANSITION  MAY  BE  DUE  TO  STARK 
EFFECT  FROM  LARGE  DENSITY  OF  IONS  AND 
ELECTRONS  CREATED  FROM  AUTOIONIZATION  OF 
BARIUM  ATOMS  ABSORBING  TWO  KrF  PHOTONS  (AS 
EVIDENCED  BY  A.S.E.  AT  650nm  FROM  Ba'*’) 


—  FREQUENCY  HALVING  LINEWIDTH  T  CALCULATED 

AS  LASER  LINEWIDTH  QF  0.1  cm’*  MAY  BE  LARGER 
DUE  TO  ATOMIC  EFFECTS 
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—  WHERE  POSSIBLE.  TRANSITION  DIPOLE  MATRIX 
ELEMENTS  WERE  CALCULATED  FROM  MEASURED 
EINSTEIN  A  COEFFICIENTS.  FOR  RYDBERG 
TRANSITIONS,  CALCULATIONS  WERE  BASED  ON  THE 
COULOMB  APPROXIMATION. 


—  SUBSEQUENT  EXPERIMENTAL  MEASUREMENTS 
INDICATE  THAT  ONE  CRITICAL  TRANSITION  MOMENT 
FROM  THE  LEVEL  AT  24192  cm‘^  TO  THE  LEVEL  AT 
39335  cm-l  MAY  BE  SUBSTANTIALLY  LOWER  THAN 
CALCULATED.  IT  WAS  NOT  SEEN  IN  AN  ABSORPTION 
EXPERIMENT  THAT  DID  MEASURE  TRANSITIONS  TO 
OTHER  NEARBY  RYDBERG  LEVELS.  THE 
CALCULATED  TRANSITION  MOMENT  WAS  BASED  IN 
PART  ON  A  PUBLISHED  M.Q.D.T.  CALCULATION 

SHOWING  MIXING  BETWEEN  THE  6s9d  ^02  AND  THE 
5d8s3D2  LEVELS. 
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